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ABSTRACT 


The  principal  objective  of  this  study  vu  to  perform  bade  research  Investigations  into  an 
innovative  power  conversion  concept  for  trans- atmospheric,  beamed-energy  propulsion:  a  new  class 
of  External  Surface  Impulse  (ESIl  thrusters.  This  advanced  thruster  principle  could  be  used  for 
atmospheric  VTOL,  high  acceleration,  and  lateral  flight  (e.g.,  short-term  ‘cruise”)  propulsion  of 
SIngle-Stage-To-Orbit  (SSTO)  beam-powered  shuttleeraft  of  the  next  century. 

Three  dams  of  ESI  thrusters  were  initially  examined:  l)simple  thermal,  2)electrostatic, 
and  3)electro magnetic.  Beam  power  wavelengths  from  10cm  (microwave)  to  0.3um  (laser)  were 
considered.  The  subsequent  research  effort  concentrated  on  the  simple  thermal  repetitively-pulsed 
ESI  thrusters,  energized  with  laser  power  and  using  air  as  the  working  fluid.  Laser  frequencies 
were  selected  because  of  the  relative  wealth  of  experimental  data  and  theoretical  research  on  "laser 
impulse  coupling”  which  exists  in  the  literature. 

The  first  year  analytical  effort  has  proven  conclusively  that  such  an  engine  can  deliver  high 
levels  of  thrust-to-beam-  power  at  liftoff  (e.g.,  at  least  an  order-of-magnitude  greater  than  beam- 
powered  hydrogen-fueled  rockets),  with  infinite  specific  impulse  (decreased  only,  perhaps,  by  ablation 
of  the  thruster  surface).  Later  along  an  orbital  trajectory,  the  primary  propulsion  function  would 
transition  to  other  modes;  upon  leaving  the  atmosphere,  the  SSTO  vehicle  would  continue  in  a  pure 
rocket  mode  with  a  specific  impulse  (/«p)  of  1000  seconds  or  more.  The  prindpal  motivation  behind 
this  variety  of  "combined-  cycle"  tranaatmospheric  propulsion  is  the  potential  for  "more  payload  to 
orbit”  for  a  given  laser  power  level  -  than  given  by  beam-powered  rockets  alone. 
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CHAPTER  I 


INTRODUCTION 

The  fint  yew  effort  began  with  an  initial  investigation  of  several  alternative  approaches 
to  External  Surface  Impulse  (ESI)  thrusters:  1)  simple  thermal,  2)  electoetatic,  and  3)  electro¬ 
magnetic  varieties.  The  primary  focus  of  the  subsequent  investigation  wu  on  the  simple  thermal 
class  of  thrusters:  specifically,  the  External  Radiation-Heated  (ERH)  thruster.  For  the  purpose  of 
comparison,  the  latter  two  thruster  types  will  be  described  first. 

A)  EUetoeUtk  ESI  Thruster 

The  electrostatic  ESI  thruster  concept  (shown  in  Fig.  1-1)  utilises  a  remote  projected  mi¬ 
crowave  or  laser  power  beam  to  accelerate  and  eject  relativistic  streams  of  electrons  to  some  remote 
location  away  from  the  thruster  surface.  This  leaves  the  vehicle  exterior  with  a  large  resultant 
positive  charge.  The  ejected  electrons  preferentially  attach  to  previously  neutral  oxygen  or  water 
molecules  (in  the  ambient  air)  thereby  creating  a  negatively  charged  cloud.  Impulse  is  then  gener¬ 
ated  by  electrostatic  forces  occurring  between  the  charged  cloud  and  vehicle  surface.  As  the  vehicle 
moves  through  the  negative  ion  cloud,  its  surface  charge  is  neutralised;  then  the  entire  process  is 
repeated. 

B)  Electromagnetic  ESI  Thruster 

The  electromagnetic  type  of  ESI  thruster  uses  microwave  or  laser  energy  to  fint  ionise  the 
ambient  air  working  fluid  near  the  thruster  surface  (e.g.,  see  Fig.  1-2).  Once  ionised,  the  driver  loop 
(an  electromagnet)  is  pulsed  and  an  electric  field  is  induced  in  the  gas.  The  gas  is  then  repelled 
from  the  thruster  surface  by  means  of  electromagnetic  forces,  (i.e.,  F  =  J  x  B),  thereby  generating 
thrust.  Presumably,  a  powerful  and  lightweight  on-  board  eleetrie  power  supply  would  be  required, 
one  which  could  be  recharged  with  beamed  energy.  Refer  to  Fig.  1-3  which  shows  the  electrical 
circuit  model  for  one  exemplary  electromagnetic  ESI  thruster. 

Both  electrostatic  and  electromagnetic  ESI  thrusters  are  technologically  complex  and  require 
additional  onboard  power  conversion  equipment  -  with  which  simple  thermal  engines  would  not 
be  burdened.  For  such  complicated  thrusters,  it  is  presently  unclear  as  to  which  flight  velocities 
and  altitudes  would  yield  superior  performance,  and  how  they  could  be  integrated  into  energetic 
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Figure  1-1  Puleed  Electroetatlc  ESI  Thruater 
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Figure  1-2  Pulsed  Electromagnetic  ESI  Thruster 


combined-cycle  engines. 

On  the  other  hud,  the  reults  of  the  present  etudy  have  shown  that  simple  thermal  ESI 
thrusters  are  ideally  suited  for  the  initial  liftoff  boost  period  (e.g.  Mach  0  to  3)  of  advuced  air- 
breathing  shuttlecraft,  are  easily  integrated  into  a  variety  of  combined-cycle  engine  concepts,  ud 
require  minimal  energy  conversion  (on-board)  equipment.  For  these  reasons,  the  engine  holds  much 
promise  for  use  with  beam-  powered  transatmospheric  vehicles  of  more-or-less  ‘conventional*  design 
(i.e.,  heavier  than-air,  ud  relatively  compact  in  size). 

C)  Simple  Thermal  ESI  Thruster 

The  External  Radiation-Heated  (ERH)  thruster  concept  is  pictured  in  Fig.  H.  The  reac¬ 
tion  surface  is  integrated  with  the  underside  exterior  structure  of  the  vehicle.  The  actual  thruster 
geometry  can  vary  greatly  depending  upon  the  desired  vehicle  configuration.  ERH  surface  shapes 
considered  here  were  circular  fnutrums  ud  fiat  circular  plates,  which  are  appropriate  for  vehiclea 
of  radial  symmetry  (e.g.,  the  Apollo  Command  Module).  Other  shapes  would  be  more  suitable  for 
delta- winged  abuttlecnft,  like  the  space  shuttle  orbiter. 

To  create  thrust,  high  intensity  electromagnetic  radiation  (either  microwave  or  laser  fre¬ 
quencies)  must  be  projected  radially  inward  from  the  vehicle’s  perimeter,  ud  brought  to  focus  near 
the  geometric  center  of  the  ERH  thruster  surface.  This  radiation  ignites  and  sustains  %  number  of 
air-plasma  “waves”  which  propagate  radially  outward  (ud  parallel  to)  the  thruster  surface  as  shown 
in  Fig.  1*4.  These  laaarsustained  plasmas  produce  blast  waves  which  subsequently  expand  over  the 
ERH  surface  -  generating  impulsive  thrust. 

The  process  of  thrust  production  is  described  more  thoroughly  in  the  following  sequence: 
i)  Depending  on  the  frequency  and  intensity  of  focused  radiation,  either  a  Laser  Supported  Det¬ 
onation  (LSD)  or  Combustion  (LSC)  wave  is  initiated  by  the  mechanism  of  inverse  Brems- 
strahiung.  In  these  radiation-sustained  plasma  waves,  electromagnetic  energy  (e.g.4aser 
light)  is  converted  to  thermal  energy  -  resulting  in  a  high  temperature  zone  with  elevated 
pressure.  In  the  case  of  u  LSD  wave,  a  very  high  pressure  blast  wave  is  formed.  However  if 
u  LSC  wave  is  used,  the  resultut  blast  wave  pressure  is  several  orders  of  magnitude  lower, 
but  furthermore  the  plasma  temperature  is  much  higher. 

U)  As  the  blast  wave  decays  to  the  local  ambient  pressure,  a  reaction  impulse  is  transmitted 
to  the  thruster  surface.  Some  portion  of  this  impulse  is  converted  to  usable  upward  thrust, 
depending  on  the  inclination  of  the  ERH  thruster  surface  to  the  flight  axis, 
iii)  Once  the  blast  wave  has  decayed,  the  region  adjacent  to  the  thruster  surface  still  contains 
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Figure  1-4  Simple  Thermal  ESI  Thrueter 
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hot  |bm.  That*  “spent*  |bm  moat  be  removed  tad  replaced  with  unprocessed  ambient 
air  baton  a  new  thrust  cycle  can  be  started.  Expressed  another  way,  the  surface  most 
be  “refreshed".  This  process  is  equivalent  to  the  had  rejection  part  of  a  conventional 
thermodynamic  cycle.  Unlike  most  current  propulsive  engines  which  utilise  “pulaejet*  heat 
addition  in  a  Brayton  cycle,  the  ERH  thruster  is  closest  to  the  “constant  volume”  class  of 
enfines  (e.g.,  Bke  the  V-l  “Buss  Bomb*  used  by  Germany  afainat  Enfiand  in  World  War 
D). 

Note  that  the  ERH  thruster  is  a  repetitively-pulsed  engine.  Its  Pulse  Repetition  Frequency 
(PRF)  is  loverned  by  the  blast  wave  decay  time  and  the  thruster  surface  ssSaA  tiffifi-  Therefore, 
the  total  cycle  time  is  dependent  on  the  ambient  air  preawre  and  density  in  which  the  ERH  thnater 
operates.  This  implies  that  the  thruster  performance  varies  with  both  altitude  and  ffifht  velocity. 
The  maximum  engine  PRF  varies  throughout  the  entire  Sight  envelope;  a  typical  cycle  time  is  on 
the  order  of  a  millisecond. 

Although  not  specifically  examined  in  the  first  year,  it  is  important  to  consider  the  use 
of  microwave  frequencies  with  the  "simple  thermal  dam*  of  ESI  thrusters.  Less  literature  exists 
on  microwave-sustained  plasmas,  but  sufficient  understanding  has  now  been  accumulated  to  begin 
research  on  microwave  beam  energised  ESI  therms!  engines. 

These  long  wavelengths  (microwave)are  critically  Important  for  the  future  beamed-powered 
SSTO  ahuttleenft.  Air  Force  scenarios  generally  dictate  all-weather  operation,  and  only  microwave 
beams  can  penetrate  doude/rain/fog  (all  forme  of  precipitation)  without  significant  attenuation. 
By  contrast,  laser  beams  must  “burn-through”  clouds,  exploding  and  evaporating  water  droplets 
-  a  brute  force  solution  that  Is  highly  wasteful  of  beamed  power.  A  prefenrred  solution  is  for  the 
shuttlecraft  to  perform  a  “pop-up*  maneuver,  wherein  It  climbs  vertically  through  the  doud  layers 
(e.g.,  up  to  30-40,000  ft)  under  a  low-power  microwave  beam,  then  rotates  to  engage  a  high  power 
laser  beam  link  for  the  final  boost.  It  has  become  apparent  that  both  modes  will  be  important  for 
future  SSTO  beam-powered  vehicles. 

D)  Report  Organisation 

Chapter  II  covers  the  reference  vehicle/platform  configuration  assumed  in  the  present  re¬ 
search  effort.  Since  the  method  of  beating  the  engine  working  fluid  (air)  requires  the  use  of  optics 
and  an  air  refresh  scheme  (lateral  and/or  axial  modee),  it  is  clearly  platform  dependent. 

Chapter  III  reviews  the  basic  phenomenology  of  inverse  bremastrahhug  absorption  of  laaor 
energy  into  tbe  air  working  fluid. 
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Chapters  IV  and  V  introduce  two  modda  for  LSD- wave  heated  ERH  thrusters.  These  models 
an  of  tbo^latkl  variety,  in  that  they  do  not  directly  take  into  account  a  particular  “refresh*  scheme 
for  convicting  a«w  atr  Into  tbo  active  thruster  volume.  Chapter  V  ako  analyse*  the  “dynamic*  per 
fonnance  of  a  conical-  fruatrum  thruater  surface*. 

Chapten  VI  and  VH  describe  anniyaea  of  an  LSC  wave  heated  EM?  thruater.  Chapter  VI 
review*  the  impulse  feneration  model,  and  aufgeata  additional  refinements  which  an  neceaaaiy.  A 
radiation  heat  tranafer  model  waa  aiao  created,  and  uaed  to  aaaeaa  the  aeverity  of  the  surface  heatinf 
problem.  The  temperature  hiatory  reeulta  an  pnaented  in  Chapt.  VO. 

Chapter  VIS  provider  an  analyaia  of  the  far  field  eound  level*  produced  by  the  reference  ERH 
thruater  configuration.  In  Chapt.  VIS,  this  analytical  model  is  used  to  simulate  the  eound  natter™ 
(La,  frequency,  amplitude,  doppler  shift,  etc.)  produced  by  an  ERH  thruster  as  it  accelerates  aloof 
representative  boost  trajectories  -  as  perceived  by  obeerven  on  the  Earth’s  surface. 

Finally,  Chapt.  X  fives  a  summary  of  principal  conclusions  reached  in  the  first  year  fundiiif 
segment. 
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n  V LIGHT  PLATFORM  REFERENCE  CONFIGURATION 


The  primary  task  of  this  "bask  research*  imrotigation  was  to  develop  tad  exercise  analytical 
models  for  eatfmating  the  performance  of  advanced  beamed-energy  propulsion  concepts.  However, 
s  reforeace  flight  platform  geometry  was  adopted  in  order  to  introduce  the  combined  requirements 
for  primary  and  secondary  optica,  impulse  feneration,  surface  contours,  and  "refresh”  modes  for  the 
ERH  thruster.  In  addition,  the  pursuit  of  a  reference  configuration  helped  to  clarify  how  the  ERH 
thruster  concept  mifht  be  incorporated  into  variable-cycle  engines  optimised  for  acceleration  along 
an  Earth- to-orbit  trajectory,  wherein  rockets  an  used  for  final  orbit  injection. 

These  innovative  airbreathing  engines  must  incorporate  the  functions  provided  by  conven¬ 
tional  inlet  and  nossle  (exhaust)  components.  It  is  interesting  to  note  that  for  hypersonic  spends, 
inbreathing  inlets  can  become  very  large  -  large  enough  to  house  the  entire  payload.  Hence,  the 
following  discumion  uses  the  terms  “vehicle*  and  “engine*  almost  interchangeably. 

A)  Radar— ie  Platform 


The  reference  eogfoe/vehkie  shape  is  that  of  a  five  meter  diameter  axi -symmetric  body 
largely  composed  of  conic  (rostrums,  as  illustrated  in  Fig.  H-l.  The  body  can  be  divided  into  three 
bask  components:  forebody,  cowl  sad  afterbody. 

The  forebody  consists  of  an  external  isen tropic  “spike*  difueer  with  a  conical  tip  of  15 
degrees  (semi-vertex  angle)  and  a  35  degree  turn  Into  the  cowl  It  provide*  the  dual  role  of  an 
external  compreeeton  Inlet  and,  equally  as  important,  the  parabolic  primary  receiver  optics  for  the 
incoming  power  beam.  Roughly  half  of  the  cross  sectional  area  represented  by  the  inlet  centerbody 
Is  covered  with  the  primary  receptor  mirror. 

The  forward  inner  rim  of  the  cowl  contains  twelve  small  circular  secondary  optica  (ref. 
conflg.)  which  art  positioned  dose  to  the  focal  point  of  the  large  parabolic  primary  receiver  mirror. 
These  secondary  mirrors  then  project  the  high  power  beam(s)  parallel  to  the  afterbody  external 
surface  -  at  an  intensity  of  5  x  10*  Wfcm*  (see  Fig.  D-2).  The  cowl  itself  provides  the  function  of 
(Greeting  the  com  pirns  id  inlet  airflow  so  that  it  "blows*  across  tns  afterbody  -  for  the  purpose  of 
"refreshing*  the  spent  air  in  the  detonation  tone.  The  cowl  ia  assumed  able  to  translate  axially  so 
that  optimal  inlet  flow  characteristics  can  be  maintained  throughout  the  flight  envelope;  the  cowl 
may  alao  have  a  deformable  forward  lip  and  trailing  edge.  This  variable-geometry  cowl  would  ideally 
keep  the  conical  bow  shock  always  attached  at  the  cowl  lip  -  above  some  design  mach  number  (e.g., 
mach  3). 
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The  afterbody  fonts  tba  ba—  of  Um  vahide  tad  behaves  u  an  extsmal-axpa— ion  conical 
pln(  nook  of  4#*  half  anyk.  ft  ft  bara  that  tbs  LSC/LSD  waves  dethrar  their  impufts,  tinea  the 
piuy  afto  aenrao  at  tha  ERH  thnatar  aurfaca.  la  tha  rate—  coaflyuratioo,  twelve  propeyatiny 
LSD  wsvae  (a— 1  drcnlar  croaa>  taction,  1  cb  radhs)  product  twelve  Uyk  pwa— ra  (600  atm),  15 
a  bag  cyflndrteai  phana  val —  (haraaftar  ea&ed  *ftnyere*),  which  —beequ— tiy  expand  laietalhr 
ova  tha  pfan  —la  aurfaca  to  pm mala  th«u at.  Fmh  air  ft  fad  axially  to  tha  soak  afterbody,  u 
diractad  by  tha  cowl;  tha  iataral  rafreah  martini— i  afto  acrraa  —  Important  role. 

As  mentioned  above,  a  ba—  aagla  of  45*  waa  aalactad  for  tha  LSD  mods  to  tha  rafaranca 
eonflyuratioa.  A  ba—  aagta  of  0*  waa  pkkod  for  tha  LSC  mods  tinea  amt  of  tha  impufts  may  ba 
transmitted  to  a  Marty  horisootal  surface  ana  rapraaaatad  by  tha  lowar  (bsmi-aphsrieal)  and  of  tha 
tnmeatad  com  afterbody. 

B)  Engine,  Vahiela  aad  Power  System  Archdtaetww 


Aa  maatioaod  abova,  tha  eaatral  throat  of  thft  tovaatiyntion  ft  a  book  raaaarch  afort  to  to 
tha  phytics  of  ERH  thr—lara.  Howvar,  ft  ordar  to  compota  tha  p  wfnn— ari  of  tha  dovtea  one  Bam 
— adfy  iBjftllii  a  number  of  ‘boundary  eoodWo— *  -  which  an  Intimately  rrtatad  to  tha  sup 
daporttioa  moda  (m*i  LSD  va  LSC  wava),  optical  fee— toy,  aad  j—dymmir  a— mptin—  (tochidtof 
*rrtnah*).  to  othar  words  tha  investigation  ft  built  upon,  or  daiaad  around,  a  apedftc  *wh— to 
conflforation*. 

to  tha  ponuft  of  raaaarch  on  CW  biamid  — fy  cafilBB  oaa  m— t  — adfy  tha  baam  powar, 
degree  of  focutiay,  focal  foomstry  (ay.,  point  va  Baa),  absorption  chambar  di  mentions,  notaia  p- 
omatry  and  working  fluid  conditio—  (pm— re,  tamparatun  and  mam  flow  rata).  Ail  thorn  ‘system- 
rakted”  com idaratio—  hava  an  impact  on  baak  raaaarch  tovmtiyatio— .  On  tha  othar  hand,  such 
ftaoaa  aa  tha  aha  of  propellant  tankage  or  dim— Inna  of  receiving  optica  art  of  littla  concarn  to  tha 
propuftfon  maarchar.  Obvious ly,  all  tha  propaUant  ft  carried  onboard  for  tba  rockets,  and  optic 
dim— lo—  an  driv—  by  tba  miming. 

In  contrast,  —  airbnathtoy  anytoa  has  a  somewhat  difarant  aat  of  *drivtay  ft— aa*.  Ptrbapa 
tha  dominant  oaa  ft  tha  naad  for  aa  fflicfta t  iaftt  tinea  tha  anytoa  worfctoy  fluid  ft  eoHscted  from 
tba  atmoophartc  environment,  hrouyb  which  tha  anytoa  propeift  itself.  For  hiyh  math  numbers 
(e.y.,  hiyh  —parsonic  and  hypanonic),  the  inlet  aka  oft—  dwarb  tha  raft  of  tha  anytoa.  Why? 
Baca—  of  tha  naad  to  maximise  total  praaaura  recovery  within  tha  inlet.  In  fact,  for  hypanonic 
anytoaa  proponed  far  the  "Orient  Expram*  (ay.,  acramjate),  tha  tolet/no— k  cans idaratio—  (fietate 
tha  entire  lowar  airframt  yao metry.  In  —  ca,  all  tha  payload  — d  much  of  tha  fusi  will  ba  boused 
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within  the  external  compr— Ion  inlet  volume. 

Thle  Hut  gives  grant  Incentive  for  future  ‘designers"  of  beam-powered  transatmoepheric 
aeroapaeecraft  to  integrate  the  primary  receptive  optica  with  the  engine  inlet,  vehicle  airframe,  or 
exhanet  oosale.  Depending  on  the  choice  of  beam  wavelength,  traaemittiaf  antenna  diameter  and 
range-to-the-powar  eouree  (or  relay),  the  receiving  antenna  rise  can  be  a  significant  leant.  The 
reference  engine/piatform  configuration  adopted  (or  the  preaent  atudy,  aaaumea  that  the  primary 
optic  ie  integrated  with  a  radially-  aymmetric  ieentropk  tpike  inlet;  an  external-expanrion  nocale 
aleo  doubiea  aa  the  vehicle  afterbody.  Furthermore,  it  la  aaenmed  that  the  intake  io  oullciently 
large  to  contain  the  payload;  the  engine  and  vehicle  are  merged  into  a  ebgle  entity.  The  following 
dlanurion  attempta  a  broad  (fiacuerion  of  the  ehokae  in  engine /optica/ airframe  integration  which 
face  future  deaignen  of  beam-booeted  tranaatmoepberk  ahuttlocraft  -  that  utilise  giiH  IhTtil 

Ten  baric  "platform  conrideratioai"  have  been  identified  (there  may  be  men)  which  influ¬ 
ence  the  detailed  derign  aapecta  of  ERH  thruster*.  In  tom,  thana  conridarariona  afact  the  choice 
of  "boundary  coodWons*  that  muet  bo  apodflod  in  order  to  analyse  the  performance  of  any  given 
ERH  thrurier  concept.  They  are: 

1)  Power  eourte  location  (ground-baead  vs.  pm  based) 

2)  Power  beam  wavelengths  (frequency-twitching) 

3)  Ovaml  vshkk/airfram*  conflguntlon 

4)  night  vector  orientation  (relative  to  beam) 

6)  Airbreathing  inlet  geometry 

6)  Shroud/eowl  configuration 

7)  Receptive  optical  train  geometry 

8)  Thnmter  aurfhee/noarie  geometry 

9)  Thnmt  vectoring  control 

10)  Energy  deposition  aodea  for  thnmt  generation 

11)  Air  "refreah*  acbemea 

Each  of  theae  ERH  thrueter- related  coariderations  will  now  be  (flacuaaod  In  more  detail 
1.  Power  Source  Location  (GB  va.  S3) 

The  beam  power  aource  may  be  baaed  either  in  apace  (SB)  or  on  the  ground  (GB),  and 
utifiae  low  altitude  relay  station.  The  preaent  reference  ayatem  aaaumea  6  to  10  GW  aatolilta  aoiar 
poweretatiooa  in  gaoatationary  orbit  and  multiple  low  altitude  relay  eateltttaa  (e.g.,  at  IIS  km).  The 
system  amumm  the  urn  of  "croea-linka*  and  "down-finks",  but  no  "up-finks”  of  beam  power  from 
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the  ground. 


2.  Power  Bonn  Wavelengths  (frog— acy  switching) 

A  sing)*  wavelength  for  the  power  betm  Is  usually  invoked  in  beainsd-energy  propulsion 
studies.  The  capability  for  frequency  twitching  has  been  mentioned  on  occaakm,  usually  in  connec¬ 
tion  with  enhancing  the  eflldency  of  power  beam  tnnsmbtioa,  or  convenion  to  propolehre  throat. 
For  example,  ahorter  (Le.  laser)  wavelaogthe  an  advantageooa for  longer  tranambrioa  Dnka  between 
receiver  and  tranmnitter  antennae  of  a  given  else.  Laser  power  tranambslon  b  usually  limited  to 
fair  weather,  but  a  high  power  beam  could  also  “burn  through”  clouds,  evaporating  all  moisture  in 
the  path.  Ultraviolet  wavelengths  shorter  than  0  J  pm  are  too  badly  dbalpilsd  by  molecular  and 
particulate  scattering  -  to  be  of  much  use  at  the  bottom  of  the  atmosphere.  However,  if  "down¬ 
linked"  hum  space,  ultraviolet  beams  could  still  be  used  to  advantage  (by  a  beam-  powered  vehicle) 
in  the  upper  reaches  of  the  atmoephere  and  of  course,  throughout  the  vacuum  space  environment. 
In  contrast,  microwave  and  millimeter  radiation  exhibit  the  lowest  atmospheric  attenuation  through 
weather  (clouds  and  all  forms  of  precipitation),  but  wavelengths  much  longer  than  10  cm  yield 
transmltter/recehrer  apertures  that  become  oaweidiy. 

The  present  reference  system  aseumm  that  frequency  ewitchiaghavaiiabla in  the  bandwidth 
from  10cm  through  0.3  pm.  Microwave  wavelengths  (lOcm-O.lcm)  are  used  for  near-vertical  “popup” 
maneuvers,  wherein  the  ehuttlecraft  climbs  (at  low  acceleration)  up  through  weather.  A  low  altitude 
(186km)  transmissive  fresnel-type  microwave  relay  satettte  may  be  used  to  refocus  a  large  diameter 
microwave  beam  (*4.,  7  km),  down  to  the  small  dimensions  of  a  ehuttlecraft.  Once  abovethe  clouds, 
the  shuttle  (and  primary  optic)  b  rotated  to  engage  a  laser  link”  -  by  way  of  another  laser  relay 
satellite  at  186km  -  which  it  chases  all  the  way  into  orbit  (the  Anal  transatmoepheric  acceleration 
run). 

The  characteristics  of  microwave  radiation  sharply  restrict  the  maximum  beam  intensity 
level  which  can  be  collected  upon  the  vehkk/engine  primary  optica.  Thb  “plasma  maintenance* 
inteanity  threeheld  scales  directly  with  wavelength  to  the  negative  0.46  power.  As  a  result,  microwave 
frequencies  are  only  good  for  low  power  accelerations;  alternately  UftOC  frequencies  can  easily  enable 
vehicle  accelerations  beyond  human  tolerance. 

It  b  known  that,  radiation-sustained  plasmas  producsd  by  opticsl  wavelengths  (i.e.,  laser) 
regioa  have  two  distinct  structures,  depending  on  the  Incident  radiation  flux  denrity.  These  struc¬ 
tures  resemble  conventional  detonation  and  deflagration  wavee  generated  by  chemical  reactions.  High 
beam  intensities  crests  laser  supported  detonation  (LSD)  waves;  lower  intensities  laser  supported 
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combustion  (LSC)  wivm.  A  critical  question  k  *D©  timilar  absorption  wave  structures  exist  with 
microwave  radiation?* 

The  analytical  expresrions  derived  for  laser  supported  detonation  wave  velocity  and  peak 
pressure  (by  Raiser)  show  no  frequency  dependence.  As  a  result,  one  might  expect  that  “detonation* 
waves  could  also  be  sustained  by  microwave  power.  However,  the  electric  breakdown  threshold  of 
air  appears  to  scale  with  the  inverse  square  of  the  beam  wavelength,  so  the  threshold  is  much  lower 
for  the  microwave  frequencies  than  for  optical.  In  the  microwave  melon,  breakdown  at  atmospheric 
pressure  occurs  at  a  flux  density  of  1*2  MWfait,  as  compared  with  4-5  GW/m i3  for  10.9  pm 
radiation.  Raiser  indicates  out  that  microwave-  triffered  breakdown  normally  occurs  before  a 
detonation  wave  can  be  generated  at  sea  level  pressure.  The  breakdown  threshold  is  simply  lower 
than  that  intensity  which  can  maintain  a  detonation. time  structure. 

However,  as  will  be  demonstrated  in  Chapt.  VII,  laser-  sustained  “deflagration- type*  waves 
can  indeed  produce  blast  waves,  the  eeeential  ingredient  for  generating  propulsive  force*  with  ERH 
thrusters.  So  there  is  good  reason  to  (aspect  that  microwave  eustiined  ‘deflagration*  waves  can 
hi  i  — flilljr  be  used  ae  well  Mkro wave  radiation  can  sustain  plasma  temperatures  which  are  much 
lower  than  thoee  generated  by  optical  frequandae.  For  microwaves,  reeultant  air  plasma  temperatures 
can  be  on  the  order  of 2000  to  4000K,  whwees  can  it  exceed  10000  to  20000  K  for  optical  frequencies. 
Therefore,  radiation  hast  transfer  to  ERH  thruster  surfaces  might  be  reduced  by  nting  microwave 
power  during  low  acceleration  maneuver*  (Le^  Bft-of  and  landing)  and  hover. 

One  muat  be  careful,  however,  not  to  compare  apples  and  oranges.  As  noted  above,  the 
ERH  thruster  ie  a  renedtivalY  nuhad  engine  which  generate*  thrust  by  expanding  blast  waves  of  the 
impulse-coupling  surface.  The  pulse  repetition  frequency  and  duty  cycle  obviously  must  enter  into 
this  heat  transfer  problem.  Hence  one  cannot  just  examine  peak  microwave- and  laser-heated  plasma 
temperatures,  but  should  carefully  asses  the  convective,  conductive  and  radiative  beat  transfer  to 
the  ERH  thruster  surface. 

9.  Overall  Vehicle/ Airframe  Configuration 

Clearly  there  ahould  exist  an  infinite  variety  of  airframe  choices  for  transatmospheric  aero- 
spacecraft  Historically,  however,  there  have  been  only  two  broad  categories:  i)  delta-  wing,  and  U) 
exi-iymmetric.  The  first  category  is  exemplified  by  the  space  shuttle  orbiter,  and  perhaps  by  several 
earlier  lifting-body  reeearch  vehicles  (glide  re-entry  mistion).  Of  more  recent  vintage  are  numerous 
delta-winged  concepts  for  the  "Orient  Express*,  propelled  by  hypertonic  acramjet  engines. 

The  other  major  clam  ie  exemplified  by  the  radially  symmetric  Apollo,  Gemini  end  Mercury 
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■pace  capsules,  and  certain  axi-symmetric  experimental  engines  like  the  NASA  Hypertonic  Research 
Engine  (HRE).  The  present  reference  configuration  la  ammtd  to  be  asd-qrnunetric,  and  has  a 
paraboBc  primary  receptive  optica  intimately  integrated  with  the  large  external-compression  inlet 
centerbody. 


4.  night  Vector  Orientation  (relathre  to  beam) 

In  theory,  a  beam-powered  shuttlecraft  could  be  propelled  in  any  fight  direction  relative 
to  the  projected  energy  beam.  For  transatmoepheric  mbtions,  some  fixed-geometry  laser-heated 
rocket  schemes  have  been  propoeed  to  utilise  a  large  range  of  beam  intercept  angles  (e.g.,  up  to 
CO*);  these  concepts  do  not  require  focusing  optics  in  the  propulsion  convertor,  however.  For  those 
that  &  require  precWon  receptive  optics,  two  alternatives  exist:  i)  fixed  geometry,  and  ii)  variable 
geometry. 

A  good  example  of  the  second  class  la  the  articulated  “telescope"  receptive  optical  train 
propoeed  by  TRW  for  their  laser-heated  Ht  rocket-powered  orbit  transfer  vehicle  (OTV)  concept. 
Clearly,  in  the  vacuum  of  space,  this  telescope  arrangement  gives  complete  freedom  in  the  fight 
vector  orientation  (+/  - 180*).  With  no  atmospheric  drag  penalty  in  a  too  “G*  environment,  these 
structures  could  be  both  large  and  fight-weight. 

Worlds  apart  from  the  OTV  minion  is  that  of  an  airbreathing  TAV  powered  by  beamed 
energy.  In  the  TAV  misnon,  articulated  primary  optic  geometries  for  beam  reception  must  be 
abandoned,  unless  small  enough  to  be  placed  beneath  a  transparent  high-power  dome/window  - 
that  also  has  low  drag.  The  remaining  fixed-geometry  alternative  leaves  two  options  for  filght/beam 
vector  orientation:  Ijaxial,  or  lijlateral  In  the  former,  the  vehicle  fight  rods  is  aligned  with  the 
beam,  and  the  engine  “swallows*  the  beam  as  it  climbs  to  orbit.  In  the  latter  option,  the  vehicle 
accelerates  at  right  angles  to  the  beam  (or  at  tome  fixed  angle). 

Both  option*  might  need  to  invoke  a  low  altitude  relay  satellite  in  order  to  negotiate  the  pitch 
change  schedule  required  along  the  boost  trajectory.  A  viable  alternative  is  to  switch  through  several 
relay  stations  in  strategically  placed  orbits,  to  accomplish  a  series  of  greater  vehicular  pitch  changes. 
It  should  be  noted  that  a  large  primary  optic  which  is  permanently  affixed  to  a  TAV  shuttlecraft  can 
probably  not  sustain  a  beam  misalignment  (pointing  accuracy)  of  more  than  1-3  degree*.  However, 
moat  studies  of  rocket  engine  thrust  vector  data  gathered  on  current  launch  vehicle*  (including  the 
•pace  shuttle  orbiter)  indicate  that  2-3  degrees  is  typically  what  is  demonstrated. 

The  present  study  ammo*  that  the  vehicle  axis  of  symmetry  is  aligned  with  tbs  fight 
direction  to  an  "accuracy  of  +  /-  2  or  3  degrees*  -  to  permit  thrust  vector  control  of  the  fight  path 
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angle. 


5.  AirfarMtUag  Inlet  Geometry 

There  to  no  guarantee  that  a  revolutionary  abbreathing  engine  concept  such  a*  the  ERH 
thrueter  would  necessarily  derive  benefit  from  the  addition  of  a  eupenonk  inlet  of  conventional 
derign  (or  geometry).  This  toeue  to  especially  worth  contidering  becauoe  it  to  the  iow  euperaonic  and 
eubeonk  mach  number*  where  the  ERH  thrater  to  expected  to  exeell  (e.gM  Mach  0  to  S  for  the 
reference  configuration).  However,  other  factors  related  to  multi-mode  or  combined-cycle  operation 
can  be  expected  to  enter  into  the  requirement*  for  defining  the  ultimate  airbreathing  inlet  geometry. 
For  example,  in  the  reference  configuration  the  primary  propulsive  role  may  transition  from  ERH 
thruster  to  ramjet  at  Mach  3,  then  to  scramjet  at  Mach  ( (or  quite  possibly  into  scramjet  at  Mach  3, 
skipping  the  ramjet  mode  altogether).  These  latter  cycles  may  require  a  cowl  for  best  performance, 
and  the  ERH  thrater  should  benefit  from  the  additional  air  supply  for  axial  “refresh*. 

Early  studies  by  BUlig  (3],  and  others  have  proven  that  seramjets  with  conventional  inlets 
and  cowls  to  contain  ths  “combustor”  section,  yielded  superior  propulsive  efficiencies  to  un-cowled 
or  “external- burning  ramjet*  concepts.  However  these  studies  examined  the  use  of  conventional 
chemical  fuels  in  hypersonic  cruise  modes.  Hence,  the  conclusions  may  require  some  interpretation 
when  applied  to  radically  different  advanced  engine  concepts  designed  for  the  acceleration  role  - 
utilising  beamed  energy.  Furthermore,  one  must  consider  that  a  beam-powered  TAV  engine  will  be 
in  this  propulsive  mode  for  a  very  short  time  (e.g.,  less  than  60  seconds). 

Similar  statements  could  be  made  for  the  ERH  thruster  since  it  to  an  unusual  repetitively- 
pulsed  engine,  and  fails  in  a  completely  different  darn  than  the  more  familiar  ones  using  Brayton 
“constant  pressure”  cycles.  One  thould  note  that  the  compression  process  (in  the  ERH  thruster)  to 
provided  by  a  propagating  LSD  or  LSC  wave,  and  not  a  mechanical  compressor.  Also,  the  spent 
air  must  be  vented  and  replaced  with  new  air,  and  two  mechanisms  have  been  identified  that  could 
provide  this  function:  i)  axial  refresh,  and  ii)  lateral  refresh  (discussed  below,  in  another  section). 
Hence,  the  unusual  engine  cycle  exhibited  by  the  ERH  thruster  concept  demands  that  one  take  a 
completely  new  look  at  the  available  options  for  air  inlet,  power-injection,  and  exhaust  mechanisms. 

As  to  commonly  known,  all  fixed-geometry  supersonic  inlets  have  a  design  Mach  number 
(Mo)  whereupon  the  bow  shock  to  exactly  attached  at  the  inlet  cowl  lip.  Off-design  operation  forces 
the  inlet  into  a  sub-critical  or  super-critical  condition  wherein  the  total  pressure  recovery  can  be 
expected  to  decline. 

Variable-geometry  inlets  attempt  to  extend  tbs  range  of  optimum  performance  beyond  just 
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Mo.  Clearly,  the  approaches  to  varying  the  inlet  geometry  are  different,  depending  upon  whether 
the  configuration  baa  2-dlm enaional  (2*D)  or  radial  symmetry  (iaentropic  spike,  bi-conic,  tri-conic, 
etc.).  The  present  study  will  focus  will  be  upon  the  latter,  since  it  is  used  on  the  reference  platform. 

There  appear  to  be  at  least  six  methods  (see  Fig.  II-3)  which  (singly  or  jointly)  have  the 
potential  for  preventing  a  conical  bow  shock  from  becoming  lodged  (or  swallowed)  into  the  inlet 
cowh  i)  translating  the  shrould  fore  and  aft  (e.g.,  like  NASA's  HRE  scramjet),  ii)  deforming  the 
cowl  leading  edge  (Le.,  like  retractable  "high  lift  slats”  on  wings),  ill)  replacing  the  conical  tip  with 
an  ogive  shape,  hr)  extending  a  pointed  spike  forward  out  the  cone  tip,  scheduled  with  increasing 
flight  mach  number,  v)  "ticker  platen”,  and  vi)  blunted  hemispherical  note  caps. 

The  first  option  has  a  potential  problem  in  that  the  secondary  optic  may  be  forced  to 
translate  fore  and  aft  with  the  shroud,  thereby  removing  them  from  the  primary  optica  focal  region. 
(Alternatively,  shroud  translation  can  be  used  to  engage  a  sequential  series  of  secondary  optics.)  The 
second  option  may  increase  the  mechanical  complexity  of  the  cowl,  but  should  be  feasible  for  cowls 
having  a  large  diameter  and  relatively  short  length.  The  third  option  may  spill  a  lot  of  air,  causing 
increased  "spillage  drag”  for  the  inlet.  It  ia  not  clear  whether  or  not  the  fourth  option,  extending  a 
spike,  will  deflect  enough  air  at  supersonic  speeds  to  do  any  good.  The  fifth  option,  ticker  plates, 
has  been  luccemfuliy  employed  with  2-D  inleta  in  almost  all  current  Mach  2+  jets.  The  sixth  option 
of  blunting  the  noae  on  the  conical  tip  (e.g.,  with  a  hemisphere),  would  result  in  a  detached  shock 
forming  just  off  the  noae.  Although  the  flow  could  be  expected  to  re-attach  later  on  the  ogive  or 
conical  Inlet  surface,  some  lose  in  total  pressure  recov«y  would  be  incurred. 

0.  Shroud /Cowl  configuration 

At  least  three  mechanisms  (see  Fig.  11-4)  can  be  invoked  for  turning  the  engine  air  working 
fluid  around  a  large  external-compreaaion  inlet  centerbody  -  so  that  the  ERH  thruster  surface  maybe 
“refreshed”  axially:  1)  a  physical  shroud  with  a  variable  exit  area  noule,  ii)  a  shroudlese  configuration 
that  leavee  Prandtl-Meyer  expansion  fans  to  turn  the  flow,  and  iii)  a  thin  mini-shroud  “ring”  which 
solely  consists  of  a  small  annular  secondary  (reflective)  optica. 

The  flrat  option  is  the  traditional  approach  (moet  conservative),  and  could  be  counted  upon 
for  enabling  good  total  pressure  recovery  throughout  the  entire  Mach  0-3  envelope  (required  of  the 
ERH  thruster  in  the  reference  platform).  Clearly,  this  jet  of  pressurised  air  leaving  the  shroud  could 
supply  moet  of  the  ERH  thruster  surface  with  “refreshed”  air  at  densities  much  above  ambient.  It 
would  also  permit  the  engine  to  develop  thrust  levels  beyond  that  which  could  be  demonstrated  in 
the  static  (l.e.,  motionless)  condition  at  any  given  altitude.  However,  this  favorable  situation  is  not 
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Figure  U-3  Methods  for  incorporating  variable  geometry  in  radlally-eymmetric  lupereonlc  Inlets 
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secured  without  suffering  the  penalties  of  in  crewed  dreg  produced  by  the  annular  shroud  end  redial 
support  strata,  and  spilled  ram  air. 

The  second  option  avoids  shroud  drag  altogether,  by  eliminating  the  shrould  itself.  However, 
the  penalties  for  this  approach  will  appear  in  several  other  areas.  During  high  supersonic  speeds,  the 
ERH  thruster  surface  will  be  immersed  in  "refresh  air*  at  a  lower  density  than  with  the  shrouded 
option.  Hus  condition  translates  into  a  reduction  in  maximum  attainable  thrust  Abo,  with  small 
"a"  plug  noules  (Le.,  nearly  flat-  bottomed  configurations),  flow  separation  could  occur  immediately 
aft  of  the  maximum  engine/vehicle  cross-section.  Hence,  the  resultant  low  densities  behind  the 
vehicle  would  not  only  produce  greatly  increased  drag,  but  adversely  affect  thrust  generation  as 
well.  Another  potential  problem  b  the  complex  flow  field  that  would  result  from  separation  and 
vortex  shedding,  off  the  aft  ERH  thruster  surface.  This  may  have  a  deleterious  effect  on  the  LSD  wave 
ignition  and  propagation  processes.  Finally,  there  may  be  other  significant  losses  due  to  vorticity  in 
the  exiting  flow  as  it  expands  away  from  the  vehicle. 

In  the  third  option,  the  mini-shroud  b  simply  a  small  annular  secondary  mirror,  having  an 
aspect  ratio  not  much  different  from  a  "hulla-hoop*  -  only  with  a  bi-convex  cross-  section.  Although 
thb  thin  miniature  airfoil  could  induce  some  "aerodynamic  lift”  flow-turning,  the  bulk  of  the  work 
would  be  accomplished  by  high  velocity  jet  flow  fields  which  trail  behind  propagating  LSD  waves  (2] 
-  as  shown  in  Fig.  H-5.  Pro-compressed  air  that  enters  the  "cowl*  will  have  been  turned  outward 
(e.g.,  38*  to  40*  from  axial)  by  the  external  inlet  forebody.  Next,  hypersonic  LSD  waves  would  chop 
repeatedly  through  thb  annular  "duct”  flow.  After  much  turbulent  mixing,  the  annular  cowl  air 
flow  would  be  vectored  back  past  the  axial  Erection  (perhaps  aided  by  an  "ejector”  process),  and 
expanded  against  the  engine/vehicle  afterbody  plug  noule. 

The  obvious  penalty  incurred  here,  will  be  due  to  the  additional  beam-power  invested 
in  longer  laser  pubes  (le.,  to  turn  the  flow  with  LSD  waves).  Normally,  the  laser  pulse  would 
be  terminated  just  as  the  LSD  waves  reach  the  uppermost  perimeter  of  the  ERH  thruster  "plug”, 
right  before  they  enter  the  annular  shroud  air  gap.  Clearly,  thb  type  of  flow-turning  mechanism 
would  be  accessible  only  to  receptive  optical  trains  which  incorporate  "ring-type*  secondary  optics. 
In  contrast,  a  slightly-scalloped  primary  optics  which  links  to  a  set  of  small  circular-type  secondary 
optics,  could  obviously  aai  invoke  the  LSD  wave  flow-turning  mechanbm. 

7.  Receptive  Optical  Train  Geometry 

Taken  as  a  whole,  the  receptive  optical  train  comprises  the  receiving  "antenna”  which  collects 
the  remote-beamed  power  and  delivers  it  to  the  ERH  thruster  surface.  In  general,  the  train  requires 
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only  Primacy  Optica  (PO)  and  Secondary  Optica  (SO),  but  in  aomo  eaaaa  may  include  Tertiary 
Optica  (TO).  Por  brevity,  tin  following  dictation  b  reetricted  to  engbe/veMci*  configurations  with 
radiai  aymmetry. 

Aa  mentioned  in  the  previous  taction,  radially -symmetric  engines  cam  i)  be  shrouded,  H) 
be  shroud-baa,  or  ifi)  uae  a  mlal-shroud.  Figure  II-4  displays  repraaentative  receptive  optical  tralna 
which  are  compatible  with  theae  three  categories;  it  b  not  meant  to  be  aa  exhauerive  treatment 
of  the  aubject.  For  the  purpoee  of  comparison,  the  reference  platform  b  ahown  in  F1|.  D-fla.  The 
large  receptive  primary  optic  b  a  paraboloid  of  revolution.  If,  for  example,  aa  mulit  high-intensity 
beam  muat  be  delivered  acroe*  the  ERH  thruater  surface,  then  the  secondary  optic  muat  alao  be 
a  rmitwiuaiM  fine  mounted  to  the  inner  surface  of  the  shroud,  just  short  of  the  annular  line  focua 
(i.e.,  of  the  parabolic  PO).  If  a  large  number  of  cylindrical  detonation  wave  sources  (Le.,  Ungers*) 
b  desired  screw  the  ERH  thruater,  the  primary  optic  surface  muat  be  slightly  scalloped.  Thb  will 
concentrate  and  spfit  the  power  beam  onto  a  aerie*  of  email  circular-shaped  secondary  optics  - 
which,  in  tun,  will  reflect  the  beam  in  a  multi- fingered  fashion.  Thb  option  may  be  referred  to  w 
a  multi-point  focua. 

Finally,  it  should  be  mentioned  that  If  the  ehroud  b  allowed  to  translate  fore  and  aft, 
—ratal —is  of  the  email,  lightweight  SO  mirrors  could  be  engaged  -  either  at  dberat*  umgyiiK 
Sight  Mach  no.’i  (perhaps  during  transition  to  Afferent  engine  modea),  or  poaaibly  aa  a  result  of 
power  beam  frequency-shifting.  In  the  latter  case,  larger  SO  apertures  may  be  required  with  longer 
beam  wavelength*  in  order  to  counteract  diffraction  effects.  In  general,  the  SO’e  may  be  designed 
to  project  colQmnated  beams  with  roughly  constant  intensity  along  the  ERH  thruater  plug. 

The  unshrouded  engine/vehide  configuration  shown  in  Fig.  D-db  takes  a  slightly  Afferent 
approach  to  the  optical  train  design.  Here,  a  high  power  window  b  placed  over  the  parabolic 
reflective  primary  optic.  Ciearty,  there  will  be  losses  suffered  by  the  beam  as  it  passes  through  thb 
transparent  aeroehelL  Upon  reflecting  off  the  PO,  the  beam  b  focused  towards  the  SO,  then  traveb 
onward  to  the  ERH  thruster  surface.  Either  annular-line  or  multi-point  foeii  can  be  invoked  for  use 
with  truncated  plug  nozzles.  Alternately,  one  might  use  primary  lenaea  instead  of  the 

reflective  variety,  and  eliminate  one  set  of  optical  surfaces  (i.e.,  the  high  power  'planar”  window). 

Figures  H-6c  and  6d  portray  the  mini-strut  geometry,  which  is  a  variation  of  the  mini-ahroud 
theme  shown  earlier  in  Tigs.  ll-4c  and  II-S.  If  the  reflective  PO  ia  slightly  scalloped,  then  the  annular 
SO  may  be  replaced  by  a  set  of  circular  SO’s  affixed  to  the  end  of  radial  mini-struta.  ineidently,  theae 
external  struts  could  be  incorporated  into  other  engine  modea  and  thereby  serve  multiple  functions. 

Fig.  U-6e  illustrates  one  final  option  for  the  receptive  optical  train  geometry,  one  which 
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Note 

PO  can  be  parabolic 
■urface  to  produce 

annular  Dm  focus, 
or  Scalloped”  to 
produce  multi-point 
focus  -  at  the  rite 
of  the  secondary 
optk(s). 
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d)  urfat-otrat  (top  rWw) 


Figure  Q-4J  Optloaa  for  Racaptiva  Optical  Traiaa,  coat’d. 


dsm oneleaass  tktwtf  lrt»]  optic*,  b  this  example,  the  pmktlc  PO  Ini  Uww  the  power 
bwa  out  towirdiM  ihmIk  Bne  (beta;  too  tgijlggfldJIBnllESO  intercepts  to  benm  tad  rtoete 
U  towards  m>  irf  small  graatog-inddanc*  TO’a.  As  Indkstod  in  Fig.  B-4e,  to  flan  top  id  ton 
between  nek  SO  and  TO  (see  Fig.  0-W)  would  b«  concentrated  in  only  sm  dktetba.  Kadi  TO 
would  be  of  approxlma My  cyfiadrkal  contour  and  rtAtcl  a  circular  beam  across  to  ERE  thruster 
ph * 

Tbia  geometry  would  enable  to  um  of  multiple  detonation-  Bne  aourcea  (Le.,  “plasma 
flngert*),  yet  not  demand  to  uae  of  a  acalloped  PO  auiface.  Why  uoe  thk  three  element  optical 
train  when  it  appeaia  more  eompfieated  than  to  two-element  system?  Because  other  combined- 
cycle  engine  modes  may  require  a  diferent  combinatioo  of  leaaaa  -  aecemitathu  an  optical  train  with 
variabte-feometry.  For  example,  the  TO'i  in  Fig.  II-4b  are  aufldentiy  email  to  bo  sadly  retracted. 

Finally,  Fig.  Q-7  portrmyi  a  top  view  of  the  PO  aurface  for  to  reference  point  (iedgi  Figure 
D-7a  ahowa  to  iaaer  illumination  pattern  uaed  to  generate  to  nmultaneoaa  12  Bne  ( Unger')  source, 
and  Fig.  II- 7b  for  to  rotating  angle  Bne  so  urea 

S.  Thruster  SwrCnco/Nanab  Qwomatey 

b  general,  to  vehicle  afterbody  on  future  beam-powered  reusable  SSTO  shuttlecraft  will 
need  to  serve  multiple  functions,  b  to  reference  radlalty-eymmetrlc  engbe/ vehicle,  to  afterbody 
ia  used  aa  a  nhif  flftltk  type  of  external-exp  and  on  aurface  for  to  entire  combined-cycle  engine.  This 
multi-mode  propuldon  system  ia  composed  of  the  ERB  thruster  from  Mach  0  to  3,  several  other 
airbreathing  modes  from  Mach  3  to  2S,  and  finally  rockets  for  orbit  Insertion.  The  afterbody  alao 
acts  aa  a  re-entry  beat  shield.  Hence  it  is  Ukeiy  to  be  an  uncooied  surface,  and  probably  fabricated 
from  advanced  refractory  materials  like  to  carbon/carbon  composites  used  on  the  Space  Shuttle 
Orb  iter  note  and  leading  edges. 

a)  Plug  nottle  varieties 

It  ia  logical  to  aaaume  that  the  plug  nouie  should  be  combined  with  the  ERH  thruster 
surface  itself  because  of  the  need  for  “axial  refresh*  -  a  replenishment  of  the  'spent'  air  heated  by 
the  previous  LSD  wave.  The  specific  variety  of  plug  noaile  used  in  the  reference  design  is  shown 
in  Fig.  O-fc  Le.,  a  parti  el -in  tern  il/ partial -external  expansion  type.  For  tide  geometry,  a  portion 
of  the  flow  expand  on  will  oecur  within  the  cowl;  the  minimum  cross  sectional  duct  area  would  bo 
upstream  of  to  eowi  exit  station. 

Figure  11-9  graphically  demonstrates  the  efect  on  nouie  (and  cowl)  contours  aa  to  per 
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Figure  n-7  Top  view  of  PO  «urf»c#  for  ref.  pots*  dotigB 
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eentage  of  internl  exposition  is  increased.  Not*  that  the  initial  angle  which  the  flow  makes  with 
the  central  axis  of  eymmetry  (l.eM  a)  gets  larger  with  increased  internal  expansion.  Eventually, 
with  100  %  internal  expiation,  the  noule  becomes  inverted’,  end  ill  flow  expansion  takes  place 
against  the  inner  surface  (see  Fig.  0-10).  This  variety  of  noule  is  sometimes  referred  to  as  an 
Expansion-Deflection  (ED)  noule. 

Since  laser  power  can  only  be  projected  along  straight  lines,  the  external  plug  noule  surface 
must  be  approximated  by  a  simple  cone.  For  example,  Fig.  0-11  shows  the  outline  of  a  50%  internal 
expansion  plug  noule  and  demonstrates  how  the  complex  isentropic  spike  contour  can  be  replaced 
with  a  variety  of  conical  surfaces.  Upon  further  examination,  it  becomes  obvious  that  one  must  select 
the  larweat  indicated  cone  apex  angle,  in  order  to  prevent  high  intensity  laser  beams  (projected  by 
the  SO  mirrors)  from  impinging  upon  the  plug  noule  surface.  For  the  case  given  in  Fig.  11-11,  "a* 
would  roughly  equal  40*,  which  is  also  the  angle  that  cowl-ducted  gas  would  emerge  upon  the  plug 
noule.  (Inddently,  since  the  reference  engine/vehicle  has  an  a  of  45°,  one  might  expect  that  25% 
to  50%  of  the  flow  expansion  would  take  place  within  the  shroud.) 

b)  Importance  of  base  angle 

One  important  parameter  to  consider  in  the  design  of  such  conical  plug  nouies  (Le.,  for  use 
as  ERH  thruster  surfaces),  is  the  base  angle  a  indicated  in  Fig.  D-12.  As  the  base  a  is  increased 
so  does  the  length  *L*  of  the  conical  ERH  thruster  surface  (i.e.,1  =  xj  cos  or  ).  This  creates  a 
longer  “blast  tone"  for  the  input  of  energy,  as  well  as  additional  surface  area  over  which  the  blast 
wave  can  expand.  However,  the  additional  area  does  not  always  produce  more  propulsive  force, 
because  paly  the  vector  component  aligned  with  the  vehicle  flight  direction  results  in  “thrust”  (i.e., 
usable  thrust  =  T  cos  a).  By  increuing  the  base  angle  a  the  efficiency  of  the  ERH  thruster  will 
actually  decrease,  since  more  laser-generated  impulse  is  required  for  an  equivalent  propulsive  thrust. 

Neverthelem,  when  designing  the  overall  engine/ vehleie  configuration,  many  compromises 
come  into  play,  and  the  above  discussion  is  just  one  example.  Ooe  should  note  that  a  decrease  in 
thruster  efficiency  (due  to  large  a)  may  to  some  extent  be  offset  by  a  decrease  in  aerodynamic  drag, 
and  an  improvement  in  “refresh”  function. 

In  the  search  for  an  ideal  thruster  base  angle,  a,  one  should  also  consider  plug  noule 
performance  in  the  rocket  mode  -  which  is  an  important  isaue  for  combined-cycle  engines.  Note  that 
the  extreme  am  ratio  plug  noule  shown  in  Fig.  11-13  would  give  excellent  expansion  characteristics 
for  rocket  engines  initiated  above  200,000  ft  altitude.  Shallow  thruster  angles  (e.g.,  or  *  30*) 
would  give  exceptional  performance  in  the  ERH  thruster  mode  as  well.  Throughout  the  intervening 
airbreathing  modes  from  Mach  3  to  25,  ducted  air  flow  must  be  prevented  from  separating  at  the 
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Figure  n-9  Noiile  contour*  r*.  percent  Internal  expansion 
(From  Graham,  Ret  #6) 
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Figure  H-10  3-D  Prandtl-Meyer  expuutou  (constant  area  ratio,  con* taut  thrust) 
|  (From  Graham,  RsL#6) 
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Figure  H-IJ  Effect  of  baee  angle 


extreme  cowl  «odt  aagi**,  of  court*. 

c)  Blast  wave  overlap  problem 

In  the  aerothermodynamic  design  of  ERH  thruster  surfaces,  it  is  important  to  consider  the 
degree  of  overlap  by  adjacent  cylindrical  blast  waves  as  they  expand  down  to  local  ambient  pressure, 
against  the  conical  plug.  Expressed  another  way,  one  must  consider  the  ratio  of  “active*  to  “passive* 
thruster  am  available  over  the  entire  thruster  plug.  Here,  “active”  am  is  defined  as  that  plug  am 
which  lies  immediately  above  an  LSD  wave  as  it  races  up  the  entire  length  of  the  nozzle.  “Passive" 
area  is  that  lateral  am  over  which  the  resultant  cylindrical  blast  wave  expands,  in  the  process  of 
delivering  impulse  to  the  plug. 

The  maximum  lateral  radius  to  which  a  cylindrical  blast  wave  expands  (when  the  pressure 
within  it  reaches  ambient)  is  a  critical  design  parameter  for  efficient  “lateral  refresh*.  If  significant 
overlap  of  adjacent  blast  waves  occurs,  then  the  lateral  refresh  process  will  become  impaired. 

As  the  engine/vehicie  climbs  in  altitude,  this  lateral  “refresh  impairment*  will  first  occur  at 
the  lower  point  of  a  complete  conical  plug  nozzio.  Hence,  ERH  thruster  efficiency  can  be  enhanced 
by  truncating  the  cone.  Inddently,  if  the  ERH  thruster  surface  must  also  serve  as  a  re-entry  heat 
shield,  the  cone  tip  be  rounded  to  a  large  spherical  radius  anyway  (see  Fig.  H-12). 

One  logical  way  to  prevent  the  lateral  refresh  impairment  problem,  (even  for  truncated  cone 
plugs)  would  be  to  utilize  a  conically  focused  “multi-finger”  beam  geometry  across  the  ERH  thruster 
face  -  instead  of  the  cylindrical  “multi-finger*  beam  geometry  used  in  the  reference  design.  This 
means  that  the  plasma  fingers  will  actually  be  tapered;  hence,  the  plasma  will  relax  to  ambient 
pressure  faster  at  the  lower  end  of  the  plug  (i.e.,  at  the  finger  tips),  than  at  the  top.  Since  more 
room  exists  for  lateral  expansion  at  the  top  of  the  truncated  plug  than  at  the  bottom,  this  would 
seem  to  be  an  optimum  solution. 

Then  in  order  to  maintain  a  constant  beam  intensity  upon  propagating  LSD  waves  (required 
for  a  homogeneous  detonation  pressure),  beam  power  must  be  increased  (.i.e.,  ramped  upward) 
during  the  laser  pulse.  This  will  result  in  a  “sawtooth-shaped”  pulse  energy  distribution  vs.  time. 
Since  the  primary  receptive  optics  has  a  fixed  area  over  which  to  receive  the  power  beam,  a  “saw¬ 
toothed”  laser  pulse  would  more  severely  load  the  PO  at  the  sod  of  the  laser  pulse,  than  at  the 
beginning.  This  condition  might  necessitate  a  reduction  in  received  time-averaged  beam  power, 
below  that  which  is  assumed  for  the  reference  design.  (Incidently,  a  homogeneous  detonation  pressure 
field  might  be  advantageous  from  an  engineering  design  perspective,  because  the  entire  plug  surface 
would  have  the  same  load  per  unit  area.) 


d)  MhcWlaa* >■  «w 

The  available  receptive  PO  mirror  area  and  peak  beam  power  can  be  easily  computed  for  a 
given  E&H  thruster  plug  -  once  the  entire  enfine/vehide  geometry  has  been  defined.  (This  subject 
is  treated  thoroughly  in  Section  C  below,  and  will  not  be  repeated  here.)  It  bears  mentioning, 
though,  that  the  available  PO  area  can  be  directly  related  to  the  choice  of  forebody  Inlet  geometry 
(he.,  as  dominated  by  inlet  the  cone  apex  angle),  inlet  design  Mach  number  (Mb),  and  angle  a  of 
the  afterbody  plug  noale. 

In  designing  an  ERH  thruster,  one  would  be  Ul-advised  to  ignore  the  relationship  of  the 
shroud  to  the  expanding  blast  waves  shed  by  the  ERH  thruster  plug.  With  some  thought,  one  could 
imagine  a  shroud/phig  geometry  wherein  the  shroud  helps  to  direct  the  rapidly  expanding  plasma 
aft,  yielding  increased  thruster  efficiency. 

As  a  final  note,  ERH  thruster  surfaces  will  be  exposed  to  radiation-sustained  plasmas  at 
temperatures  of  roughly  4000  K  with  LSD  waves,  and  greater  than  10,000  K  with  LSC  waves.  Such 
elevated  temperatures  may  cause  the  thruster  surface  to  ablate  (a  potential  materials  problem). 
Furthermore,  LSD  waves  will  produce  very  high  peak  pressures  (e.g.,  600  atm  in  the  reference  design) 
upon  very  localised  areas,  resulting  in  extremely  high  stresses.  Moot  conventional  high  temperature 
materials  (e.g.,  ceramics)  are  unable  to  withstand  such  high  localised  stresses.  In  contrast,  ductile 
materials  (refractory  metals)  which  can  survive  large  non-uniform  stress  distributions,  lose  strength 
when  heated  to  extreme  temperatures. 

Actual  material  selection  for  the  reference  ERH  thruster  configuration  was  not  attempted. 
However,  such  advanced  materials  are  likely  to  become  available  In  the  foreseeable  future,  if  they 
are  not  already  present:  e.g.,  carbon/carbon  re-  entry  heat  shields,  or  whisker- reinforced  ceramic 
composites. 


9.  Thrust  Vectoring  Control 

The  final  topic  of  this  section  involves  the  need  for  “thrust  vectoring”,  a  necessary  element  in 
the  control  of  vehicle  flight  path  angle.  The  ERH  thruster  utilises  a  “segmented*  type  of  operation, 
which  automatically  offers  a  unique  mechanism  for  thrust  vector  control  -  without  the  need  of  a 
gimbal  system.  By  increasing  and  decreasing  the  time-averaged  pressure  across  various  segments, 
the  necessary  forces  for  controlling  the  thrust  vector  may  be  generated.  The  time-averaged  pressure 
over  any  segment  angle  of  the  plug  may  be  varied  simply  by  changing  the  laser  intensity  delivered 
to  the  LSD  waves  propagating  across  that  segment.  This,  in  turn,  varies  the  detonation  pressure 
[PD w),  and  therefore  impulse  delivery. 
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Tha  thruat-vactoring  control  procM  can  bo  accomplish  ad  in  thm  additional  ways.  First, 
the  radial  symmetry  of  the  blast  waves  may  be  disrupted  to  asymmetrically  push  one  side  of  the  pluf 
more  often  than  the  other.  Second,  at  indicated  in  Fig.  11-14,  the  LSD  wave  ignition  point  may  be 
vailed  up  or  down  the  plug  to  decrease  or  Increase  the  cylindrical  length  of  the  blast  waves.  Third, 
the  laser  pulse  duration  could  be  shortened  in  order  to  restrict  impulse  delivery  (preferentially)  to 
the  lower  end  of  the  thruster  plug. 

It  may  be  necessary  to  configure  the  plug  geometry  such  that  the  impulse  generated  by  a 
single  full-length  “line  source*  blast  wave  is  directed  along  a  vector  that  pasees  exactly  through  the 
vehicle  center-of-gravity.  This  condition  would  enable  the  vehicle  to  propell  itself  «i  right  angles  to 
the  laser  power  beam,  perhaps  a  necessary  control-related  function  for  a  beam-riding  engine/vehicle. 

In  Fig.  11-15,  it  is  seen  that  the  time-averaged  pressure  across  the  high  pressure  side  of  the 
plug  is  defined  as  p«*  and  that  on  the  low  pressure  side  as  pw,  which  acts  over  the  segment  angle  V* 
of  the  plug.  In  general,  when  p<*  is  increased,  Pei  is  decreased  in  such  a  way  as  to  maintain  constant 
total  thrust.  In  normal  Sight,  the  entire  thruster  plug  would  operate  with  the  same  time-averaged 
pressure,  defined  as  pM. 

When  thrust  vector  control  is  obtained  by  this  technique,  the  vector  can  be  made  to  tilt 
over  through  the  angle  *6"  and  also  shift  laterally  through  a  distance  “d".  Normally,  the  direction 
of  vector  angle  and  vector  shift  would  be  additive,  in  so  far  as  generating  a  turning  moment  about 
the  vehicle  center  of  gravity.  Finally,  the  vector  angle  can  be  expected  to  depend  upon  at  least  three 
parameters:  i)the  percent  overpressure,  (p<*  -  pM)fpm',  ii)operating  pressure  ratio,  pM/p  (where 
Pt  a  ambient  atmospheric  pressure),  and  iii)piug  noule  geometry  itself. 

10.  Energy  Deposition  Modes  for  Thrust  Generation 

This  section  offers  a  broad  perspective  of  the  various  energy  deposition  modes  which  can  be 
used  by  ERH  thrusters.  In  describing  these  modes,  the  following  physical  principles  come  into  play: 

a)  Performance  vs.  incident  beam  intensity 

b)  Noise  generation 

c)  Asymmetric  thrust  considerations 

d)  Scaling  laws  for  various  blast  wave  geometries 

e)  Absorption-wave/thruster  geometry 

f)  Passive  refresh  of  planar  vs.  line  sources 

g)  Wavelength  scaling  effects  (LSC/MSP  regime) 

h)  Ignition  requirements  for  radiation-sustained  plasmas 
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Figure  11*14  RdiriowMp  of  ladlfUMl  Hm  mpm*  tkmt  victor  aad  nUcb 
C.G.  (eeater  at  gravity) 


Figure  11*15  T bruit  rector  control  reUtlonihipa 


i)  Atmospheric  smosoi/pfikulatc  contamination 


a)  Performance  v*.  Incident  beam  Intenrity 

The  ERH  thruster  is  a  repetively-pulsed  engine  which  produces  thrust  by  expanding  blast 
waves  over  an  Impulse-coup llnf  surface  especially  desifned  for  that  role.  These  blast  waves  are,  of 
course,  created  by  pulsed  radiation-sustained  plasmas  ignited  below  the  reaction  surface.  As  will  be 
demonstrated  in  a  later  chapter,  maximum  coupling  coefficients  (Newtons  of  thrust  per  megawatt  of 
beam  power)  are  realised  for  the  greatest  beam  intensities  delivered  to  the  propagating  absorption 
waves.  Whyt  Because  the  highest  intensities  create  the  largest  blast  wave  overpressures. 

Clearly,  this  relationship  has  a  strong  wavelength  dependence,  which  also  limits  the  clear  air 
electrical  breakdown  threshold.  As  indicated  in  Fig.  11-16,  the  highest  beam  intensities  are  accessible 
only  for  the  shortest  wavelength*.  The  shortest  wavelength  of  interest  is  0.3  microns,  for  two  reasons. 
Fust,  it  is  difficult  to  construct  a  mirror  which  will  be  highly  reflective  in  the  ultraviolet  regime. 
Second,  much  below  0.35  microns,  atmospheric  (molecular  and  particulate)  scattering  loesee  become 
intolerable.  Hence,  one  must  confine  ERH  thruster  operation  to  “atmospheric  window”  wavelengths 
between  0.3  microns  and  10  cm. 

This  is  not  meant  to  imply  that  lower  beam  intensities  and  lower  blast  wave  overpressures 
typical  of  the  LSC  wave  regime  are  not  of  real  interest.  Quite  to  the  contrary,  this  regime  may  be 
ideally  suited  to  low  vehicle  accelerations  and  VTOL  flight  with  reduced  noise  levels.  It  is  apparent 
though,  that  the  highest  vehicle  accelerations  are  attainable  only  with  the  shortest  practical  beam 
wavelengths. 

b)  Noise  generation 

Since  the  greateet  blast  wave  overpressures  yield  the  highest  coupling  coefficients,  much  at¬ 
tention  should  be  given  to  mechanisms  for  noise  supresslon/attenuatioa/alleviation-  especially  when 
addrearing  human  payloads.  The  problem  of  noise  generation  by  an  ERH  thruster  goes  beyond  con¬ 
siderations  of  environmental  impact  to  the  immediate  external  surroundings.  One  must  also  consider 
health  and  safety  issues,  especially  with  regard  to  the  crew  who  mus*  reside  within  close  proximity 
to  an  extremely  powerful  noise  source. 

We  should  aote,  here,  that  the  human  ear  is  sensitive  to  regular  compression  waves  in  the 
range  of  40Ht  to  20  kHt;  and  for  young  people,  16  Ht  to  20  kHe.  In  contrast,  dogs  can  hear  up  to 
30  kHe  -  40  kHe  and  bats  up  to  100  kHt.  Regular  compression  waves  at  frequencies  below  those 
which  we  perceive  u  sound  can  often  be  perceived  as  vibration:  one  can  still  sense  frequencies 
lower  thin  those  which  we  consider  ‘heiring’.  Ordinarily,  one  can  distinguish  three  characteristic 
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audible  sound*  pitch,  vohuas  (tetsokty),  and  ton*  quality.  For  example,  regular  wnusoidal  varying 
compress ion  waves  will  hav«  a  different  tone  quality  than  a  sawtooth  shaped  wave. 

Figure  n-17  Illustrates  the  normal  equal  loudness  contours  for  pure  tones.  Note  that  if 
ERH  thruster  operation  is  restricted  to  either  the  isfrarsonic  regime  (for  low  acceleration,  PRF 
*  40  kHs),  or  the  ultra-sonic  regime  (for  maximum  acceleration,  PRF  a  15  kHs),  the  apparent 
Impact  on  human  hearing  can  be  minimised.  Incidentally,  when  operating  an  ERH  thruster  within 
the  Lnfra-sonic  regime,  one  should  avoid  the  potentially  damaging  resonant  frequencies  of  the  heart, 
lungs  and  nerve  synapses. 

Another  obvious  way  to  reduce  thruster  noise  is  to  minimis  energy  delivered  in  a  single  laser 
puke.  The  required  thrust  level  and  time-averaged  power  can  be  maintained  amply  by  increasing 
the  frequency  of  laser  pulses.  This  approach  would  result  in  LSD-  wave-heating  of  the  smallest 
possible  volume  of  air,  with  the  minimum  feasible  laser  puli'  energy;  then  this  “detonation*  source 
would  be  rotated  sequentially  about  the  thruster  surface.  More  detaik  on  thk  concept  will  be  given 
below  in  another  section  devoted  to  the  “Rotating  Detonation  Wave  Engine*  (RDWE). 

Inddeotly,  the  reference  ERH  thruster  energy  deposition  mode  of  12  LSD  waves  (puked 
■Htmanushr  with  a  total  beam  energy  of  approximately  1.6  MJ)  was  chosen  without  regard  to  noise 
generation.  A  more  ideal  “source”  would  utilise  48  smaller  LSD  waves,  each  one  puked  individually 
(with  an  energy  of  about  34  kJ  each)  to  create  a  rotating  One  source.  Thk  would  result  in  a  pulsing 
frequency  48  times  greater  than  the  reference  design.  When  producing  maximum  thrust  with  a 
time-averaged  beam  power  of  2  GW,  the  ideal  thruster  could  accelerate  a  5666  kg  spacecraft  at  30 
G’s,  nft»— >*— jy  tt  an  ultra-eonic  frequency  of  60,000  Ha.  In  hover  (above  “ground  effect”  altitude), 
when  all  48  waves  are  puked  simultaneously  at  40  Hl  the  thruster  will  operate  almost  silently  within 
the  infra- sonic  regime. 

c)  Asymmetric  thrust  considerations 

During  low  PRF  hover  operation,  caution  must  be  exercised  to  avoid  unstable  rotation  rates 
for  ofigl  thrust  vectors,  of  course.  This  situation  could  be  alleviated  by  using  symmetric  thrust 
geometries  of  two,  three  or  four  “lift  points”,  which  rotate  over  the  thruster  surface  (at  the  cost  of 
increased  noise  generation).  Figure  11-18  shows  these  options. 

d)  Scaling  laws  for  particular  blast  wave  geometries 

Tabit  ffl-1  contains  a  listing  of  the  various  available  blast  wave  geometries  and  applicable 
scaling  laws  (e.g.,  planar,  powered  planar,  cylindrical,  powered  cylindrical,  spherical  and  powered 
spherical).  First  order  estimates  of  the  impulse  coupling  efficiency  for  a  given  ERH  thruster  model 


Normal  equal  loudness  contours  lor  pure  tones.  They  can  be  applied 

when: 

a)  The  source  ol  sound  Is  directly  ahead  ol  the  listener. 

b)  The  sound  reaches  the  listener  in  the  form  ol  a  free  progressive  plane 
wave. 

c)  The  sound  pressure  level  is  measured  In  the  absence  of  the  listener. 

d)  The  listening  Is  binaural. 

e)  The  listeners  are  otologically  normal  persons  in  the  age  group  18  to  25 
years  inclusive. 


Figure  n-17  Normal  equal  loudneee  contour*  for  pure  tone*  (from  Ref. 55) 
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Table  n-l-Sedov’i  Scaling  Laws,  (See  Sedov  (22|) 
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Figure  n-18  ERH  thrust  vector  geometry  for  rotating  line  source(s) 
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can  b*  derived  by  application  of  th«aa  sealing  law*,  one*  th*  detailed  Source*  geometry  and  energy 
deposition  mods  have  been  established. 

e)  Absorption-wave/thruster  geometry 

Figure  11-19  displays  seven  representative  energy  deposition  modes  for  radially-symmetric 
thrusters  utilizing  truncated  conical  plugs  (like  the  point  design).  Variations-on-a-theme  for  the 
line/cvllndrical  source  are  presented  in  Flgi.  0-19a  and  II- 19b;  whereas,  similar  modes  for  the 
nlanar  source  are  given  In  Figs.  Q-19c  through  U-19g.  These  illustrations  present  bottom  views  of 
the  ERH  thruster  plug  with  the  secondary  optics  (SO)  wrapped  about  the  upper  plug  circumference. 
Ignition  of  the  radiation-  sustained  plasma  wave  is  assumed  to  take  place  around  the  lower  end  of 
the  truncated  plug.  For  the  “line  source*  cases,  the  SO's  are  assumed  to  project  high  intensity  power 
beams  towards  a  point  focus;  for  the  “planar  source*,  a  line  focus. 

All  the  modes  displayed  in  Fig.  Q-19  assume  the  use  of  rapidly-propagating  (i.e.,  hypersonic 
velocity)  LSD  waves,  probably  falling  in  the  “super  detonation  wave  regime*  (according  to  the  no¬ 
menclature  of  Raizer  (2|).  The  plasma  expansion  time  is  generally  a  factor  of  10  longer  than  the 
laser  pulse  duration  (Le.,  the  beam  energy  “deposition*  time). 

The  following  describes  the  process  of  pulsed  energy  deposition  for  each  mode  in  Fig.  D- 
19.  (Note  that  the  basic  principles  behind  the  “simultaneous-line”  source  and  the  “rotating-line” 
source  were  already  covered  in  Fig.  11-18.)  Two  choices  for  PO/SO/TO  optical  train  geometries 
exist:  i)annular  PO  with  a  scalloped  annular  SO  linked  to  a  set  of  small  grating-  incidence  TO’s, 
or,  ii)acalloped  PO  and  a  set  of  ^Tular  SO’s.  For  *  simultaneous-line  source,  the  entire  PO  must 
be  illuminated  at  once,  with  a  ar  beam  (i.t  having  a  hollow  center).  For  the  rotating-line 
source,  discrete  portions  of  the  exposed.  This  is  accomplished  with  a  single  (smaller)  circular 
beam,  or  with  multiple  circular  -ms,  distributed  around  the  annular  PO  surface.  After  delivery 
of  a  laser  pulse,  the  power  beam  must  then  be  slewed  to  the  next  “targeted”  location  on  the  PO,  as 
indicated  earlier  in  Fig.  II- 7. 

The  other  modes  portrayed  in  Figs.  Q-19c  through  Q-19g  would  use  continuous  annular  PO 
and  SO  optical  trains.  The  Fig.  II-19c  mode  pulses  the  whole  engine  lower  surface  at  once;  hence, 
the  entire  plug  surface  is  “active”  and  the  only  unilluminated  surface  results  from  shadows  cast  by 
the  shroud  support  struts. 

Figures  D-19d  and  D-19e  show  the  concept  of  a  “rotating-wedge”  planar  source,  an  alterna¬ 
tive  mode  obtained  by  illuminating  some  smaller  portion  of  the  PO  (i.e.,  during  a  given  laser  pulse). 
A  sequence  of  pulses  is  delivered  to  adjacent  segments  of  the  annular  PO  “race  track”.  Figure  II-19f 
indicates  a  variation  on  the  Fig.  II-19d  mode,  wherein  the  laser  beam  placement  on  the  PO  “race 
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I 

C)  SIMULTANEOUS  PLANAR  D)  ROTATING  "WEDGE"  PLANAR 

I 


Figure  D-19  Energy  deposition  mode*  -  bottom  view  (un  powered  cylindrical 
and  planar  aourcea) 
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E)  ROTATING  "MULTI -WEDGE"  F)  ROTATING  SPIRAL  WEDGE 


G)  ROTATING  SPIRAL  MULTI-WEDGE 


Figure  H-19  Energy  deposition  mode#  -  bottom  view  (unpowered  cylindrical 
and  planar  source*  ),  cont’d. 


track*  u  «Uwd  during  the  Iiht  pula*  -  yi aiding  %  apirmi-ahapad  plaama  wedge.  Tha  apiral  moda  may 
enable  a  reduction  in  the  time-avenged  heat  load  sustained  by  an  annular  SO  mirror,  by  spreading 
the  incident  beam  ‘foot-print’  over  a  larger  SO  surface  area  (Le.,  during  a  given  laser  pulse).  One 
additional  advantage  may  occur  from  simplifying  the  pointing/tracidng  requirements  for  the  laser 
transmitter.  The  transmitter  antenna  can  be  sieved  constantly  in  a  circular  pattern  about  the  PO 
“race  track”  without  having  to  start  and  stop  (Le.,  re-target)  between  laser  pulses. 

One  additional  comment  should  be  made  with  regard  to  the  "rotating-line*  energy  deposi¬ 
tion  mode  of  Fig.  U-  19b.  If  pulsed  at  sufficiently  frequent  intervals,  the  thruster  could  operate  aa  a 
continuous  Rotating  Detonation  Wave  Engine  (RDWE)  mode  as  shown  in  Fig.  D-20.  In  this  mode, 
frmh  sir  which  resides  just  ahead  of  the  advancing  blast  wavefront  would  be  “detonated"  just  prior 
to  the  arrival  of  the  front.  Hence,  In  a  time-averaged  sense,  the  energy  deposition  process  would 
simulate  that  of  a  detonating  line  source  that  rapidly  wipes  around  the  thruster  plug. 

f )  Passive  refresh  of  planar  vs  line  sources 

One  final  point  should  be  made  with  regard  to  the  “lateral  refresh”  mechanism  and  its 
suitability  for  use  with  line  (finger)  source  and  planar  source  geometries.  The  line  source  seems 
Ideally  suited  to  this  refresh  mechanism.  However,  the  planar  source  does  not  appear  aa  appropriate, 
especially  when  the  entire  lower  surface  is  pulsed  (see  Fig.  H-19c).  Why?  Because  with  the  planar 
source,  the  ratio  of  active  to  passive  thruster  surface  is  so  small  (Le.,  roughly  equal  to  1),  that  the 
most  expedient  “refresh"  must  take  place  at  right  angles  to  the  thruster  surface  -  as  one  would 
expect  with  a  1-D  planar  source.  As  a  result,  the  maximum  frequency  at  which  planar  sources  can 
be  pulsed  will  be  different  than  that  for  multi-line  sources.  The  implications  of  this  on  ERH  thruster 
performance  has  yet  to  be  examined  in  any  detail. 

g)  Wavelength  scaling  effects  (LSC/MSP  regime) 

Up  to  now,  energy  deposition  modes  suitable  for  high  intensity  LSD-waves  have  been  em¬ 
phasized  largely  because  of  their  greater  impulse  coupling  efficiency.  However,  thrust  can  also  be 
generated  using  lower  intensities  typical  of  the  LSC  wave  regime  (albeit  with  reduced  efficiency),  if 
the  plasma  is  still  ignited  within  close  proximity  to  the  impulse  surface.  Figure  11-21  indicates  two 
feasible  energy  deposition  modes:  i)  powered  cylindrical,  and  ii)  powered  spherical.  In  the  first  case 
(Fig.  Q-21a),  an  annular  power  beam  is  projected  by  a  ring-type  SO  mirror  towards  a  liflsJgqia  - 
which  is  oriented  perpendicular  to  the  thruster  surface.  (This  thruster  mode  is  analyzed  in  Chapt. 
VII.)  Figure  H-21b  shows  a  variation  of  this  theme,  wherein  the  annular  power  beam  is  again  pro¬ 
jected  by  a  “ring-type"  SO  mirror  towards  a  central  focal  region  -  but  this  time  a  point  focus  is 
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uiumed,  not  *  lino.  The  resultant  optical  gaometry  con  perhaps  crooto  o  powered  aphoricol  bloat 
wove  below  the  ERH  thruster  surface. 

For  the  geometries  displayed  in  Fig.  Q-21,  microwave  radiation  could  have  a  significant 
advantage;  it  can  produce  plasma  temperatures  which  are  much  lower  than  those  generated  by 
optical  frequencies.  Whereas  electron  temperatures  will  still  be  10000  to  15000  K,  ions  and  neutrals 
can  be  as  low  as  2000  to  4000  K  -  so  the  resultant  air  plasmas  can  be  in  a  non-equilibrium  condition. 
In  contrast,  LSG  wave  plasmas  will  produce  equilibrium  temperatures  of  10,000  -  30,000  K.  Since 
radiant  heat  transfer  to  the  ERH  thruster  is  proportional  to  oT4,  a  factor  of  3  reduction  in  plasma 
temperature  (T)  can  lower  overall  heat  transfer  rates  by  several  orders  of  magnitude.  Subsequently, 
the  thermal  stress  problem  (for  ERH  thruster  materials)  should  be  greatly  alleviated  by  using 
microwave-sustained  plasmas  for  low  acceleration,  “pop-up*  maneuvers  and  hover  flight  modes. 
During  these  times,  the  plug  surface  would  be  exposed  to  potential  “overheat*  conditions  caused  by 
ingestion  (i.e.,  undesired  recirculation)  of  hot  exhaust  gases. 

It  should  be  noted  that  microwave-induced  electrical  air  breakdown  normally  occurs  before 
detonation  waves  can  be  produced:  e.g.,  flux  densities  of  1  •  2  MW /cm?  at  sea  level  pressure.  How¬ 
ever,  radiation-sustained  plasmas  typical  of  the  “deflagration”  regime  are  easily  ignited.  Hence,  all 
evidence  points  favorably  to  the  potential  for  thrust  production  with  repetitively-pulsed  microwave 
power.  Since  inverse  bremastrahlung  absorption  coefficients  for  microwave  radiation  are  exceedingly 
low,  the  physical  size  of  powered  cylindrical  and  spherical  plasmas  (shown  in  Fig.  H-20)  can  grow  to 
be  much  larger  than  with  laser  wavelengths.  Significant  blast  wave  overpressures  (10  atm  or  more) 
could  be  generated  over  a  sizeable  portion  of  the  ERH  thruster  surface. 

h)  Ignition  requirements  for  radiation-sustained  plasmas 

In  gll  the  above  mentioned  energy  deposition  modes,  inverse  bremsstrahiung  absorption 
is  invoked  as  the  radiation  absorption  mechanism.  However,  a  “breakdown”  mechanism  is  needed 
to  start  the  cascade  process  in  which  enough  free  electrons  are  available  to  ignite  an  LSC  wave, 
LSD  wave,  or  Microwave-Supported  Plasma  (MSP)  wave.  Several  breakdown  mechanisms  exist:  i) 
laser-induced  clean  air  breakdown,  ii)  particulate  “seeding”,  iii)  aerosol  or  molecular  “seeding”,  iv) 
electric  spark,  v)  metallic  target,  and  vi)  thermionic  or  field  emission.  Options  ii)  and  v)  present 
potential  optical  fouling  problems.  Option  iii)  is  limited  to  absorption  over  a  finite  wavelength  range. 
Options  iv)  and  vi)  require  an  on-board  electric  power  source;  of  these,  option  iv  would  probably 
consume  the  least  power. 

Overall,  option  i)  is  the  preferred  breakdown  mechanism.  However,  if  the  PO/SO  optical 
focusing  geometry  is  not  able  to  produce  sufficient  intensity  to  induce  electrical  breakdown,  an  elec- 
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A.  POWERED  CYLINDRICAL 


B.  POWERED  SPHERICAL  ^ 

LASER  laser  beam  has  spherical 


Figure  D-  21  ERH  thruster  energy  deposition  modes  -  bottom  view  (powered  cylindrical 
and  powered  spherical  modes) 
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trie  (park  device  would  bo  4  useful  alternative.  A s  mentioned  earlier,  several  strategically-located 
"spark  plugs*  distributed  along  the  "detonation-line*  sources  would  have  the  valuable  ability  to  trig¬ 
ger  breakdown  at  numerous  locations;  this  mechanism  should  enable  a  thrust-vectoring  capability. 

1)  Atmospheric  aerosol/particulate  contamination 

Several  final  comments  should  be  made  regarding  the  data  presented  in  Pig.  11-22.  Note 
that  the  clean  air  breakdown  threshold  is  beyond  VfiW/cm2.  For  the  reference  design  which  uses  a 
uniform  intensity  of  S  x  1 O^W/cm3,  breakdown  can  be  triggered  upon  a  metal  target  in  roughly  20 
nanoseconds  -  with  an  energy  cost  of  10  Jjcm2  integrated  over  the  beam  cross-section.  Ignition  on 
contaminating  particles  (0.3  micron  diameter)  would  be  triggered  in  100  nanoseconds,  or  an  energy 
cost  of  6X  larger.  Clean  air  breakdown  cannot  be  triggered  at  all  for  the  reference  system  Intensity 
of  5  x  10 *W  I  cm2.  Hence,  the  metallic  igniter  plug  or  the  electric  spark  would  be  the  preferred 
ignition  mechanism. 

Let  us  now  examine  another  optical  train  and  ERH  thruster  design  that  incorporates  a  multi¬ 
point  focus  SO  (Le.,  tapered  plasma  "fingers”),  and  uses  the  laser-induced  breakdown  mechanism. 
With  5x  10*W/cm?  at  the  ignition  point,  it  is  evident  that  roughly  40  nanoseconds  must  be  invested 
at  an  energy  cost  of  200  J/cm2  -  but  of  course,  the  focused  beam  area  (at  the  Ignition  point)  would 
be  reduced  very  much  below  the  afore-mentioned  example.  However,  as  seen  in  Fig.  H-22,  ignition 
on  contaminating  particles  (which  accidently  enter  the  beam)  could  occur  at  even  earlier  times  -  and 
perhaps  at  an  undesired  location  along  the  ERH  thruster  plug,  with  potentially  hazardous  results. 
This  may  cause  a  severe  vectoring  of  the  thrust  in  an  unpredictable  fashion.  Hence,  the  super  high 
thrust/acceleration  abilities  of  the  ERH  thruster  should  probably  not  be  invoked  immediately  at 
liftoff.  The  certainty  of  encountering  particulate  contamination  near  the  ground  may  necessitate  low 
acceleration  pop-up  maneuvers,  wherein  the  engine/vehicle  climbs  to  a  safe  altitude  (perhaps  as  low 
as  tree-top  level)  -  before  really  "blasting  off*  at  high  rates  of  acceleration. 

11.  Air  "Refresh*  Schemes 


During  normal  repetitively-pulsed  ERH  thruster  operation,  there  comes  a  time  when  the 
heated,  high  pressure  air  resulting  from  a  given  radiation-sustained  plasma  ("detonation"  or  “defla¬ 
gration"  )  wave  eventually  decays  to  ambient  pressure.  This  volume  of  hot,  ambient-pressure,  “spent” 
air  must  be  converted  out  of  the  actively-irradiated  thruiter  region  -  before  it  can  be  "pulsed”  again. 
The  entire  cyclic  procem  can  be  likened  to  the  intake/compression/ignition/exhaust  mechanisms  in 
conventional  4-cycle  internal  combustion  engines.  Clearly,  if  the  time  required  to  resupply  the  "deto- 
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THRESHOLD  INTENSITY  (W/cm 


PULSE  TIME -TIME  TO  BREAKDOWN  (sec) 


Figure  IT- 22  Measurements  of  air  breakdown  in  the  pretence  of  particulate 
matter,  (from  Reft.  37  •  39) 
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nation  sono*  can  somehow  ba  reduced,  than  U  will  ba  poaaibla  to  incraaaa  tha  angina  puiaa  rapatition 
frequency  and  hence,  time-averaged  throat. 

a)  Baaic  modes  for  refreahment 

It  is  useful  to  examine  the  mechanisms  by  which  this  ‘refresh”  process  might  be  enhanced. 
The  reference  platform  utilises  an  LSD  wave-heated  (line/cylindrical)  source,  so  it  will  serve  as  the 
exemplary  ‘straw-man*  ERH  thruster  in  this  section.  Note  in  Pig.  11-23  that  there  are  only  two 
baaic  modes  for  “refreshment”  of  the  ERH  thruster  surface:  i)  passive,  or  ii)  active.  In  the  former 
“passive”  mode,  sufficient  time  must  be  allowed  for  the  heated  cylindrical  hot  plasma  “bubble”  to 
“burst”,  then  collapse  back  upon  itself  -  allowing  cool  ambient  air  to  rush  back  into  the  “detonation 
tone".  For  the  latter  “active*  mode,  some  type  of  forced  convection  mechanism  is  employed  to  drive 
fresh  air  back  into  the  “detonation  tone”  -  often  in  leas  time  than  that  required  for  a  complete 
“passive  refresh*  of  the  surface. 

b)  Passive  refresh  enhancement 

The  various  mechanisms  available  for  “passive*  refresh  enhancement  can  be  incorporated 
Into  a  given  thruster  simply  by  making  judicious  "design  decisions"  in  the  overall  engine  geometry. 
A  number  of  these  principles  were  first  introduced  earlier  in  this  chapter.  For  example,  the  total 
time  required  for  refresh  of  the  reference  ERH  thruster  (with  its  twelve  line  source  and  truncated 
conical  plug)  is  the  sum  total  ofr  1)  LSD  wave  ignition  time,  ii)  time  for  LSD  waves  to  propagate 
across  the  entire  plug  length,  iii)  time  for  blast  waves  to  expand  laterally  to  the  local  static  pressure, 
and  iv)  time  for  plasma  bubbles  to  cool  and  collapse  back,  so  that  fresh  air  can  rush  in.  Clearly,  the 
first  two  time  intervals  can  be  minimized  by  selecting  the  maximum  beam  intensity  (within  some 
safety  margin)  that  can  be  projected  upon  the  propagating  LSD  waves  by  the  SO  mirrors;  this  in 
turn,  minimizes  the  necessary  laser  pulse  duration. 

The  third  time  interval  can  be  minimized  in  two  ways.  First,  for  a  given  incident  beam 
intensity  (upon  the  LSD  wave),  one  can  minimize  the  blast  wave  expansion  time  by  selecting  the 
smallest  “detonation”  volume.  For  a  given  thruster  plug,  one  would  then  select  the  minimum  laser 
beam  diameter  consistent  with  the  limits  of  diffraction,  turbulence,  thermal  blooming  and  jitter 
of  the  SO  (or  TO)  mirrors.  This  small  diameter,  high  pressure  cylindrical  volume  of  plasma  can 
certainly  relax  more  quickly  to  the  local  static  pressure  than  a  larger  volume.  Secondly,  there 
are  ways  to  immerse  the  entire  thruster  plug  (and  its  detonation  cones)  in  static  pressure  levels 
significantly  above  ambient;  hence,  lees  time  elapses  during  the  plasma  expansion  process.  Finally, 
as  mentioned  earlier  for  the  multi-fingered  reference  thruster,  care  must  be  taken  to  minimize  the 
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overlapping  of  adjacent  cylindrical  blast  war*.  Thi*  condition  will  hinder  tho  subsequent  passive 
refresh  process,  u  well  as  optimum  thrust  production. 

Similar  comments  can  be  made  with  refard  to  the  fourth  time  interval,  which  is  proportional 
to  th#  maximum  blast  wavs  expansion  radius  (i.s.,  when  plasma  pressure  first  equals  local  static), 
and  the  sonic  velocity.  Since  cool  air  can  rush  in  only  at  the  speed  of  sound,  one  would  want  to 
minimize  the  expanded  plasma  volume,  and  maximize  the  local  sound  speed.  The  local  sound  speed 
can  be  increased  to  some  extent,  with  the  aid  of  compressed  cowl  air,  as  will  be  discussed  later. 

c)  Active  refresh  enhancement 

Figure  0-24  displays  eight  schemes  for  “active"  refresh  enhancement.  Over  half  have  already 
been  introduced.  In  Fig.  Q-24a,  the  engine/vehicle  is  shown  in  supersonic  axial  flight,  utilizing  the 
mechanism  of  forced  axial  convection  for  “active  refresh*  of  the  ERH  thruster  plug.  The  degree  to 
which  fresh  air  can  be  convected  across  the  plug  surface  is  dictated  by  the  specifics  of  engine/vehicle 
configuration:  i.e.,  inlet  forebody  shape,  shroud  contour  (if  not  shroudless),  and  nozzle/afterbody 
geometry.  As  will  be  analyzed  more  thoroughly  in  Chapt.  VI,  a  properly  designed  cowl  and  inlet 
forebody  enables  the  plug  to  be  “immersed”  in  axial-refresh  air  with  a  local  pressure  (and  sound 
speed)  quite  above  ambient. 

Depending  on  the  inlet  desip  Mach  number  (Afo)  and  Bight  speed,  airflow  emerging  from 
the  annular  cowl  duct  can  either  be  subsonic  or  supersonic.  If  ducted  air  expands  at  supersonic 
velocity  across  the  plug  nozzle,  it  can  be  most  effective  in  the  “active”  refresh  process.  Even  so, 
during  the  time  interval  required  to  clear  the  entire  length  of  the  ERH  thruster  plug  of  hot  spent 
air,  several  new  detonations  will  have  been  triggered.  However,  note  that  blast  waves  ignited  at 
the  lower  end  of  the  plug  have  the  most  time  to  expand  (i.e.,  as  the  LSD  wave  races  up  the  plug 
surface),  so  the  plug  aft  end  will  be  refreshed  first,  passively.  Hence,  it  is  the  upper  end  that  needs 
the  active  refresh  enhancement  most  urgently,  and  gets  it.  Incidently,  the  fully  expanded  plasma 
“fingers*  should  not  obstruct  a  significant  fraction  of  the  internal  cowl  cross-section  -  so  as  to  cause 
a  reduction  in  the  mass  flow  rate  of  axial  refresh  sir. 

Figure  D-24b  shows  another  forced  convection  refresh  mode  that  could  occur  in  subsonic 
lateral  flight.  In  this  esse,  the  local  wind  vector  would  “blow”  at  rigbt  angles  to  the  plasma  fingers  - 
helping  to  clear  out  spent  hot  air.  In  subsonic  flight,  this  mode  may  not  accomplish  anything  more 
significant  than  to  prevent  the  local  static  temperatures  (of  “passive  refresh”  air)  from  becoming 
too  large. 

d)  “RDWE*  active  refresh  mechanism 
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A.  SUPERSONIC  AXIAL  FLIGHT 
(FORCED  AXIAL  "REFRESH") 


B.  SUBSUN  I C/SUPERSUN  I C  LATERAL  FLIGHT 
(FORCED  LATERAL  "REFRESH) 


HORIZONTAL 
FLIGHT  DIRECTION 


Figure  11-24  Mechanisms  for  active  refresh  enhancement 
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The  active  refresh  mechanism  displayed  in  Fig  II-24c,  may  enhance  the  parformanea  of  Ro¬ 
tating  Detonation  Wave  Engines  (RDWE)  using  a  symmetric,  multi-line  ERH  thruster  geometry. 
Beneficial  interaction  is  anticipated  between  two  adjacent,  sequentially-pulsed  sources.  One  possi¬ 
bility  is  that  the  pre-cursor  cylindrical  shock  wave  shed  by  a  given  detonating  One  could  convect 
“refresh*  gases  laterally  into  the  rite  of  the  next  “line*  to  be  pulsed.  It  might  even  able  a  pre- 
pressuriring  effect.  A  second  possibility  Involves  a  similar  beneficial  exchange  between  adjacent 
detonations,  but  only  for  the  flai  and  iul  line  sources  -  as  the  RDWE  “overlaps”  after  a  complete 
360*  sweep  of  the  ERH  thruster  plug.  A  third  possibility  involves  beneficial  interaction  between  the 
trailing  jet  behind  one  spiral-source  LSD  wave,  and  the  adjacent  spiral  detonation  cone  (e.g.,  see 
the  engine  geometry  in  Fig.  U-24g).  These  opportunities  for  refresh  enhancement,  have  yet  to  be 
explored  in  any  depth. 

e)  “porous  plate”  active  refresh  mechanism 

Figure  H-24d  portrays  another  scheme  for  active  refresh  enhancement,  one  which  forcibly 
convects  “refresh"  gases  in  a  direction  perpendicular  to  the  ERH  surface.  This  concept  assumes  that 
the  thruster  “plate*  is  porous  and  that  the  refresh  fluid  is  pumped  out  with  an  on-board  auxiliary 
power  supply.  The  refresh  “fluid"  would  be  stored  internally  in  a  liquid  state,  then  evaporate 
through  the  ERH  thruster  as  it  emerges  from  the  gaseous  plate.  (Note  that  a  similar  function  could 
be  provided  by  an  ablatively-cooled  ERH  thruster  surface.) 

Two  technical  problems  are  evident  for  such  a  scheme.  First  is  the  mass  loss  due  to  refresh- 
fluids,  which  would  penalise  the  ERH  thruster  with  a  finite  specific  impulse.  Clearly,  one  would 
want  to  minimize  both  the  refresh-fluid  mass  flow  rate  and  the  total  porous  plate  surface  area. 
Perhaps  the  only  porous  regions  could  be  restricted  to  just  the  actively-irradiated  plug  nozzle  areas 
(i.e.,  immediately  below  the  line  sources).  The  second  problem  area  relates  to  the  pressure  ratios 
which  must  be  applied  across  the  porous  plates,  remembering  that  peak  blast  wave  pressures  reach 
600  atm.  in  the  ref.  platform.  Also,  it  is  not  obvious  that  the  refresh  gases  should  be  admitted 
continuously:  perhaps  a  repetitively-pulsed  mode  would  give  better  performance.  If  the  former  is 
chosen,  the  porous  plate  surface  might  act  as  a  self-operated  acoustic  “valve* ,  anyway. 

One  final  note  should  be  made  regarding  tbe  porous  plate  option.  This  mechanism  might 
permit  the  thruster  surface  to  perform  as  a  rocket  within  a  near  vacuum  environment.  However, 
the  gaseous  working  fluid  must  be  carefully  expelled  and  uniformly  heated  to  a  high  temperature 
before  the  rocket  expansion  process  is  completed.  This  sequence  of  events  might  present  a  technical 
problem  with  the  rocket-mode  scheme. 
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c.  ROTATING  DETONATION  WAVE  ENGINE  (RDWE) 

(BLAST  WAVES  CONVECT  "REFRESH"  GASES  LATERALLY) 
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Figure  II- 24  Mechanim*  for  active  refreeh  enhancement,  cont’d. 
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f)  “Rotating  plate"  active  refresh  mechanism 

Illustrated  to  Fig.  Q-24e  Is  smother  concept  for  active  refresh  of  an  ERH  thruster,  which 
invokes  a  rapid  rotation  of  the  radially-symmetric  plug  nozzle.  In  this  scheme,  the  thruster  surface 
is  spun  up  to  a  maximum  velocity  of  Mach  I  (or  2)  at  the  outermost  rim  of  the  conical  plug;  high 
subsonic  velocities  might  also  be  usefuL  The  goal  would  be  to  “drag”  local  refresh  air  laterally 
through  the  detonation-  zones  (line  source  sites),  by  means  of  viscous  forces  within  the  turbulent 
boundary  layer.  As  shown  in  Fig.  Q-24e,  the  detonation  zones  must  be  placed  extremely  close  to 
the  thruster  surface  anyway,  in  order  to  maximize  impulse  coupling.  Hence,  line  sources  would  be 
located  completely  inside  the  boundary  layer,  probably  within  one  centimeter  (or  leas)  from  the 
ERH  thruster  surface.  It  is  evident  that  the  resultant  blast  wave  would  cause  separation  of  the 
boundary  layer,  and  sufficient  time  must  be  allowed  for  it  to  re-attaeh  before  the  thruster  could  be 
pulsed  again. 

This  refresh  mechanism  could  possibly  be  enhanced  by  increasing  surface  roughness,  and 
perhaps  a  well-used  ablatively-cooled  re-entry  heat  shield  might  just  make  the  perfect  ERH  thruster 
surface.  Exaggerated  surface  roughness  will  certainly  produce  additional  aerodynamic  drag,  and 
require  more  power  for  “spin  up".  However,  this  may  be  worth  the  trouble. 

Azimuthal  flow  circulation  might  be  further  enhanced  by  the  attachment  of  “flow  fences” 
to  the  ERH  thruster  surface  (as  indicated  in  Fig.  H-24f).  Flow  fences  would  be  high  enough  to 
efficiently  pump  fresh  air  into  the  detonation  zones,  yet  not  so  high  as  to  trigger  separation  of 
the  boundary  layer.  With  such  devices  attached,  plug  nozzle  rotation  could  generate  significant 
centrifugal  forces  within  the  immediate  boundary  layer  air  flows,  as  shown  in  Fig.  24g.  The  effect  of 
these  forces  upon  the  resultant  spiraling  air  flow  patterns  and  the  refresh  process  should  eventually 
be  examined  in  subsequent  research. 

All  the  above  mentioned  “active  refresh”  mechanisms  would  require  the  ERH  thruster  plug 
to  be  spun  up  slowly  on  the  ground,  just  prior  to  launch,  using  an  auxiliary  on-  board  power 
supply.  Following  liftoff,  no  on-board  power  would  be  needed  during  the  20  -  30  second  boost  period 
when  the  ERH  thruster  is  accelerated  from  Mach  0  to  3.  The  massive  rotational  inertia  (stored  in 
the  engine/vehicle)  would  be  gradually  be  released  and  dissipated  by  aerodynamic  drag,  along  the 
orbital  trajectory. 

g)  “Axial  fan-driven”  active  refresh  mechanism 

Figure  U-24h  portrays  another  active  refresh  mechanism  which  invokes  rotating  components. 
In  this  engine/vehicle  configuration,  a  supersonic  axial  (or  perhaps  centrifugal)  fan  is  housed  within 
the  annular  shroud  volume.  In  preparation  for  liftoff,  a  small  auxiliary  power  supply  would  (again) 
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Figure  11-24  MechanUma  for  active  refreeh  enhancement,  cont’d. 


59 


G,  CENTRIFUGAL  FORCES  ON  BOUNDARY  LAYER  FLOWS  (HOVER  CONDITION) 


SUEJia 


BOTTOM  VIEW 


H,  AXIAL  FAN-DRIVEN  "REFRESH"  AIR: 


AXIAL  FAN 
HOUSED  WITHIN 
SHROUD 


Figure  II-24  Mechanic  mi  for  active  refreih  eohancement,  cont’d. 


•pin  up  th*  *ngin*,  storing  much  rotational  »n«rgy  in  th«  available  mum  -  whil*  «x«rting  a  reaction 
force  against  the  ground.  During  the  spin-up  process,  the  variable-pitch  fan  would  be  set  at  zero 
pitch  to  minimize  input  power  requirements  and  aerodynamic  drag. 

At  liftoff,  the  pitch  would  suddenly  be  increased,  accelerating  refresh  air  across  the  ERH 
thruster  plug.  Efficient  axial  refresh  would  probably  not  occur  until  the  ducted  air  velocity  turned 
supersonic.  Even  though  the  rate  of  decrease  in  engine/vehicle  rotation  would  significantly  exceed 
that  for  the  “rotating  plate”  mechanism  in  Fig.  II-24e,  it  still  might  provide  sufficient  axial  refresh 
flow  for  the  20  -  30  seconds  duration  of  ERH  thruster  boost.  (After  Mach  3  or  so,  the  fan  blades  would 
be  feathered  and  the  propulsion  role  would  transition  from  the  ERH  thruster  to  other  airbreathing 
modes.) 

h)  Ground  resonance:  infrasonic  levitation 

Figure  H-26  illustrates  another  potential  refresh-related  mechanism  which  is  worthy  of  note. 
In  operation,  it  can  be  likened  to  the  ‘ground  effect*  mechanism  which  all  winged  vehicles  experience 
at  very  low  altitudes.  The  effect  is  associated  with  a  reduced  flight  power  requirement  and  a  region 
of  pressurized  air  being  momentarily  “trapped"  between  the  wing  and  ground  plane. 

In  this  scheme,  the  ERH  thruster  would  act  as  a  directional  sonic  generator  to  establish  a 
standing  wave  system  between  the  vehicle  lower  surface  and  a  ground-fixed  reflector  (ideally  a  hard, 
smooth  prepared  surface).  As  indicated  in  Fig.  U-25,  nodes  would  be  enforced  at  the  ground  plane 
and  the  vehicle  underside  surface.  Then,  by  simple  phase  modulation,  the  pressure  nodal  planes 
could  be  made  to  shift,  elevating  the  vehicle  into  the  atmosphere  -  to  any  height  still  within  ground- 
effect.  In  flight,  the  engine  would  rebound  small  “elastic”  air  masses  which  move  at  high  velocity 
(up  to  Mach  1)  between  the  vehicle  and  the  ground.  Sub-audible  pulse  repetition  frequencies  could 
be  employed  to  transfer  sufficient  momentum  for  hovering  flight. 

In  essence,  the  vehicle  would  be  levitated  within  the  standing  wavesystem  at  an  upper  node. 
For  example,  a  five  meter  diameter  engine  could  hover  at  5  meters  altitude  or  one-half  wavelength, 
with  an  ultrasonic  frequency  of  33.14  Hz  -  as  shown  in  Fig.  H-26.  (This  calculation  assumes  a  cold 
0°  C  day  when  the  sound  speed  is  331.4  m/s.)  If  sufficient  thrust  exists,  the  5  ra  vehicle  could  also 
hover  at  10  meters  altitude  (one  wavelength)  at  the  same  sub-audible  frequency. 

The  goal  here  is  to  produce  a  time-averaged  overpressure  across  the  engine/vehicle  lower 
surface,  sufflvient  to  support  the  gross  takeoff  weight.  For  the  5  m  diameter,  5656  kg  reference 
platform,  this  equates  to  an  overpressure  of  0.02737  atm.  Whether  it  is  accomplished  by  bouncing 
planar  shock  waves,  subsonic  aerodynamic  “ring  vortices”,  or  any  other  mechanism  involving  an 
“elastic”  air  mass,  the  resultant  effect  is  a  simple  momentum  exchange  with  a  massive  immovable 
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Figure  H-25  Infra-sound  levitation:  ground  effect 
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object:  the  Earth. 

In  order  to  generate  appropriate  infra-round  radiation  patterns,  ERA  thrusters  may  require 
specially-shaped  surfaces,  perhaps  even  flat-bottomed  plugs  (i.e.,  a  ■  0).  Note  also  that  the  re¬ 
ference  ERH  thruster  utilised  LSD-wave  sources  with  idloherts  ‘excitation*  frequencies.  In  order 
to  efficiently  “pump"  the  33  He  standing  wave  system  with  LSD  waves,  unusual  amplitude-sha¬ 
ped  macropulses  (see  Fig.  0-27)  might  be  required.  On  the  other  hand,  the  substantially  lower 
‘excitation*  frequencies  which  are  expected  from  microwave-supported  ‘deflagration*  waves,  may 
couple  well  enough  without  macro-  pulse- amplitude-shaping. 

In  summary,  by  establishing  a  standing  wave  system  between  an  ERH  thruster  plug  and 
the  Earth’s  surface,  a  ‘ground  resonance*  effect  can  be  created.  The  standing  wave  system  may  aid 
in  convecting  ‘spent*  gases  out  of  the  active  ‘detonation*  regions  and  sweeping  fresh  gases  back 
in.  When  operating  in  ground  effect,  ERH  thruster  ‘impulse  coupling*  efficiency  should  increase, 
causing  a  net  reduction  in  propulsive  power  required  for  hover  flight.  If  the  laser  macro-pulse  energy 
is  preserved  (e.g.,  1.6  MJ  for  the  12-  finger  ref.  source),  this  lower  beam  power  level  could  appear 
as  a  diminished  pulse  repetition  frequency  (e.g.,  a  reduction 

C)  Analysis  of  Primary  Receptive  Optics  Area 

This  section  analyzes  the  relationship  of  specific  engine/vehicle  parameters  to  the  primary 
optics  area  available  for  collecting  beamed  power.  The  analysis  applies  only  to  the  reference  vehicle 
configuration  which  receives  the  beam  from  an  orbiting  relay  satellite,  onto  a  primary  optical  surface. 
This  parabolic  reflecting  lens  is  intimately  integrated  with  an  external  compression  inlet  that  closely 
resembles  an  isentropic  spike.  To  avoid  the  requirement  for  a  needle- fine  inlet  spike  lip  (i.e.,  an 
extension  of  the  parabolic  PO  mirror),  an  annular  laser  beam  is  chosen;  hence,  for  simplicity,  a 
non-focusing  conical  tip  is  selected  for  the  reference  design. 

The  general  layout  of  the  engine/optics/vehicle  configuration  is  depicted  in  Fig.  11-28.  The 
present  analysis  applies  only  to  simple  PO/SO  optical  trains;  the  use  of  tertiary  optics  requires 
a  different  set  of  model  input  assumptions.  Note  that  the  angle  a  sets  the  ERH  thruster  apex 
angle,  and  determines  the  maximum  diameter  of  the  engine/vehicle  centerbody;  also,  a  sharp  flow 
turning  angle  exits  inside  the  shroud  -  created  by  the  intersection  of  PO  and  conical  nozzle  afterbody 
surfaces.  This  shrouded  duct  would  function  well  only  with  supersonic  air  flow  velocities;  in  practice, 
the  sharp  edge  would  probably  be  rounded  off  to  some  minimum  radius  -  to  prevent  flow  separation 
at  subsonic  and  low  supersonic  flight  speeds. 

A  major  emphasis  of  this  analysis  was  to  compute  the  amount  of  power  which  may  be 


64 


AMPLITUDE 


Flgur*  11-27  Amplltud«-*hap«d  macro-pul*** 


66 


I 


66 


received  by  the  craft.  This  peak  power  level  ia  directly  proportional  to  the  ana  of  tha  receiving 
optica  and  tha  beam  wavelength.  Note  that  throughout  thia  treatment,  “optica'  refen  to  the  primary 
aperture  area  and  is  independent  of  the  power  beam  wavelength.  The  fraction  of  inlet  spike  frontal 
area  available  for  collecting  and  focusing  the  beam  is  defined  as  the  “active  area  traction”.  This 
active-area  fraction  varies  with  cone  tip  angle  and  the  desired  focus  location  relative  to  the  conical 
bow  shock  wave  (produced  by  the  external  compression  inlet  tip).  The  analysis  assumes  that  the 
conical  bow  shock  angle  is  unaffected  by  the  curvature  of  the  primary  optics  (which  would  act  as  an 
additional  compression  surface). 


1.  Analysis 

The  basic  geometry  of  the  optical  trains  and  the  key  variables  used  in  the  analysis  are  shown 
in  Fig.  11-29.  The  vertical  distance,  Zj,  measured  from  the  PO  optics  base  to  the  conical  shock 
wave  ia  computed  from  inviscid  shock  data.  The  bow  shock  wave  ia  assumed  to  attach  to  the  cowl 
lip  in  order  to  minimize  “spillage*  drag;  thus,  Z*  also  represents  the  cowl  lip  location.  The  vertical 
distance  from  the  PO  optics  base  to  the  parabola  focus  ia  denoted  by  *F”.  The  quantity  D  is  the 
difference  between  Zi  and  F.  In  the  calculations,  D  was  set  equal  to  zero,  requiring  the  parabola 
focus  to  coincide  with  the  cowl  leading  edge.  This  situation  would  not  be  possible  if  the  shroud 
has  a  leading-edge  flap,  because  the  SO  ring  would  have  to  be  moved  aft  within  the  shroud.  The 
equations  presented  below  are  generalized  for  MLL  value  of  D. 

The  active  area  fraction  was  calculated  from  Eqn.  0-1,  with  the  terms  defined  in  Fig.  11-29. 
The  expression  which  describes  the  parabolic  reflecting  surface  is  given  below  In  Eqn.  H-2  [1]. 

Fraction  active  area  *  [T*?  )  |  (//  —  1 ) 

In  order  to  make  the  inlet  spike  a  continuous  surface  (i.e.,  without  ‘drag-increasing’  dis¬ 
continuities)  the  contours  of  the  conical  tip  and  parabolic  reflecting  surface  must  be  matched  by 
equating  the  first  derivatives  of  the  governing  equations(i.e.,  the  slopes).  The  slope  of  the  cone 
surface  is  given  in  Eqn.  II-3,  and  that  of  the  parabola  in  Eqn.  11-4.  The  dimensions  of  the  cone 
base  x*  and  parabola  height  y *  were  then  found  using  Eqns.  II-6  and  II-6,  in  terms  of  the  parabola 
focus  location. 
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Figure  n-2®  Geometry  for  receptive  optic* 


u  .  ,  dy  X-R 
P»nbol»;  jj  ■  ~~TF~ 


(//-  4) 


z‘  =R- 


(II  -  5) 


tan3  a 


(//-«) 


The  parameters  i*  end  I  were  then  found  from  Eqns.  H-7  and  II -8  with  the  quantity  D  given  by 
Eqn.  n-9. 


(//-  7) 


i  ~2/'  R 
tan3  a  tana 


0  *  f  (-i  -  *-)  +  *(_!-  -  *) 

v  tan3a;  vtana  tanp' 


(//-») 


(//  -  9) 


At  this  point,  the  quantity  *D”  (measured  from  the  focus  to  the  cowl)  can  be  specified.  As 
mentioned  above,  D  was  set  equal  to  zero  for  this  analysis;  thus,  Eqn.  II- 9  could  be  solved  for  F 
(the  focus  location),  given  a  cone  half-angle  and  resultant  bow  shock  wave  angle  (determined  by  the 
flight  Mach  number),  and  a  maximum  optics  radius  R.  The  resulting  expression  is  given  below. 


fife)  (//  -  10) 

After  obtaining  F  and  x*  from  Eqn.  II-5,  z,  must  be  found  by  solving  for  the  intersection 
of  the  line  through  points  (R,F)  and  (R-X,0),  and  the  parabola  (Eqn.  II-2).  The  line  location  is 
specified  by  Eqn.  11-11,  with  X  calculated  by  Eqn.  11-12.  The  angle  a  is  that  of  the  afterbody  plug 
nozzle,  as  defined  in  Fig.  11-29.  The  intersection  point  r,  was  then  found  by  equating  Eqns.  0-2 
and  0-11,  resulting  in  Eqn.  11-13  to  be  solved  using  the  quadratic  formula.  The  root  of  interest  is 
negative  due  to  the  choice  of  coordinate  system. 


Line;  y=Lx  +  F(l-j) 


(II-  11) 


X  =  F  tana 


(//  -  12) 
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X?  +  (-2 ft  -  t£)Xi  +  (ft*  -  4/*(l  -  f ))  s  0  (II  -  13) 

In  the  cue  of  a  a  90*,  the  equation  of  the  line  ia  given  by  Eqn.  11-14,  and  the  quadratic 
can  be  solved  for  *  by  using  Eqn.  11-15.  The  negative  root  is  again  the  one  of  interest. 

y  =  F  (II  - 14) 

X?  +  (~2ft)x«  +  (ft*  -  4P)  =  0  (II- 15) 

Using  ti  obtained  from  Eqns.  II- 13  or  IMS,  the  active  area  fraction  can  then  be  calculated 
from  Eqn.  II- 1. 


2.  Results 

The  active  area  fractions  were  calculated  for  cone  semi-vertex  angles  ranging  from  16  to  36 
degrees,  with  afterbody  angles  a  of  0,  30,  45,  60  and  90  degrees,  at  flight  Mach  numbers  of  3,4.6,10 
and  20.  Tabulated  values  for  the  shock  wave  angie,  0,  used  in  the  analysis  were  taken  from  Ref.  [2j 
Plots  of  the  active-area  vs.  cone  angie  and  Mach  number  are  displayed  in  Figs.  0-30  -  0.34. 

The  results  presented  here  are  for  point  designs  having  a  specific  cone  half-angle,  afterbody 
angle  and  flight  Mach  number.  These  results  assume  that  the  focus  of  the  primary  optics  moves  with 
the  cowl  to  catch  the  shock  at  each  condition:  thus,  every  point  represents  a  different  optical  design 
and  inlet  geometry.  The  temptation  to  associate  cowl  translation  to  the  values  of  focal  distance 
presented  in  the  data  tables  must  be  avoided. 

Also,  the  present  analysis  assumed  D  »  0  so  that  the  optica1  ring  focus  coincides  with 
the  cowl  leading  edge;  however,  for  proper  operation,  the  focus  should  actually  be  slightly  outside 
the  shroud  such  that  secondary  optics  could  then  intercept  and  reproject  the  beam  across  the 
notale/afterbody  as  shown  in  Fig.  11-36.  The  desired  D  value  used  in  Eqn.  U-  9  would  have  to  be 
updated  to  include  this  effect. 

From  the  results  in  Figs.  11-30  to  (1-34,  it  is  apparent  that  larger  cone  angles  yield  lower 
active  areas,  as  does  operation  at  high  Mach  numbers.  A  desip  Mach  number  of  3  and  50  %  active 
optics  area  wu  selected  for  the  reference  point  desip. 
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Friction  Active  Arm 


Cone  Half  Angie  (degrees) 


Figure  11-30  Active  area  fraction  v*.  now  cone  angle  (a  =  0*) 
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Figure  H-31  Active  area  fraction  va.  note  cone  angle  (a  =  30°) 
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Friction  Active  Area 
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Figure  11-32  Active  area  fraction  vs.  nose  cone  angle  (a  =  45s) 
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Fraction  Active  Area 
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Figure  n-34  Active  area  fraction  va.  noae  cone  angle  (a  =  90°) 
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CHAPTER m 


LASER  ENERGY  ABSORPTION  INTO  GASES 


A)  Introduction 

Much  of  the  earlier  work  on  laser  impulse  generation  has  been  found  pertinent  to  this  study 
of  the  ERH  thruster  concept.  Yet,  there  are  significant  differences  which  should  be  mentioned.  Many 
of  these  earlier  studies  are  founded  on  slightly  different  assumptions,  and  reach  different  conclusions. 
This  fact  has  been  noted  by  Boni  el  aL  (2j . 

Figure  m  i  shows  a  typical  “perpendicular  beam/surface”  geometry  and  the  resultant  blast 
wave  structure.  Most  studies  also  included  a  direct  Interaction  between  the  laser  and  the  surface 
in  their  analyses,  Le.,  a  direct  ‘burning*  of  the  surface  by  the  laser.  The  nature  of  the  resultant 
absorption  wave  is  then  dependent  on  the  type  of  surface  upon  which  the  beam  is  incident.  The 
surface  Itself  is  participating  In  the  overall  radiative  interaction  with  the  surrounding  gas  and  can 
contribute  some  mass,  via  ablation,  to  the  blast  wave.  The  latter,  however,  is  known  to  have  a 
negligible  effect.* 

Pirri{l5|  Reilly  et  ai  (21j  and  Raiser  [l7^8j,  are  among  those  researchers  who  have  attempted 
analyses  which  consider  the  interaction  between  the  laser  and  the  gas  as  the  principal  reason  for 
the  impulsive  and  thermal  coupling  with  the  surface.  Heating  of  the  surface  is  by  radiation  and 
convection  from  the  gas,  rather  than  by  direct  radiative  heating  of  the  surface  by  the  laser.  The 
latter  does  occur  but  is  considered  negligible  compared  to  the  former.  Note  that  once  the  plasma 
has  been  created,  its  absorption  of  the  laser  radiation  is  so  great  that  the  surface  may  be  considered 
as  being  covered  by  an  optically  thick  gas  and  thus,  does  not  “see”  the  laser  beam.  Radiation  is  also 
considered  to  be  the  dominant  heat  transfer  mechanism  by  which  the  high  temperature  plasma  heats 
the  surface  (i.e.,  rather  than  convection).  The  process  of  the  surface  heating  will  not  be  discussed 
further  here. 

B)  Earlier  Models  of  LSD  Wave  Propagation 

Most  prior  analyses  begin  with  the  gas  dynamic  equations  to  describe  the  motion  of  the  ab- 

*  Most  of  these  earlier  studies  are  in  fact  largely  concerned  with  enhancing  the  impulse  and 
thermal  coupling  of  the  wave  with  the  surface  for  destructive  applications.  Advanced  cooling  meth¬ 
ods,  stronger  materials  and  transversely  incident  pulsed  laser  beam,  can  perhaps  prevent  serious 
damage  to  the  ERH  thruster  surface.  In  fact,  some  research  has  indicated  that  'laser  countermea¬ 
sures’  to  the  destructive  aspects  may  be  easy  to  develop.  This  matter  is  not  discussed  further 
here. 
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sorption  wav*.  Viscous  and  therm «i  low  an  usually  n*gisct*d.  This  is  understandable  considering 
that  tbs  absorption  wave  pressure  may  be  greater  than  these  effects,  even  at  very  high  propagation 
Mach  numbers. 

Raiser  (18),  explains  the  process  of  igniting  an  LSD  wave  (in  air)  at  the  point  focus  of  a 
laser  beam,  as  follows.  Fust,  a  strong  shock  wave  is  generated  in  the  focal  region  where  the  laser 
energy  is  deposited  in  the  gas.  This  shock  travels  through  the  gas,  heating  and  ionising  it,  further 
enhancing  the  absorption  of  more  laser  radiation.  Radiative  energy  is  deposited  in  successive  layers 
of  gas  adjoining  the  shock  front,  and  these  layers  in  turn  act  as  heat  ‘sources’  which  maintain  the 
wave. 

The  wave,  thus,  moves  along  an  'optical  channel’  towards  the  laser  source  and  is  maintained 
by  the  beam.  The  wave  would  continue  to  travel  up  the  beam  as  a  planar  wave  moving  within  a 
cylindrical  regL  i,  at  constant  strength,  if  not  for  its  inevitable  lateral  expansion.  However,  most 
researchers  believe  this  lateral  motion  is  negligibly  smaller  than  the  axial,  and  have  ignored  it  until 
the  end  of  the  laser  pulse. 

Raizer  (17.18),  derived  specific  equations  for  the  velocity  and  pressure  of  the  LSD  wavefront 
from  the  integral  form  of  the  gas-dynamic  equations  using  a  control  volume  drawn  about  the  wave 
(see  Figure  HI-2).  Raiser  obtains  Equations  QM  as  given  below.  In  deriving  these  relations,  he 
apparently  assumes  that  the  wave  moves  steadily,  at  constant  strength  and  without  any  losses, 
immediately  after  its  initiation. 


PqVlsd  -  P* 


Po  +  Po  Visa  -  P  +  P*>3 

e+2«  *_  =  e  +  E  +  ^.  (///- 1) 

Po  2  poVlsd  p  2  1 

where  Vtso  is  the  LSD  wavefront  velocity,  u,  p  and  p,  the  velocity,  pressure  and  density  of  the 

plasma  behind  the  wave,  respectively.  The  subscript  ‘0’  denotes  the  ambient  values  of  these  last 

three  quantities.  *  is  the  laser  beam  intensity.  The  internal  energy  of  the  gas  is  defined  as  e  = 

RT /( 7  -  1)  =  p/[(7  -  l)p),  where  R  is  the  gas  constant  and  T  is  the  temperature. 

From  these  equations,  Raicer  next  obtains  an  expression  for  the  shock  adiabat  as: 


(///- 2) 

This  equation  may  be  combined  with  the  equation  of  state  of  the  gas  (p  =  pRT),  to  obtain  an 
explicit  expression  for  the  shock  adiabat,  p(p),  for  the  supersonic  wave,  as: 
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LASER  BEAM 


Figure  ffl-1  The  LSD  ware  far  tha  caaa  la  point 
(Altar  Plrri,  Root  and  Wn,  Ref.  36) 


Figure  m-2  The  LSD  ware  geometry  considered  by  Raiser  (see  ref.  (18)) 
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(Iff -3) 


Raizer  (18),  graphically  represents  this  shock  adiabat  as  shown  in  Figure  ID-3.  Note  the 
limiting  points  on  this  curve.  Wave  velocity  increases  from  the  initial  state  'O',  along  a  straight  line 
tangent  to  the  shock  adiabat  at  point  ‘J’,  The  Chapman-Jouguet  point  Point  ‘A*  represents  the 
maximum  possible  heating  and  compression  of  the  gas  in  the  detonation  regime,  with  the  assumption 
of  zero  ambient  pressure.  Raizer  claims  that  the  shock-wave  ionized  gas  obtains  additional  energy 
from  the  absorption  of  the  radiation,  and  expands  along  the  line  AJ  until  it  reaches  point  J  at  the 
end  of  the  energy  release.  This  action,  he  believes,  occurs  within  the  wavefront,  which  he  treats  as 
a  discontinuity.  Raizer  admits  that  this  conclusion  requires  verification  by  a  more  detailed  theory 
of  the  internal  wavefront  structure.  The  points  ‘C’  and  ‘D’  represent  the  possible  final  states  of  a 
‘superdetonation’  wave.  Raizer  implies  that,  unlike  chemical  detonation,  this  superdetonation  may 
not  be  sustainable  with  laser  energy. 

Using  Equations  (m-1),  Raizer  shows  that  the  minimum  velocity  of  a  supersonic  wave  in 
the  detonation  regime  is  given  by: 


Fzso  =  (2(T2-lA*  (III -4) 

Some  experimental  results  have  supported  this  finding  (e.g.,  see  Ferriter  et  aL  [16|  ). 

The  pressure  behind  the  wave  front  is  given  by: 


PDW 


=  f*vlse 

7+1 


(///  -  5) 


Raizer’s  analysis  was  not  concerned  with  the  internal  structure  and  temporal  evolution  of  the 
wave.  This  information  is  important  in  obtaining  any  relationship  for  the  impulse  generated.  Thus, 
one  must  utilize  the  inviscid  gas  dynamics  equations  in  their  differential  forms,  in  axi-symmetric 
cylindrical  coordinates: 


dp  ap»  i  dpv  _ 

Hi  +  W  +  7TF-° 


(///- 6) 


<9u  udu  vdu  _  1  dp 

Hi  +  ~HT  ~Hr  pHr 


dv  adv  vdv  1  dp 

Hi  +  ~HT  +  ~5F  ~  —pTr 


(III -7) 

(III  -8) 
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Figure  m-3  Shock  adiabat  of  a  laser  energy  absorption  wave  (tee  ref.  [18] ) 
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where  u  is  ths  axial  velocity,  v,  radial  velocity,  p,  pressure  tad,  p,  the  density.  The  previous  analyses 
did  not  consider  the  energy  equation  in  detail.  Holmes  et  el  (g  ),  used  a  polylropic  state  equation, 
(assuming  a  polytropic  gas),  given  by: 

£  =  Const.  (Ill  -  9) 

where  the  7  used  has  a  value  of  1.2  behind  the  wavefront.  They  assumed  that  the  LSD  wave 
propagates  in  the  following  manner. 

1.  The  wave  starts  out  at  the  strength  and  speed  given  by  Rainer. 

2.  It  propagates  up  the  beam  with  decreasing  strength,  due  to  the  axial  and  radial  rarefaction 
fans  moving  inwards  from  the  outer  shock  boundaries. 

3.  When  the  wave  reaches  a  height  equal  to  its  diameter  (which  is  assumed  to  be  about  equal 
to  the  beam’s  diameter  if  the  wave  speed  is  constant),  the  expansion  of  the  wave,  presum¬ 
ably,  becomes  two-dimensional  or  three-dimensional,  i.e.,  the  wavefront  geometry  is  either 
cylindrical  or  spherical.  This  point  in  time  has  been  denoted  by  fao,  which  is  the  time 
required  for  an  axial  rarefaction  fan  to  reach  the  surface. 

The  rarefaction  fans  travel  at  the  local  speed  of  sound,  c,  ,  which  is  not  constant.  But, 
using  an  average  etl  tip  is  defined  as: 

(///-  10) 

C, 

where  r,  is  the  laser  spot  radius.  Different  definitions  of  tip  have  been  used  in  the  past,  but  they 
are  all  about  equal  to  the  one  which  uses: 


_  Vpw 
2 


(iu-n) 


Holmes  et  al  j  8 ),  assume  that  the  radial  expansion  fans  are  an  effect  of  the  wave’s  decay 
and  hence  consider  only  the  radial  gas  dynamic  equations: 


dp  .dpv  pv_. 

-5t+W  +  T-° 


(III-  12) 


dv 

li 


+  -  i  dp 


(III  - 13) 


along  with  the  polytropic  state  equation.  Raizer  similarly  assumes  that  lateral  losses  are  largely 
responsible  for  the  wave  decay  after  the  pulse  has  terminated. 

Holmes  et  al  [8) ,  then  neglect  the  1/r  term  in  considering  wave  propagation  far  from 
the  centerline. ‘This  assumption  then  renders  the  equation  completely  one-dimensional.  Thus,  some 


*  For  the  spherical  case,  a  2/r  term  appears  instead  of  i/r. 


81 


characteristic  solutions  m ay  bs  obtain*!  for  the  fu  pressure  and  tha  partide,  fan  and  wava  va- 
lodties.  BacauM  the  1/r  term  is  neglected,  Equations  HM2  and  01-13  are  not  valid  at  r  =  0, 
the  centerline. The  solutions  of  these  equations  say  that  straight  rarefaction  fans  propagate  radially 
inwards  at  the  local,  non-constant,  speed  of  sound.  Holmes  d  si  j  g  j,  describe  the  wave  decay 
process  based  on  the  motion  of  each  rarefaction  coming  in,  beginning  when  the  first  one  reaches  the 
centerline. 


Applying  the  boundary  condition  of  tero  velocity  at  the  surface,  the  characteristic  solutions 
yield  a  surface  pressure  of: 

«50  =  [2±1|*p°* 


Table  IIM-Sedov's  Scaling  Laws,  (See  Sedov  (22j) 


UNPOWERED 

POWERED 

PLANAR 

P/Pref 

at 

<C/W2/3 

P/Pref  * 

const 

2/zref 

3 

(t/tr,£)z/3 

z/zref  * 

CYLINDRICAL 

P/PreC 

9 

P/Pref  " 

(t/W1/2 

r/rref 

9 

<■ ''W 1/2 

r/rre£  = 

SPHERICAL 

P^reC 

3 

(t/tref)-S/5 

P/Pref  = 

r/rre£ 

3 

«/W2/5 

r/rref  = 

<t/PraJ>3/S 

Holmes  et  al  [  g|,  then  calculate  the  impulse  for  every  rarefaction,  within  the  time  it  takes 
each  one  to  reach  the  centerline,  using  the  characteristic  pressure  behind  each  rarefaction  -  until  a 
radius  of  three  times  the  spot  radius  is  reached.  Next,  they  apply  the  self-similar  solutions  of  the 
gasdynamic  equations  given  by  Sedov  [22] ,  to  give  the  pressure  decay  to  the  ambient  value.  These 
‘Scaling  Laws’ are  presented  in  Table  IH-1. 
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C)  Ramalt*  at  tlM  Bulkr  Moctals 

The  unpowered  scaling  laws  were  originally  derived  by  Sedov  from  dimensional  analysis 
techniques.  The  powered  scaling  laws  were  supplemented  to  Sedov's  by  assuming  a  linear  energy 
deposition  (i.e.,  substituting  an  energy  varying  linearly  with  time  in  Sedov’s  original  form  of  these 
relations).  These  similarity  relations  are  not  without  their  problems.  One  of  them  is  that  they 
present  a  singularity  at  t  =  0  and  a  zero  value  of  absolute  pressure  u  t  approaches  infinity. 

Pirri  uses  Sedov’s  self-similar  solutions  of  the  gasdynamie  equations  without  considering  the 
motion  of  the  rarefactions.  Indeed  these  equations  appear  to  be  complete  and  valid  solutions  of  the 
much  simplified  gasdynamie  equations  (with  the  1/r  term). 

In  Pirn’s  analysis,  the  wave  is  assumed  to  be  of  constant  strength  while  the  laser  is  on,  then 
decays  according  to  the  planar  scaling  laws  until  the  time  ho  (if  hD  is  greater  than  tp).  The  wave 
then  decays  according  to  the  unpowered  cylindrical  scaling  laws  until  ambient  pressure  is  reached. 
If  t9  >  tjo,  he  applies  the  powered  cylindrical  scaling  laws  before  going  to  the  unpowered. 

Reilly  et  al  (21),  use  what  appears  to  be  a  mixture  of  the  above  analyses.  A  rarefaction 
fan  is  said  to  be  radially  approaching  the  centerline  (from  its  starting  point  at  the  spot  radius),  at 
the  avenge  local  speed  of  sound  c,.  The  surface  pressure  within  this  rarefaction  is  defined  to  be 
the  rarefaction  pressure  from  shock  tube  analysis.  They  assume  that  plsd  is  felt  across  the  surface 
from  the  centerline  to  the  radial  location  of  the  rarefaction  at  time  t.  The  impulse  is  then  calculated 
with  these  pressures  assumed  constant  over  these  areas  until  ho  ■  Next,  they  apply  either  of  Sedov’s 
scaling  laws  depending  on  how  tr  compares  with  the  time  required  for  an  axial  rarefaction  fan  to 
reach  the  surface  (ho  ),  and  the  time  for  a  radial  fan  to  reach  the  centerline. 

Ferriter  et  si  [6j ,  and  Edwards  et  si.  [4] ,  have  used  a  computer  code  called  LASNEX  to 
compute  the  impulse  including  the  interaction  between  the  laser  and  the  surface  in  their  calculation. 
They  claim  to  have  obtained  results  that  agreed  with  each  of  the  above  models.  Holmes  et  sL  (11), 
did  an  experimental  simulation  of  LSD  wave  propagation,  using  explosive  strips,  and  claim  adequate 
agreement  with  their  theory.  They  all  however,  admit  that  more  sophisticated  modeling  needs  to  be 
done. 

Each  of  these  analyses  appears  to  be  based  on  dissimilar  physical  assumptions  but  that  is 
not  entirely  the  case.  All  yield  similar  results  for  small  radii,  for  which  t2D  is  also  small.  Although 
the  characteristic  solutions  are  not  valid  for  small  r,  they  are  used  by  Holmes  et  al  (8j ,  and  Reilly 
et  si  [21]  to  give  a  better  description  of  the  wave  structure  (particularly  in  its  transition  to  a  three- 
dimensional  geometry)  than  Sedov’s  solutions  supposedly  could  in  that  region.  But  if  r  is  small,  the 
additional  assumptions  made  by  Holmes  et  al  (8  )  and  Reilly  et  al  (2I|  will  not  be  advantageous 
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bscaus*  the  rarefaction  fan#  will  not  be  of  weli  enough  defined  structure  to  make  any  difference. 
Hence,  the  characteristic  analysis  is. dot  all  that  advantageous  and  Sedov’s  scaling  laws  can  provide 
adequate  Information  on  the  wave. 

D)  Summary 

This  chapter  has  presented  some  of  the  assumptions  behind  the  earlier  analyses  on  laser 
impulse  generation  for  beams  which  are  perpendicularly  incident  to  a  surface.  The  following  chapter 
will  introduce  another  analytical  model  for  laser  impulse  generation,  one  which  applies  to  LSD  waves 
which  propagate  parallel  (but  not  necessarily  in  contact  with)  the  impulse  surface. 


CHAPTER  IV 


INSTANTANEOUS  'LINE  SOURCE*  ERH  THRUSTER  MODEL  (STATIC) 

A.  Impulse  Generation  Model  for  LSD  Ware*  Propagating  Parallel  to  a  Surface 

When  an  LSD  wave  propagate!  parallel  to  a  surface,  the  type  of  boundary  conditions  im¬ 
posed  on  the  plasma  expansion  differ  from  the  case  of  the  previous  section.  As  the  wave  races  up  the 
laser  beam,  it  is  constrained  from  expanding  in  a  fully  cylindrical  geometry  by  the  close  proximity 
of  an  adjacent  surface  (see  Figure  IV-l).  However,  this  does  not  necessarily  remove  axi-symmetry 
from  the  problem.  The  wave  front  within  the  spot  radius  may  remain  as  in  Figure  HI-1  or  Figure 
HI-2,  and  the  blast  wave  could  still  be  hemi-cylindrical  (or  partly  spherical).  Whatever  the  case,  the 
details  of  the  LSD  propagation  the  blast  wave  expansion  are  more  complicated  than  the  for  previous 
case. 

1.  Initial  ‘Line  Source*  Conditions  Resulting  From  Rapid  LSD  Wave  Heating  of  Air 

That  the  parallel  beam/surface  geometry  is  a  more  complicated  case  is  evident  from  an 
examination  of  Figures  III-l,  III-2  and  IV-l.  Thermal  and  viscous  losses  are  certainly  of  greater 
importance  due  to  the  LSD  wave  propagating  so  close  to  the  surface.  The  wave  front  velocity  can 
then  be  reduced  by  the  presence  of  the  surface.  It  has  been  shown  that  viscous  effects  can  reduce 
the  strength  of  a  shock  wave,  or  make  it  more  oblique.  Reflections  of  rarefaction  fans  shed  off  the 
LSD  wave  front  edge  may  even  be  absorbed  within  the  air/surface  boundary  layer. 

Therefore,  the  terms  representing  these  losses  have  to  be  included  in  Equations  IH-6  -  IH-8. 
However,  the  reasons  for  neglecting  them  before  could  still  apply.  If  the  intensity  of  the  incident 
laser  radiation  is  sufficiently  high,  the  velocity  at  which  the  LSD  wavefront  propagates  across  the 
plate  would  be  hypersonic.  Yet,  it  is  still  possible  for  the  losses  to  be  negligible  compared  to  the 
high  pressures  that  are  generated.  To  simplify  the  following  analysis,  these  same  assumptions  will 
be  made  here. 

Figure  IV-l  shows  a  side  view  of  the  parallely  propagating  LSD  wave’s  structure.  Note  the 
rarefaction  fans  reflecting  back  and  forth  off  the  plate  and  the  shock  boundaries.  As  indicated,  the 
wave  is  quite  a  bit  more  complicated  in  its  axial  propagation  geometry  than  in  the  perpendicularly- 
incident-beam  case.  This  implies  that  the  axial  dependence  of  the  gasdynamic  equations  may  no 
longer  be  as  negligible  as  in  the  analysis  of  Holmes  et  al.  [8] .  However,  if,  again,  the  velocity  of  the 
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LSD  wav*  aero**  the  plat*  is  hypertonic,  then  th*  rarefaction  fan  reflections  may  not  be  “felt’  in  the 
time  required  for  the  LSD  wave  ta  cross  the  surface  (i.e.,  the  laser  pulse  duration). 

In  the  case  of  the  normally-incident  beam,  the  LSD  wave  velocity  was  assumed  to  be  much 
higher  than  the  radial  expansion  velocity.  This  may  also  be  adopted  here  for  the  same  reasons.  The 
LSD  wave  propagation  and  the  blast  wave  expansion  may  then  be  considered  as  axially  independent. 
Hence,  the  blast  wave  expansion  may  be  assumed  to  stem  from  a  thin  ‘line  source*  of  energy. 

3.  Modeling  of  the  Cylindrical  Blast  Wave  Expansion  Process 

The  thin  line  source  model  of  the  cylindrical  blast  wave  expansion,  according  to  the  as¬ 
sumptions  made  above,  practically  neglects  the  axial  propagation  time,  tt,  of  the  wave  across  the 
surface.  Clearly  this  thin  line  has  a  finite  diameter  (equal  to  the  beam  diameter)  and  the  actual 
LSD  wavefront  requires  a  finite  amount  of  time  to  traverse  the  length,  L,  of  the  surface.  Hence,  the 
validity  of  this  model  rests  mainly  on  this  time,  tc,  being  very  small  relative  to  the  expansion  time, 
to. 

The  time  ti  is  also  equal  to  the  pulse  time  of  the  laser,  t,,  since  it  would  make  no  sense  to 
continue  the  LSD  wave  off  the  surface  -  where  no  suitable  impulse  may  be  obtained.  Nevertheless,  the 
plasma  "leftover'  after  the  termination  of  laser  pulse  might  still  continue  off  the  surface  boundaries 
if  it  is  not  suitably  constrained.  At  present,  this  effect  is  not  thought  to  be  a  significant  loss  factor 
in  the  impulse  generation  process;  yet,  it  could  limit  how  small  tr,  (or  L)  could  be.  Bear  in  mind 
that  the  line  source  model  requires  a  ‘long  line'  and  a  short  energy  deposition  time  to  match  the 
implicit  assumptions  of  an  infinitely  long  and  thin  line  with  an  infinitely  small  energy  deposition 
time  of  Sedov’s  equations(22] .  These  requirements  can  be  met  with  a  high  LSD  wave  propagation 
velocity. 

When  the  laser  pulse  is  terminated  (i.e.,  at  t  =  tF),  the  cylindrical  air  plasma  is  not  likely  to 
have  a  radius  much  larger  than  the  spot  radius,  r„  which  has  been  used  in  the  past.  The  LSD  wave 
velocity  across  the  surface  will  in  general  be  very  high  compared  to  the  lateral  expansion  velocity. 
Hence,  the  radius  of  the  cylindrical  blast  wave  at  t  =  tf  may  be  set  equal  to  r,. 

It  is  assumed  that  the  laser  pulse  shape,  (i.e.,  the  laser  beam  intensity  vs.  time,  .'or  t  <tr), 
is  such  that  the  LSD  wave  sees  a  ‘fairly  constant'  laser  flux  after  the  beam  has  passed  through  the 
optical  trair  o  the  thruster  surface.  Also,  for  simplicity,  it  is  assumed  that  the  beam  cross-section 
is  held  som’.  at  semi-circular  immediately  adjacent  to  the  thruster  surface  (see  Figure  IV  2)  This 
ensures  a  plasma  expansion  process  with  a  cylindrical  geometry. 

As  a  result  of  the  line  source  assumptions,  the  problem  becomes  virtually  one-dimensi  >n».  - 
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SHOCK  LSO  RAREFACTION 

FRONT  WAVE  FANS 


Figure  TV-1  Envisioned  structure  of  the  parallel?  propagating  LSD  war* 


Figure  TV- 2  Perspective  view  of  the  blast  wave  structure  for  the  tine  source 

model 


8' 


ud'iymoitriciL  Neglecting  Iomm,  the  wave  than  could  still  travel  along  the  surface  with  a  velocity, 
Vise,  at  given  by  Raizer. 


VUD-W-  Ojl*  (JV-1) 

Then  the  time  tf  =  (l,  can  be  defined  ac 


••'vhs  ('V‘J) 

This  definition  satisfies  the  requirements  of  the  fine  source  model  as  previously  mentioned,  if  the 
laser  intensity,  4,  is  large. 

The  equations  are  basically  the  same  as  the  equations  used  in  the  previous  analyses  (p  bring 
radially  dependent  only).  Equations  m-1  can  be  applied,  without  modification,  to  the  line  source 
helping  to  link  the  model  to  the  hut  that  the  blast  wave  is  bom  of  laser  energy  deposition.  Tins 
further  justifies  the  applicability  of  EquationQH  to  the  line  source  modeL  The  boundary  conditions 
of  sero  surface  velocity,  used  to  obtain  the  surface  pressure  from  the  gasdynamic  equations,  are  also 
the  same  since  the  blast  wave  is  also  radially  dependent  only.*Slnce  t,  is  small,  it  is  assumed  that 
the  surface  virtually  feels’  this  same  pressure,  plsd,  at  the  end  of  the  laser  pulse. 

PLSD  —  ^PDW  {IV  -  3) 

Note  that  plsd  is  the  pressure  with  which  the  unpowered  cylindrical  decay  is  assumed 
to  begin.  Figure  IV-3  displays  plsd  as  a  function  of  laser  intensity,  4,  for  several  altitudes  of 
interest.The  4  dependence  ht.ps  to  incorporate  the  fact  that  the  blast  wave  is  based  on  laser  energy 
deposition  into  the  fine  source  model.  However,  the  incremental  impulse  generated  during  the  transit 
of  the  LSD  wave  (t  <  tr),  is  assumed  to  be  negligible  even  though  the  initial  pressure  may  be  quite 
large.  This  assumption  is  made  because  the  LSD  wave  is  propagating  parallel  to  the  surface  and 
hence,  perpendicular  to  the  thrust  direction. 

At  time,  tao,  a  rarefaction  fan  propagates  radially  to  the  centerline  of  the  blast  wave  while  it 
is  at  the  radius,  r„  assumed.  Then,  the  blast  wave  expansion  process  assumes  a  cylindrical  geometry 
after  the  time  tap  (this  time  might  be  better  defined  as  t ,  but  is  kept  as  such  because  it  turns  out 
to  be  the  same  as  that  defined  in  the  previous  analyses). 

=  ?  {IV  -  4) 

e» 

*  Because  this  is  not  an  infinite  line  source,  the  non-ideal  3-D  effects  that  perturb  the  axi- 
symmetry  of  the  problem  may  actually  alter  the  surface  pressure. 


Figure  IV-3  Peak  aurface  (reference)  pressure,  Plsdi  ▼*•  Incident  laser 


whan  e,  again  has  an  avorage  value  oh 


(IV  -  5) 


Hence,  the  wave  has  a  itructure  and  temporal  history  that  is  largely  seif-similar  (see  Figures  VL4). 


3.  Computation  of  Impulse  Coupling  to  the  Surface 


The  impulse  on  the  surface  Is  defined  as: 

r<»  rf<() 

/  =  J  l  [p(t)-poliA*  (/V-6) 

where  dA  ■  L  dr.  The  pressure  history  used  in  this  model  from  Table  1,  is  displayed  in  Figures 
IV-4.  By  the  previously  mentioned  assumptions  and  this  equation,  from  t  »  0  to  f  *  ho,  *the 
surface  feels  the  impulse  of  a  ‘planar  powered  wave'  defined  by: 

Iff  ■  to'LpLSDhD 

During  this  time,  the  decremental  ambient  contribution  is: 

Iff  —  -Wpoho 

From  tjo  to  to,  the  unpowered  cylindrical  scaling  law  gives: 

-11  (IV  -  9) 

and  for  the  ambient  term: 

-  1)  (IV  - 10) 

The  sum  of  these  terms  defines  the  total  impulse,  1. 

From  the  unpowered  cylindrical  scaling  law,  the  decay  time  may  be  calculated  as: 

t0  *  *2£££2£  (/V  - 11) 

Po 

It  may  be  assumed  that  the  blast  cone  can  be  refreshed  with  unheated  air,  moving  radially 
inwards  with  the  last  rarefaction  fan  as  the  wave  collapses,  as  to  is  approached.  Then,  the  plasma 

*  tp  is  not  used  as  the  reference  time  here  because  it  wu  assumed  that  there  was  no  contri¬ 
bution  to  the  impulse  beforehand. 


ir.LpisohDlJj^ 


(IV -7) 

(IV-  8) 
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decay  and  the  impulse  generation  cycle  end  at  (a  to.  With  tre/ree*  essentially  aero,  the  maximum 
pulse  repetition  frequency,  PRF,  of  the  laser  is  then  defined  as: 


PRF-I 

*o 

Hie  peak  input  power  during  a  laser  pulse  is  given  by: 


(IV  - 12) 


'-S' 5* 

and  the  laser  pulse  energy  is: 

Er  =  Pip 


(IV  - 13) 
(TV -14) 


Note  that  this  definition  of  the  laser  pulse  energy  also  helps  to  incorporate  the  fact  that  the  blast 
wave  is  obtained  from  a  linear,  laser  energy  deposition  over  a  time,  t,.  The  duty  factor  of  the  cycle 

is 


DF  a  t,PRF  (IV  - 15) 

From  this,  the  tune-averaged  laser  power  deposited  into  the  gas  over  the  impulse  generation  cycles 
may  be  calculated  «P»Px  DF.  The  time-averaged  thrust  over  the  cycle  is  defined  as: 


T-IxPRF  (IV  - 16) 

A  measure  of  the  impulse  delivery  efficiency  is  given  by 


oc  3 


I  T 


(IV  - 17) 


winch  is  the  definition  of  the  coupling  coefficient,  CC.  From  the  unpowered  cylindrical  scaling  law, 
the  final  radius  of  the  wave  may  be  calculated  as: 


ro 


(IV  - 18) 


The  final  width  of  the  cylindrical  blast  wave  is  twice  this  'alue,  and  should  not  approach 
the  length  L,  or  the  assumption  of  axi-symmetiy  is  violated.  Also,  the  plasma  must  not  expand 
beyond  the  finite  width  of  the  plate  during  t .  Otherwise,  a  spherical  development  of  the  blast  wave 
must  be  invoked. 

A  computer  program  is  used  to  calculate  the  above  variables  for  an  interesting  range  of 
beam  intensities  and  atmospheric  conditions.  The  laser  radius  is  set  at  1  cm  and  the  surface  length, 
25  m.  A  value  of  1.2  is  used  for  as  in  earlier  analyses.  The  results  are  discussed  in  the  next  chapter. 
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B)  Remits  of  tbo  lias  soares  modal 

The  main  concern  of  this  analyais  ia  the  efficiency  of  impulse  delivery  to  the  thruster  surface. 
Critical  parameters  were  varied  to  observe  and  understand  their  effect.  The  remits  of  the  line  source 
calculations  are  discussed  below.  The  single-pulse  Impulse  is  plotted  in  Figure  IV-6.  Note  that  more 
impulse  is  imparted  to  the  surface  by  increasing  the  laser  intensity.  This  is  of  course  an  advantage 
which  can  be  exploited  by  the  ERA  thruster  concept.  Also  note  that  the  impulse  decreases  with 
increasing  altitude.  This  is  generally  true  for  all  airbreathing  propulsion  systems.  This  decrease 
occurs  because  lens  air  is  available  for  static  thrust  generation  as  the  air  density  and  pressure 
decrease  with  altitude.  The  time-averaged  thrust,  displayed  in  Figure  IV-6,  also  has  these  same 
trends. 

The  peak  laser  pulse  energy  and  power  absorbed  by  the  LSD  wave  are  shown  in  Figures  IV-7 
and  IV-8.  They  more  or  less  agree  with  Reilly's  observation  that  a  reduction  in  the  laser  intensity  is 
required  to  generate  LSD  waves  at  higher  altitudes  (20j.  Figure  IV-9  shows  the  average  power  going 
into  a  cycle  as  a  function  of  incident  laser  intensity.  This  figure  helps  to  demonstrate  how  the  input 
laser  power  is  distributed  throughout  the  eyrie. 

However,  the  best  measure  of  the  impulse  delivery's  efficiency  is  the  coupling  coefficient.  This 
has  been  plotted  in  Figure  IV-10  against  the  laser  intensity.  As  indicated,  coupling  coefficients  are 
quite  promising,  bring  on  the  order  of  1000  N/MW.  Note  also  that  the  thruster  efficiency  increases 
with  altitude  and  laser  intensity.  The  increases  imply  a  more  efficient  absorption  and  application  of 
laser  energy  by  the  air  at  higher  altitudes. 

Apparently,  the  coupling  coefficient  levels  off  at  the  higher  intensities  and  altitudes.  This 
may  occur  because  the  air  is  approaching  a  limit  in  its  ability  to  absorb  more  laser  energy;  a  condition 
which  can  be  linked  to  use  of  Equation  m-9  which  originates  from  the  entropy  equation  (i.e.,  the 
second  law  of  thermodynamics,  see  Holmes  et  si.  [8] ).  As  derived,  the  coupling  coefficient  really 
depends  only  on  the  gas  properties  and  the  laser  beam  intensity.  This  point  is  also  mirrored  by 
the  LSD  wave  pressure  (shown  in  Figure  IV -3)  which  is  a  stronger  function  of  laser  intensity  than 
altitude.  The  variations  with  altitude  and  intensity  is  also  displayed  in  the  maximum  blast  wave 
expansion  radius,  rc  in  Figure  IV-11.  Note  that  ro  is  also  more  sensitive  to  increases  in  beam 
intensity  than  those  in  altitude.  The  maximum  plasma  radius  is  directly  proportional  to  the  square 
root  of  the  LSD  wave  pressure  and  inversely  proportional  to  the  ambient  pressure. 

The  parameter  ro  also  defines  the  radius  over  which  the  impulse  and  the  mean  thrust  are 
delivered.  However,  these  two  quantities  do  not  appear  to  be  as  sensitive  to  the  area  increase.  This 
may  be  because  the  other  variables  they  depend  on  (like  p  and  t),  are  not  very  sensitive  to  the 
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Pulse  Energy,  (kJ) 


Flux,  (W/sqcm) 


generation 


Pignre  IV-7  Peak  laaer  power  supplied  to  the  L8D  wave  for  the  plasma 


Figure  IV-8  Laser  pulse  energy  supplied  to  generate  the  line  source 


altitude  effects  end  thus,  counterbalance  the  area  increase.  Another  panmeter  on  which  the  mean 
thrust  and  the  impulse  depend  is  the  pulse  repetition  frequency,  PRF,  which  is  shown  in  Figure 
IV- 12,  to  decrease  with  Increasing  beam  Intensity,  as  well  as  altitude.  Evidently,  the  leas  dense  air 
at  higher  altitude  takes  more  time  to  decay  from  plsd  to  po,  and  clearly,  has  the  immediate  effect 
of  decreasing  PRF,  since  PRF  =  1/to. 

Largely  because  of  the  necessity  for  low  duty  factors  (he.,  DF  =  t,/<o),  the  impulse  model 
is  valid  for  intensities  greater  than  about  5  x  WWfem2.  Note  in  Figure  IV.10  that  the  LSC/LSD 
transition  region  for  10.6  micron  laser  radiation  lies  between  10*  -  10 7WfctrP.  At  the  shortest 
wavelengths  of  interest  (e.g.,  0.36  micron)  this  transition  will  scale  to  AMxWW fcm2  at  the  upper 
limit.  (Recall  that  LSC  waves  have  subsonic  velocities,  and  that  the  model  treated  here,  applies 
only  to  LSD  waves.)  Low  duty  factors  (i.e,  about  0.1)  require  high  LSD  wave  velocities,  which  goes 
back  to  the  need  for  elevated  laser  intensities  to  secure  short  laser  pulse  times  (see  Figure  IV.13). 

At  the  highest  feasible  intensities,  for  10.6  micron  radiation,  performance  is  limited  by  the 
breakdown  condition  to  1  -  2  x  10 9W/cma.  However,  the  breakdown  threshold  can  be  extended  to 
higher  intensities,  since  this  limits  scales  as  wavelength  to  the  3/2  power  (see  Reilly ).  For  example, 
this  could  permit  operation  at  flux  levels  as  high  as  WlW /cm2  with  0.5  micron  radiation;  however, 
severe  demands  would  be  placed  on  the  secondary  optics. 

It  is  obvious  from  the  above  analytical  results,  that  the  static  performance  predicted  for  the 
‘line-source’  ERH  thruster  Is  exceptionally  promising.  An  exemplary  application  for  this  advanced 
propulsion  concept  is  discussed  in  the  following  chapter. 

C)  Summary 

In  retrospect,  it  must  be  stated  that  the  assumption  of  a  line  source  wu  somewhat  restrictive 
to  the  analysis.  As  shown  above  in  section  B,  the  validity  of  the  model  is  confined  to  &  rather  narrow 
range  of  intensities  ranging  from  0.5  -  1  x  VfiWfcm2  at  the  low  end,  to  1  -  2  x  lQ9Wfcm2  at  the 
breakdown  limit,  for  10.6  micron  radiation.  However,  the  breakdown  threshold  scales  as  wavelength 
to  the  3/2  power,  so  this  upper  limit  could  easily  be  shifted  upward  by  one  or  two  orders  of  magnitude 
-  amply  by  invoting  higher  laser  frequencies  (e.g.,  0.35  to  0.5  microns).  Thus,  a  larger  range  of 
validity  for  the  model  and  for  operation  of  the  ERH  thruster  is  possible. 

The  final  set  of  the  gasdynamie  equations  imply  that  the  LSD  wave  propagation  velocity  is 
independent  of  wavelength,  but  this  is,  in  fact,  misleading.  The  exact  physics  of  the  laser  energy 
deposition  is  not  fully  taken  into  account  in  the  model.  A  more  rigorous  solution  would  incorporate 
the  laser  radiation  as  a  source  term  in  the  energy  equations,  and  include  re-radiation  heat  transfer 


Max  Blast  Wave  Radius,  (m) 


Figure  IV- 11  Maximum  cylindrical  blast  wars  radios,  no 


Figure  IV. 13  Pulae  Repetition  Frequency,  PRF,  re.  flux  level 
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to  immediate  surroundings.  Jackaon  and  Nielsen  (©J ,  and  Su  and  Boni  (26)  investigated  tho 
effect*  of  radiation  for  LSC  waves,  but  a  similar  analysis  for  LSD  waves  has  yet  to  be  performed. 

Such  effects  were  also  ignored  in  prior  laser  impulse  studies  (by  others),  with  LSD  waves. 
These  researchers  used  the  blast  wave  scaling  laws  and  a  polytropic  equation  of  state  to  include  the 
radiation  effects. 

Despite  the  many  assumptions,  earlier  analyses  have  found  good  agreement  with  laser  im¬ 
pulse  experiments.  These  experiments  were  generally  performed  with  10.6  micron  laser  radiation 
(see  Holmes  d  tL  [8]  and  Ferriter  d  ti  (5j ).  More  experimental  work  is  undoubtedly  needed  to 
verify  performance  projections  for  the  line  source  models. 

Because  of  the  model’s  simplicity,  further  improvements  in  the  theoretical  model  is  also  in 
order  -  to  include  a  more  complete  representation  of  the  physics.  A  more  elaborate  model  winch 
could  include  the  breakdown  process  and  the  actual  dynamics  of  the  LSD  wave  propagation  along 
the  plate  would  be  significantly  better.  The  following  chapter  includes  the  latter  in  an  alternate 
approach. 


CHAPTER  V 


PROPAGATING  LSD- WAVE  ERH  THRUSTER  MODEL 

In  the  previous  chapter,  laser  energy  deposition  Into  the  working  fluid  was  assumed  to  bo 
instantaneous,  resulting  in  an  infinite-length  cylindrical  blast  wave.  Actually,  energy  absorption 
follows  the  moving  shock  front  and  creates  a  propagating  near-cylindrical  blast  wave.  Here  in 
Chapter  V,  ERH  thruster  models,  which  attempt  to  incorporate  this  dynamic  behavior  are  detailed. 
Two  models  are  developed,  one  for  a  hovering  fight  platform  (subsonic  velocities),  the  other  for 
supersonic  flight  velocities.  These  are  referred  to  as  the  "static*  and  "dynamic”  ERH  thruster 
models,  respectively. 

Both  models  are  applied  to  a  simple  engine  configuration,  consisting  only  of  a  lower  impulse 
surface,  that  is  integrated  into  the  lower  fight  platform  surface.  As  illustrated  in  Fig.  V-l,  this 
thruster  surface  is  assumed  to  be  a  circular  frustum  which  comprises  the  entire  aft  section  of  the 
craft.  Abo  note  that  an  annular  cowl  wraps  about  the  vehicle  midaection,  and  that  the  forebody 
has  a  conical  shape  which  acts  like  an  external  diffuser  (at  supersonic  velocities). 

To  generate  thrust,  a  series  of  laser-induced  cylindrical  blast  waves  are  initiated  adjacent 
to  the  thruster  surface;  as  these  blast  waves  expand,  impulse  b  delivered  to  thb  surface.  Before 
modeling  the  ERH  thruster  performance,  it  b  lint  instructive  to  examine  the  laser-  induced  blast 
wave  phenomenon  which  b  fundamental  to  actual  thruster  operation. 

A.  BASIC  OPERATING  PRINCIPLES 

When  a  high  intensity  laser  beam  b  focused  In  air,  a  detonation  waveb  formed;  the  structure 
of  thb  wave  consists  of  a  shock  followed  by  an  absorption  zone.  Within  the  absorption  cone,  laser 
energy  is  transferred  Into  the  fluid  in  the  form  of  thermal  energy,  by  the  mechanism  of  inverse 
bremsstrahlung.  Thb  Laser-Supported  Detonation  (LSD)  wave  propagates  up  the  laser  beam  at 
hypersonic  velocities  and  leaves  behind  a  high  temperature,  high  pressure  plasma  cone.  This,  in  tun, 
expands  into  the  ambient  air  much  like  a  blast  wave  generated  by  a  conventional  chemical  explorion. 
Because  LSD  waves  propagate  toward  the  laser  source  much  faster  than  the  radial  expansion  of  the 
hot  plasma  behind  the  detonation  front,  the  resulting  blast  waveb  cylindrical  in  shape  (i.e.,  for  long 
laser  pubes). 

The  ERH  thruster  utilizing  LSD  blast  waves  has  the  following  operation.  First,  several  laser 
beams  are  projected  parallel  to  the  thruster  surface,  and  focused  to  initiate  LSD  waves.  As  thsee 
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LSD  w*vw  propagate  xna  ti is  thruster  plate,  Ugh  prwan  cylindrical  blast  w avss  w»  formed 
and  subsequently  expand.  The  portion  of  the  Mart  ware  in  contact  with  the  thruster  surface  exerts 
a  pressure  force  upon  the  surface.  This  Impulsive  loading  results  in  thrust  Once  the  Blast  wave 
pressure  decays  to  the  local  ambient  level,  the  thruster  surface  must  be  cleared  of  the  spent,  hot 
plasma  and  replaced  with  ambient  air  before  another  set  of  LSD  waves  can  be  initiated.  This 
exchange  of  the  hot  plasma  with  cool,  unprocessed  air  Is  equivalent  to  the  heat  rejection  portion  of 
a  conventional  thermodynamic  cycle.  In  the  present  analysis,  this  exchange  process  is  referred  to  as 
“surface  refresh*. 

From  tite  above  <fiseussion,  it  is  evident  that  the  ERH  thruster  is  a  repetitively-pulsed 
engine,  in  which  its  Pulse  Repetition  Frequency  (PRF)  is  governed  essentially  by  the  Mast  wave 
decay  time  and  the  thruster  surface  refresh  time;  the  total  cycle  time  is  dependent  on  the  local 
ambient  air  pressure  and  density.  This  Implies  that  thruster  performance  varies  with  altitude  and 
flight  velocity.  Before  the  ERH  thruster  performance  is  analysed,  the  structure  of  the  LSD  wave  is 
investigated  and  a  static  model  for  the  ERH  thruster  is  developed. 

B.  LSD  BLAST  WAVE  ANALYSIS 

A  laser  supported  detonation  (LSD)  wavecan  be  analysed  as  if  it  were  a  chemical  detonation; 
the  chemical  heat  of  reaction  is  replaced  by  the  absorbed  laser  heat  flux.  Therefore,  the  ZND 
detonation  wave  model,  developed  by  Zel'dovich,  Von  Newman,  and  Doring  for  chemical  detonations, 
can  be  adapted  and  used  to  analyse  the  LSD  wave  structure.  In  the  ZND  model,  chemical  detonation 
wave  structure  is  modeled  as  a  one- dimensional  shock  front  followed  by  a  high  speed  deflagation;  for 
LSD  waves  the  deflagration  is  replaced  by  a  thin  absorption  cone.  In  LSD  waves  there  is  a  coupling 
between  the  fluid  dynamics  and  the  absorption  kinetics.  The  leading  shock  ionises  the  air  thereby 
allowing  the  formation  of  the  absorption  sons,  and  likewise,  the  laser  energy  absorbed  within  this 
tone  drives  the  shock.  Because  of  this  coupling,  the  detonation  wave  structure  occurs  only  when 
laser  intensities  are  above  10*  -  10 rW/CM7  (for  10.6  pm  radiation).  Raiser(17],  using  the  ZND 
wave  structure,  was  able  to  derive  the  following  wave  velocity  expression  for  a  LSD  wave  propagating 
into  quiescent  gas. 

VD  =  Mt2  -  l)£|l/8  (V-l) 

PD 

Tide  velocity,  is  the  rate  at  which  the  detonation  front  travels  with  respect  to  a  stationary  observer; 
or  from  laboratory  coordinates. 

During  the  derivation  of  Eqn.  V-l,  the  detonation  wave  has  been  viewed  from  an  Eulerian 
frame  of  reference.  It  is,  however,  generally  easier  to  examine  the  detonation  wave  etructure  from  a 
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Lagrangian  frams.  In  the  Lagrangian  frun*,  the  oh— wr  io  unmod  to  bo  traveling  on  tho  shock 
front,  typically  referred  to  u  shock  coordinates  To  convert  from  Eulerian  coordinates  to  shock 
coordinates  the  following  transformation  is  applied: 

n-nf-Vo  (V  -2) 

In  Fig.  V-2  the  structure  of  the  LSD  wave  is  shown  in  shock  coordinates.  Observe  that  a  radial 
rarefaction  fan  begins  to  propagate  inward  toward  the  blast  wave  crater  directly  behind  the  absorp¬ 
tion  sons  Also  note,  that  regions  of  steady  pressure  and  expansion  velocity  exist  ahead  and  behind 
the  rarefaction  fan.  In  an  actual  blast  wave,  the  converging  expansion  tan  is  reflected  at  wave's 
centerflne  and  forms  a  second  rarefaction  wave  which  propagates  outward  from  the  wave’s  crater. 

In  this  analysis,  only  the  first  rarefaction  wave  will  be  of  interest.  Once  the  first  fan  reaches 
the  centerflne,  at  location  X*  in  Figure  V-2,  the  blast  wavs  decay  Is  approximated  using  seif-similar 
theory.  Previous  Investigators  [8,15]  have  successfully  appfied  seif-similar  blast  wave  theory  away 
from  the  LSD  detonation  front  Here,  the  blast  wave  will  be  analysed  in  two  sections:  first,  using 
the  method  of  characteristics  in  the  region  between  the  absorption  front  and  the  plane  where  the 
first  fen  reaches  the  crater  fine;  and  thereafter,  using  self-similar  blast  wave  theory. 

In  this  model,  axial  relaxation  of  the  LSD  blast  wave  is  neglected,  and  radial  expansion  is 
assumed  to  dominate  the  preeeure  decay  process.  Abo,  the  axial  velocity  b  amumed  email  when 
viewed  from  the  Eulerian  reference  frame.  These  assumptions  permit  the  use  of  the  following  one- 
(fimenaionai  unsteady  flow  equations  to  describe  the  blast  wave  expansion: 


(V-3a) 

(V  -34) 

(V  -  3c) 

For  the  present  development,  the  last  term  of  the  continuity  equation  b  dropped.  Using  an  order  of 
magnitude  analysis,  it  can  be  shown,  that  the  error  introduced  by  neglecting  thb  terms  b  small;  its 
removal  simplifies  the  development  of  the  characteristic  equations.  Using  the  method  of  characteris¬ 
tics,  the  shove  system  of  partial  differential  equations  are  recast  into  ordinary  differential  equations, 
such  that  the  fluid  properties  are  described  along  some  characteristic  path.  Hence,  the  continuity 
and  momentum  equations  can  be  combined  and  manipulated  to  give  the  following  reeult: 


<7-1 

»  J*  m  Const 

(V  -4a) 
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a  - - r 

7-1 

*  /“  »  Const. 

(V  -  44) 
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These  tqntitoi  in  the  well  known  Riemaim  invariants  which  in  valid  alone  the  following  paths 
in  the  length-time  plana. 

^  «  +  C ;  C+  Characteristic  (K-5a) 

^  -  a  -  C ;  C~  Characteristic  (V-5*) 

The  above  four  expressions  a re  fundamental  to  the  analysis  of  LSD  blast  waves  near  the  detonation 
front. 

As  previously  mentioned,  the  pressure  decay  process  for  the  rest  of  the  blast  wave  is  ap¬ 
proximated  using  self-similar  cylindrical  blast  wave  theory.  Expressions  predicting  this  type  of  wave 
decay  are  easily  derived  from  dimensional  analysis. 

Assume  that  the  column  of  high  pressure  plasma  generated  by  a  LSD  wave  can  be  approx¬ 
imated  by  a  cylindrical  ’  set  wave  resulting  from  an  intense  line  source  explosion.  See  Fig.  V-3. 
The  dependent  variables  of  this  blast  wave  expansion  are  velocity  (U),  pressure  (P),  and  density  (p), 
while  the  independent  variables  are  blast  wave  radus  (E),  and  time  (t).  The  physical  constants  of 
this  decay  problem  are  the  energy  released  during  the  explosion,  as  well  as,  the  initial  density  and 
pressure.  The  dimensions  of  these  variables  are: 


(K-6.) 

<•}-* 

(V-») 

w-m 

(V-6c) 

w** 

(V-64) 

It  can  be  shown,  that  the  initial  pressure  is  not  important  in  strong  shock  waves;  thus,  the  only 
parameter  constants  are  the  energy  released  and  the  initial  air  density.  The  objective  of  dimensional 
analysis  is  to  form  nondimensional  parameter  groups.  Observe  that  the  two  independent  variables 
can  be  combined  to  form  the  following  similarity  parameter 


Const. 


(V-7) 


Hence,  r  is  proportional  to  t-1/3.  Using  the  strong  shock  assumption,  following  expresrion  can  be 
obtained: 

/’-(—y)**8  (V-8) 

Differentiating  Eqn.  V-7  and  substituting  it  into  the  above  expression  results  In  an  expression  for 
pressure: 

P  m  (*'-•) 
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In  the  above  equation  the  bracketed  Ura»  in  constant  and  H  it  possible  to  form  tho  followinf  ratio: 


(V  - 10) 


An  additional  ratio,  can  be  derived,  uainf  Eqn.  V-7: 


r 

*W 


(V  - 11) 


Hence,  if  blast  wave  condition*  are  known  for  a  given  time,  then  the  above  expressions  determine 
■ubeequent  blast  wave  pressor*  decay. 

To  determine  the  fluid  conditions  directly  behind  the  absorption  tone,  in  the  Chapman- 
Jouguet  plane,  the  following  equations  developed  for  the  ZND  model  can  be  used: 


I  (•'-») 

Co,-{  i+T)*1  (V-  124) 

Ccj  *  1*Ci  (V  ~  12 c) 


Recall  that  the  Chapman-Jouguet  plane  is  located  Just  behind  the  reaction  tone;  in  the  case  of  LSD 
waves  it  is  behind  the  absorption  tone.  Chapman-Jouguet  theory  predicts  that  stable  detonation 
waves  eodst  only  when  the  relative  velocity  of  the  flow  in  the  Chapman-Jouguet  plane,  with  respect 
to  the  wave  front,  is  equal  to  the  local  sound  speed;  therefore: 

Vo-«cj*Cc/-|2£±1*  (V  - 13) 

'->Ci 

Combining  the  definition  of  the  ideal  gas  sound  speed,  the  expression  for  polytropic  expansion  and 
Eqn*.  V-12c  and  V-13  the  pressure  behind  the  detonation  is  given  by: 

Assuming  that  the  axial  velocity  is  sero,  the  Initial  Impulse  pressure  and  sound  speed  are  given  by 
the  following  expressions,  respectively: 

(V  - 15) 
[V  - 16) 

Using  the  above  results,  the  Riemann  invariants  and  the  self-  similar  expansion  expressions,  it  is 
now  possible  to  determine  the  expansion  of  the  blast  wave  and  subsequent  thrust  generation.  The 
pressure  decay  near  the  detonation  front  is  now  examined. 
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From  Fig.  V-4  it  is  bmh  that  ths  first  rarefaction  fen  is  a  simpls  backward  feeing  wavs; 
thus,  flow  conditions  behind  it  can  be  calculated  using  the  invariant  along  the  C+  characteristic. 
Because  the  radial  velocity  ahead  of  the  expansion  fen  is  equal  to  aero,  the  invariant  is  written  as: 

=  ** +  ^  “ 17) 
Combining  this  with  the  perfect  gas  relations,  the  above  equatton  can  be  written  as: 


w  >  ^(C,  -  <?,)  (V  - 18.) 

«-:rTT<7r)‘<>-<&)*1)  (r-u.) 

This  expression  can  then  be  solved  numerically  with  the  following  well  established  shock  tube  result, 
to  determine  the  pressure  behind  the  rarefaction  fen  (Pri). 


(V  - 19) 


Substituting  this  result  into  equation  18c,  the  particle  velocity  behind  the  refraction  fen  can  be 
determined.  Knowing  the  particle  velocity,  it  is  now  possible  to  find  the  location  of  the  radial  shock 
front.  At  X1,  the  radius  of  the  blast  wave  is  defined  as  &.«/■  This  reference  radius  is  used  to  start 
the  seif-similar  analysis. 


C.  STATIC  THRUST  MODEL 


The  impulse  loading  on  the  thruster  surface  due  to  an  expanding  cylindrical  plasma  blast 
wave,  is  given  by  the  integral  relation: 


I 


// 


PiAdt 


[V  -20) 


Using  the  above  expression,  the  impulse  delivered  to  the  thruster  piate  can  now  be  found. 

To  determine  the  impulse  contribution  of  Region  1  (see  Fig  V-2)  note  that  its  geometry 
remains  fixed  as  the  LSD  wave  propagates  over  the  thruster  surface.  The  location  of  X '  can  be 
determined  by  integrating  Eqn.  V-5b,  solving  for  the  integration  constant  and  multiplying  the 
result  by  Vo. 


X'*2  R, 
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{V  -21) 


SHOCK  ISO  RAREFACTION 

FRONT  WAVE  FANS 


Figure  V-4  LSD/Blast  Ware  Expanatoa  (from  Ref.  8) 


ThwWore  the  impulse  contribution  of  this  region  b  simply: 

/i-2r.(r.Pi+rf.,ft)t  (V  -  22) 

where  the  t  b  equal  to  the  time  necessary  for  the  detonation  front  to  propagate  aeroea  the 
thruater  surface.  Since  the  laaer  pube  b  terminated  when  the  detonation  front  reachea  the  end  of 
the  thruater  plate,  t  b  equal  to  the  laaer  pube  duration  time  tp. 

The  impube  contribution  for  the  region  X  >  X'  in  Fig.  V-2,  can  be  determined  from  Eqs. 
10, 11  and  20,  which  are  combined  to  give: 

[**  (V  -23) 

Jt„f  *r*f  *rtf 

where  to  b  the  time  neceaaary  for  the  blaet  wave  to  decay  to  the  local  ambient  premure,  and  Utj  b 
the  time  required  for  the  first  rarefaction  fan  to  reach  the  blast  wave  center.  To  account  for  varying 
flow  properties  along  the  plug  noule,  the  impube  In  the  radial  expansions  regime  b  numerically 
integrated  over  the  thruster  surface,  using  a  given  step  sbe  of  delta  x.  Thus,  the  total  impulse 
generated  by  the  LSD  waveb  the  sum  of  the  contributions  from  the  two  expansion  regimes,  aa  given 
below: 

/total  «/i+jP/<  “  24) 

Note  that  the  impube  force  acts  along  the  thruster’s  surface  normaL  For  the  ERH  thruster  engine 
configuration  examined,  all  of  the  impulse  loading  predicted  by  the  above  expression  does  not  con¬ 
tribute  to  engine  thrust;  only  the  impube  vector  component  aligned  with  the  vehicle  flight  direction. 
Thus,  magnitude  of  the  impube  component  resulting  in  upward  thrust  is: 

/usable  ■  /total  cos  a  (V  -  28) 

where  a  b  the  base  angle  of  the  truncated  conical  plug. 

In  trajectory  calculations  b  convenient  to  use  time  averaged  thrust.  Thb  mean  thrust  b 
defined  as: 

T.jJLr"*/^^  (V  -  26) 

‘CYCLE  Jt  ‘CYCLE 

where  I  b  the  impube  delivered  per  laser  pube,  and  b  the  blast  wave  period.  The  blast  wave 
period  b  simply  the  total  expansion  time  of  the  LSD  blast  wave,  plus  the  refresh  time;  as  gtven  by: 

f  cycle  =  =  <o  +  (refresh  (V  -  27) 

The  refresh  time  b  the  time  necessary  for  the  hot  expanded  gas  to  clear  the  thruster  surface  and 
be  replaced  by  the  ambient  air.  Currently,  an  adequate  model  for  the  refresh  process  does  not  exbt 
and  in  thb  analysb  the  refresh  time  b  set  equal  to  aero. 


Ill 


(V  -28) 


Input  energy  p«r  laser  pula*  is  mloiUd  as  follows: 


In  ths  derivation  of  this  equation  it  is  assumed  that  the  laser  beam  has  a  semicircular  shape  and 
uniform  Intensity.  Recall  that  the  laser  energy  is  absorbed  in  a  thin  region  behind  the  detonation 
front.  Thus,  tf  is  the  time  required  for  this  LSD  wave  to  transverse  the  ERH  thruster  plate.  To 
convert  the  pulse  energy  to  time  average  laser  power  the  following  equation  is  used: 


Finally,  an  important  measure  of  thruster  efficiency  is  the  coupling  coefficient,  which  is  the  ratio  of 
usable  single-  pulse  impulse  to  laser  input  energy  per  pulse.  Therefore,  the  coupling  coefficient  is 
defined  by: 

CC  *  a  (V  -  30) 
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D.  DYNAMIC  THRUSTER  MODEL 

The  above  analytical  model  does  not  consider  the  effects  of  flight  platform  velocity  on  ERH 
thruster  performance.  This  model  Is  sufficient  for  static  and  low  subsonic  flight  velocities,  but  when 
flight  velocities  are  supersonic,  ram  drag,  compreanbUty  effects  and  flow  expansion  over  the  aft  plug 
nozzle  (the  thruster  surface)  become  important.  For  an  airbreathing  engine  net  thrust  is  the  force 
resulting  from  the  addition  of  all  pressure  and  viscous  forces,  excluding  the  external  drag  forces; 
that  is,  gross  thrust  minus  ram  drag.  Note  external  drag  forces  caused  by  the  "installed"  engine  are 
typically  included  with  vehicle  drag.  As  a  result,  it  wu  necessary  to  develop  a  dynamic  model  for 
ERH  thruster  performance. 

In  supersonic  flight,  an  attached  oblique  shock  forms  at  the  conical  spike  tip.  The  air  is 
compressed  after  passing  through  this  oblique  shock,  and  it  is  turned  to  flow  tangent  to  the  forebody 
surface;  the  entire  forebody  of  the  Sight  platform  is  essentially  an  isentropic  spike  inlet.  The  annular 
cowl,  which  circumscribes  the  craft's  midsection,  then  redirects  the  flow  parallel  to  the  thruster 
surface.  In  the  ERH  thruster  operation,  the  cowl  only  turns  the  flow,  and  does  not  decelerate  or 
accelerate  it.  Additionally,  the  cowl  is  assumed  to  have  a  sufficiently  large  capture  area  such  that  no 
spillage  occurs  beyond  the  design  mach  number  of  three.  With  these  cowl  characteristics,  no  normal 
shock  forms  at  the  cowl  entrance  and  the  exit  mach  number  is  equal  to  that  at  the  cowl  entrance. 
Upon  leaving  the  cowl,  the  ducted  air  expands  along  the  thruster  surface  as  if  It  where  emerging 
from  a  free  expansion  type  plug  nozzle. 
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Therefor*  the  only  stagnation  significant  preesure  low  occtus  as  the  sir  pum  through  tho 
oblique  shock.  To  calculate  such  losses,  an  experimentally  determined  preaeure  recovery  inlet  sched¬ 
ule  developed  by  Marquart(29]  for  the  XRJ59-MA-3  inlet  wu  used.  In  the  present  analysis,  it  is  also 
assumed  that  the  air  flow  through  the  cowl  remains  undisturbed  by  the  LSD  wave  propagation  and 
subsequent  blast  wave  expansion.  Such  an  assumption  appears  to  be  justified,  when  the  maximum 
expansion  of  the  blast  wave  is  small  compared  to  cross  sectional  flow  area  of  the  cowL 

In  supersonic  flight,  pressure  variations  along  the  thruster  surface  become  pronounced  and 
must  be  considered  when  calculating  the  LSD  blut  wave  expansion.  When  calculating  this  pressure 
variation,  it  is  assumed  that  the  compressed  air  flow  leaving  the  cowl  expands  as  shown  in  Fig.  V-6. 
To  calculate  the  pressure  for  a  given  thruster  surface  location  (s)  it  is  necessary  to  determine  the 
local  cross  sectional  area  of  the  expanding  gases. 

Typically,  for  a  non-truncated  plug  nozzle,  the  pressure  at  the  end  of  the  plug  is  always 
equal  to  J%,  the  ambient  pressure.  Using  this  condition  and  assi uming  isentropie  flow  over  the  plug 
the  following  expression  can  be  solved  for  the  maeh  number  at  the  plug  end  (JW«). 

(V  -  31) 


where  Pu»  denotes  the  stagnation  pressure  at  the  cowl  exit.  As  previously  mentioned  this  is  deter¬ 
mine  using  the  Marquardt  inlet  data.  Once  Af«  is  known,  the  croas-sectional  flow  area,  A,,  can  be 
obtained  using: 


(V  -32) 


where  A*  and  are  the  cross-sectional  flow  area  and  the  mach  number  at  the  cowl  exit  respec¬ 
tively. 


Prom  the  nozzle's  geometry: 


(V  -33) 


If  the  croas-sectional  flow  area  is  assumed  to  vary  linearly  along  the  thruster  surface,  then: 


A  =  rjK*  -  (— )] 
1  ^  'cosor1 


(V-34) 


In  this  equation,  a  is  the  plug  base  angle  and  s  denotes  an  arbitrary  position  along  the  plug  surface. 
For  a  given  s,  the  corresponding  mach  number  can  be  found  implicitly  using: 


m{$) = i^-S-y3 


(V  -36) 


Now  the  local  expansion  pressure  can  be  calculated  from: 


Pfsl  - 

1  (l+2ylA/(*)»)^ 


(V  -  36) 
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Li  til*  dynamic  model,  the  LSD  bloat  waveamuet  expand  to  thia  praaanra. 

To  calculate  the  net  thrust,  a  control  volume  is  constructed  as  shown  in  Flf.  V-6.  By 
choodng  the  control  volume  boundary  adjacent  to  ERA  thruster  surface,  the  complex  flow  conditions 
due  to  the  interaction  between  the  blast  wave  expansion  and  the  air  flow  Greeted  by  the  cowl  can 
be  neglected;  this  simplifies  the  drag  analysis. 

Recall  that  in  the  development  of  the  LSD  blast  wave  model,  the  ambient  pressure  was 
already  subtracted  from  the  integrated  Impulse,  hence  it  is  necessary  only  to  consider  the  pressure 
difference  across  the  cowL  The  vertide  force  due  to  this  pressure  difference  is  given  by: 

Pressure  Drag  *  (Po  -  P««)Am  sin  a  (V  -  37) 

Note  that  there  is  a  force  due  to  the  momentum  change  of  the  incoming  air  flow  is  redirected  to 
convect  over  the  thruster  surface.  This  force  is  calculated  using  the  following  expression: 

Cowl  Drag  =  m^V^sm  a  {V  -  38) 


Therefore  the  net  thrust  is  given  by: 

7m»  *T- Pressure  Drag -Cowl  Drag  [V  -  39) 


E.  RESULTS 

Performance  maps  for  the  ERH  thruster  engine  were  generated  by  incorporating  the  pre¬ 
ceding  static  and  dynamic  models  into  a  comprehensive  computer  simulation.  Using  this  computer 
model,  the  effects  of  flight  mach  number  and  altitude  on  net  engine  thrust,  as  well  as  several  other 
performance  parameters  was  investigated. 

As  previously  mentioned,  the  reference  ERH  thruster  surface  configuration,  is  a  truncated 
conical  plug  which  comprises  the  entire  aft  section  of  the  vehicle.  This  plug  has  a  semi-vertex  angle 
of  45  degree*,  and  numerous  LSD  waves  are  ignited,  equally  spaced,  around  the  circumference  of 
the  plug.  The  exact  number  of  LSD  waves  which  were  studied  ranced  from  12  to  48.  In  each  caae, 
the  laaer  spot  diameter  was  varied,  such  that  the  peak  laser  power  supplied  to  the  ERH  thruster 
remained  constant.  For  this  analysis,  the  beam  intensity  used  to  initiate  the  LSD  waves  wss  held 
constant  at  6x10*  Wfcni *,  over  the  semi-circular  spot  of  0.5  cm  or  1.0  cm  for  the  48  and  12  line 
systems,  respectively. 

Before  time  avenge  performance  parameters  can  be  calculated,  it  is  first  necessary  to  deter¬ 
mine  the  total  thruster  cycle  time.  Because  the  refresh  time  has  been  neglected,  the  total  cycle  time 
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ia  ilitwmimi  by  the  maximum  npurioo  time;  Le^  the  tim*  mewiwy  (or  the  entire  bloat  wav* 
to  expand  to  ambient  pressure.  At  low  fight  velocities,  the  local  praaaure  over  the  conical  nozzle 
(the  thruster  surface)  ia  roughly  uniform,  and  the  blast  wave  expansion  time  remains  constant  over 
the  entile  thruster  area.  For  supersonic  fight  velocities,  however,  high  pressure  air  is  eonveeted 
from  the  cowl  across  the  thruster  surface.  Therefore,  the  local  ambient  pressure  at  the  base  of  the 
conical  thruster  surface  is  neater  than  that  at  its  truncated  tip.  U  the  difference  la  the  blast  wave 
expansion  times  between  the  base  and  tip  of  the  conical  plug  is  greater  than  the  transit  time  of  the 
LSD  wave  front  (this  is  always  the  case),  the  max  expansion  time  is  expected  to  occur  at  the  lower 
thruster/ noule  end  (where  the  LSD  waves  are  initiated).  Fig.  V-7  shows  the  max  expansion  time 
and  expansion  radius  for  a  thruster  at  an  altitude  of  10km  and  a  fight  mach  number  of  3.  Notice 
that  the  max  expansion  time  occurs  near  the  lower  end  of  the  plug  nossie. 

In  Fig.  V-8  net  thrust  variation  is  presented  as  a  function  of  Sight  mach  number  and 
altitude  for  the  48  beam  ERH  thruster  configuration.  Recall  that  net  thrust  is  the  time-averaged 
thrust  minus  ram  drag.  For  a  given  altitude,  observe  that  thrust  drops  off  at  some  critical  mach 
number.  This  is  due,  in  part,  to  the  increased  ram  drag  experienced  by  the  ERH  thruster  engine  at 
higher  mach  numbers. 

Fig.  V-8  also  shows  that  mean  thrust  decreases  with  increasing  altitude.  This  result  is  not 
unexpected  since  at  higher  altitudes  atmospheric  pressure  and  density  decrease.  The  ERH  thruster, 
like  other  sirbresthing  engines,  relies  on  momentum  transfer  of  incoming  ambient  sir  to  generate 
thrust;  therefore,  as  the  density  decreases,  so  does  engine  thrust. 

Another  factor  which  causes  a  reduction  in  net  thrust,  for  this  ERH  engine  configuration, 
is  illustrated  in  Fig  V-9:  the  maximum  pulse  repetition  frequency  (PRF).  Note  that  PRF  decreases 
with  increasing  altitude.  During  ERH  thruster  operation,  LSD  blast  wave  pressure  must  expand  to 
local  ambient  levels  before  the  thruster  surface  can  be  refreshed  and  new  LSD  waves  initiated.  At 
higher  altitudes,  the  local  static  pressure  is  lower,  so  the  blast  wave  expansion  time  (to)  is  generally 
longer.  Since  the  high  internal  pressure  of  the  plasma  blast  wave  decreases  rapidly  as  the  wave 
expands,  most  of  the  engine  thrust  is  generated  during  the  Initial  blast  wave  decay.  Therefore,  an 
increase  in  the  LSD  blast  wave  expansion  time  (to)  does  not  greatly  Increase  impulse  delivery  to  the 
thruster  plate,  and  time-  averaged  thrust  must  decrease. 

Closely  related  to  the  PRF  »  the  time-averaged  laser  power;  its  variation  with  altitude  and 
flight  mach  number  can  be  seen  in  Fig.  V-10.  Also,  in  Fig.  V-ll,  the  predicted  coupling  coefficient 
(cc)  performance  is  plotted,  with  the  dotted  line  indicating  the  theoretical  limit  for  airbroathing 
engines.  Since  the  calculated  coefficients  at  high  altitudes  (e.g.,  30km)  so  closely  approach  the 
theoretical  limit,  it  is  suspected  that  this  first  order  analysis  of  ERH  engine  performance  may  be 
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The  pdbntuM  maps  presented  heretofore  an  for  the  48  beam  ERH  thruster  configuration. 
Comparison  of  these  with  thoae  for  a  12  beam  thruster  configuration  (Figs.  V-12,  IS,  14  and  16) 
shows  that  for  a  given  operating  condition,  the  48  beam  configuration  fhree  better  throat  and  ita  PRF 
is  greater.  Therefore,  while  the  laaer  energy  per  pulse  is  the  same  in  both  configurations,  the  time- 
averaged  power  for  the  12  beam  configuration  is  less.  Higher  time-average  power  produces  increased 
thrust  for  the  48  beam  ERH  thruster.  Although,  time-averaged  thrust  increases  with  number  of 
beams  for  this  engine  concept,  then  obviously  is  a  limit  to  the  minimum  beam  diameter;  e.g.,  u  set 
by  optical  parameters.  This  trend  should  not  be  generalized  to  all  ERH  thruster  configurations  at 
this  point.  Each  thruster  configuration  should  be  analysed  as  a  separate  entity. 

E.  Future  Extensions 

Then  an  a  number  of  reasons  for  the  high  values  of  predicted  performance  for  the  ERH 
thruster.  In  the  current  model,  as  mentioned  above,  the  refresh  time  (tnitMh)  is  set  to  zero.  In 
an  actual  engine,  tncrwi  would  have  some  finite  value,  which  would  cause  a  decrease  in  the  PRF 
and  hence,  the  time  averaged  thrust.  A  second  model  simplification,  which  probably  contribute*  to 
producing  optimistic  results,  is  the  omission  of  axial  rarefaction  waves  within  the  LSD  blast  wave 
simulation.  Inclusion  of  axial  expansion  effects  in  this  blast  wave  model,  would  result  in  a  reduction 
in  the  impulse  delivery  to  the  thruster  surface.  At  the  present  time,  the  combined  effects  of  these 
simplifications  on  ERH  thruster  performance  have  yet  to  be  quantified.  Future  development  of  these 
ERH  thruster  models  should  address  the  above  deficiencies. 
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CHAPTER  VI 


LSC  WAVE  ERH  THRUSTER  MODEL  (STATIC) 

A  pubed  LSC  wave,  ipecifle  to  ERH  thruster  impulae  generation,  b  characterized  by  high 
plums  temperatures  0(11000  K)  sod  moderate  pressures  (  according  to  the  intend ty),  sad  low  to 
moderate  propagation  velocities.  The  high  pressures  are  somewhat  uncharacteristic  of  the  usual 
combustion  process.  The  term  Laser  Supported  Detonation  (LSD)  waveb  used  because  it  hu  the 
characteristic  high  pressures  of  the  upper  end  of  a  pressure  vs.  specific  volume  (inverse  of  density) 
plot  known  u  a  Chapman-  Jouguet  (C  J)  diagram  based  on  their  combustion  theory  ( e.g.,  see  Holmes 
et«L)|8j, 

This  point  b  an  upper  C-J  point,  the  term  LSC  implies  a  lower  pressure,  weaker  wave,  le., 
a  deflagration  wave  which  it  really  isn’t.  C-J  theory  assumes  instantaneous  energy  deposition,  not 
continuous  to  some  finite  pube  time  which  known  LSC  wave  theories  assume.  Theee  may  be  why 
some  prefer  the  term  Laser  Sustained  Pluma  (LSP)  instead  of  LSD  or  LSC. 

An  LSC  wave  eodsta  between  laser  intensities  of  about  2.5  X  104  to  7  x  10*  W far?  at  a 
wavelength  of  10.6  microns.  The  low  end  of  thb  range  represents  the  thruhold  for  air  pluma 
generation  with  a  laser  wavelength  of  10.6  microns,  whereu  the  upper  limit  b  the  transition  region 
between  LSC  and  LSD  wave  propagation  for  that  same  wavelength.  The  lower  intensity  required  for 
LSC  wave  generation  maku  it  more  attractive  than  the  LSD  wave  from  the  standpoint  of  reduced 
noise  during  hover  and  VTOL  operations. 

The  baric  elements  of  an  LSC  wave  are  a  hot  pluma  none,  the  moving  rarefaction  fans 
produced  in  the  pluma  zone,  and  the  shock  outride  of  the  pluma  (Fig.  VM).  The  physics  of  LSC 
are  not  u  well  understood  u  that  of  LSD.  An  attempt  wu  made  to  use  Sedov’s  blast  wave  scaling 
laws  to  model  the  time  history  of  the  LSC  wave  pressure  and  physical  dimension  (e.g.  pluma  cone 
radius)  in  predicting  its  effectiveness  for  low-aititude  ERH  thruster  propulrion.  Sedov's  theory  u 
also  meant  for  instantaneous  energy  deposition.  It  hu  been  used  in  a  modified  form  to  include  some 
cases  of  linear  energy  deposition  (with  time)  which  b  the  form  used  fee  a  wave  powered  by  a  laser, 
while  its  original  form  is  used  when  the  laser  pube  hu  ended  (i.e.,  unpowered). 

An  LSC  wave,  being  a  luer  generated  and  sustained  pluma,  should  supposedly  be  optically 
thin  so  that  the  laser  radiation  diffuses  throughout  the  pluma  for  its  initiation,  growth  and  support 
by  the  luer  energy.  It  hu  yet  to  be  verified,  however. 

A.  Description  of  the  Model 
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In  the  fin*  baton*  of  the  process,  laser  Ugh*  b  focwd  a*  tha  can  tar  of  tha  cnfk’a  lowar 
surface,  forming  a  cylindrical  plasma  tons  which  is  assumed  to  expand  radially  at  a  velocity,  Vise, 
given  by  the  foUowing  expression: 

Vuo  ~  C' +  '>1(7,  +  iIi+'&'m- HV)1* 

where  7  and  tc  are  the  ratios  of  specific  heats  In  the  plasma  and  the  ambient  air  respectively,  • 

is  the  received  laser  intensity  (assumed  85air,  and  W  is  the  percentage  difference  between  the  wave 

velocity  above  and  the  particle  velocity  through  the  LSC  wave.  Pirri  et  aL  give  0.04  as  a  good 

value  for  W.  7  »  1.4  and  7  **  1.2  are  taken  for  these  respective  parameters.  The  average  radius  of 
0 

the  cylinder  may  be  approximated  by  multiplying  this  velocity  by  the  time: 

H-Vtsct  (V  /  -  2) 

During  this  phase,  the  impulse  delivered  to  the  craft  is  due  to  the  pressure  inside  of  the 
plasma  cone.  The  initial  plasma  cone  pressure  is  given  by: 

- 11  -  («"-’) 

This  LSC  pressure  equation  was  first  derived  by  Pirri  at  sl|16]and  has  been  assumed  for  this 
analysis.  However,  it  was  derived  for  the  case  of  a  laser  beam  perpendicularly  incident  to  a  surface, 
not  transversely.  This  equation  (plotted  in  Fig.  VI-2),  gives  rather  large  values  of  the  preoure  at 
the  higher  intensities  which  as  previously  stated  may  be  the  nature  of  an  LSC  wave  and/or  possibly 
indicative  of  the  transition  to  an  LSD  wave  or  because  of  the  different  boundary  conditions. 

As  the  transition  occurs,  the  value  of  W  should  increase  to  that  for  an  LSD  wave,  according 
to  Pirri  et  aL  These  LSC/LSD  transition  effects  are  not  considered  in  this  analysis.  However,  higher 
’intensity  calculations  have  been  made,  despite  this,  to  observe  the  effect. 

This  LSC  pressure  is  the  reference  pressure  used  initially  in  the  appropriate  scaling  law 
for  the  changes  In  the  Internal  plasma  pressure  in  this  analysis.  Once  the  LSC  wave  is  initisted, 
a  rarefaction  fan  travels  from  the  top  of  the  plaama  tone  to  the  surface  of  the  craft  at  the  plasma 
sound  speed  given  by: 


(V/-4) 

Once  the  rarefaction  fan  has  reached  the  surface,  the  pressure  inside  of  the  plasma  tone 
begins  to  decrease  rapidly  even  with  the  laser  on.  ft  was  concluded  that  the  beet  usage  of  laser 
power  would  result  by  setting  the  laser  pulse  time,  t, ,  equal  to  the  plasma  tone  height  divided  by 
plasma  sound  speed.  Hence: 


12S 
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(IV  -5) 


This  would  increase  the  laser  pulse  repetition  rate  and  improve  the  propulsion  system’s 
coupling  coeffldent  (i.e.  thrust  per  input  power). 

The  casesstudied  utilized  cylinder  heights  of  6, 12.5, 25, 50,  and  100  cm.  The  sir  inside  of  the 
cylindrical  plasma  zone  can  reach  temperatures  in  excess  of  15,000  K,  well  beyond  the  temperature 
limit  of  any  known  materials.  Because  of  the  ehort  pulse  duration,  it  was  assumed  that  the  surface 
(either  by  active  or  passive  cooling)  can  maintain  its  structural  integrity  throughout  all  operational 
modes.  The  following  section  describes  the  various  scaling  laws  used  and  the  physical  significance 
of  each  one. 

During  the  pulse  time,  the  pressure  inside  of  the  plasma  zone  is  constantly  changing  with 
time  due  to  the  rarefaction  wave  and  the  plasma  zone’s  radial  expansion.  The  scaling  laws  as 
presented  by  Sedov  were  used  to  account  for  this  behavior.  Because  of  the  actual  physical  geometry 
of  the  plasma  region,  the  unpowered  cylindrical  scaling  law  wu  used  with  all  reference  parameters 
taken  at  the  end  of  the  pulse  duration  ( i.e.,  t,e/  =  tf,pr*/  =  Vise  )•  This  is  given  by: 


Vt,i 


(VI  -6) 


Thus,  once  the  rarefaction  wave  reaches  the  surface  and  the  laser  beam  is  shut  off.  The 
plasma  zone  then  decays  in  a  manner  described  by  the  following  relationship: 


jL  =  (_L)f 

Sr*/  tr«/ 


(VI  -  7) 


Two  conditions  were  examined.  The  first  condition  assumed  that  the  plasma  zone  would 
continue  to  expand  in  a  manner  which  could  still  be  modeled  with  the  unpowered  cylindrical  scaling. 
The  expansion  occurs  to  ambient  pressure  with  the  reference  radius  smaller  than  the  plasma  zone 
height. 


In  the  second  condition,  the  plasma  zone  needs  to  further  expand  even  when  its  radial 
dimension  is  equal  to  its  height.  Because  of  the  physical  geometry  involved  (fig.  VM),  unpowered 
cylindrical  scaling  would  not  be  an  appropriate  scaling  law.  The  additional  expansion  is  described 
by  the  unpowered  spherical  scaling  law. 


Prtf  fr  */  ’  Rf*f  U  tj 


(VI-  8) 


Full  expansion  to  ambient  is  described  with  this  law.  Once  the  plasma  zone  has  expanded 
to  ambient  (by  either  condition),  air  is  assumed  to  refresh  the  lower  surface  at  about  perhaps  the 
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craft  speed.  One*  new  air  in  in  place,  the  cycle  b  really  to  begin  again. 

A  computer  program  war  written  to  model  the  proposed  cycle.  Laser  intensity  and  cone 
height  are  entered  into  the  computer  velocities,  pressures  and  time  scales  are  all  computed.  Ambient 
data  were  all  taken  for  sea  level  conditions.  The  program  checks  whether  condition  one  or  condition 
two  will  prevail  and  run  through  the  appropriate  scaling  laws.  Reference  values  are  also  generated 
by  the  computer  program.  The  program  calculates  the  total  Impulse  delivered  to  the  craft  by  the 
equation: 

/-  r  (*'/-») 

JO  JA(  0 

where  to  b  the  time  it  takes  for  the  pressure  to  decay  to  the  ambient  value  po-  to  and  Rm*i  are 
calculated  by  the  scaling  laws.  The  program  calculates  Impulse  by  dividing  the  integral  into  domains 
of  separate  scaling  laws  and  then  adds  these  domains  to  obtain  the  total  impulse. 

The  impulse  for  each  phase  (or  domain)  and  each  condition  b  then: 
unnowered  cylindrical: 


/  *  *(Vcsctr )2(^Wc(tr  “  t»)  ~ 

(V7-10) 

or  unnowered  spherical: 

/  =  -  *1)  -  §P0tK*(t#  - t|) 

(VI  - 10) 

where:  /to  =  Vtsct,(^)* 

P *  *  Ptso(p)** 

-  (£>-•*« 

The  total  cycle  period  b  equal  to  to  plus  the  refresh  time  if  any  significant  refresh  time 
exbts.  Thb  period  b  used  to  get  the  maximum  repetition  rate  of  the  laser  pulse: 

PRF  =  i  (VI-  12) 

Then  the  mean  thrust  may  be  calculated  by: 

T  =  I  x  PRF 
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(VI  - 13) 


Time-avenged  power  k  calculated  by  multiplying  the  incident  laser  intensity  by  the  average 
puked-time  LSG  wave  surface  area. 

P  =  2»R(to)f  (V/-14) 

The  puke  energy  k  calculated  ic 

E,  =  2t  J'  *VCsctHdt  a  t Visctjff *  (V/  -  15) 

The  coupling  coefficient  k  calculated  in  the  usual  manner  by  dividing  the  impulse  by  this 

energy. 


ce 


(V7-16) 


B.  Analytical  Results 

The  results  obtained  from  the  LSC  wave  computer  model  are  illustrated  in  Figures  Vl-3 
through  VI-10.  Figure  VI-3  plots  mean-thrust  versus  intensity,  with  6  curves  corresponding  to  5 
different  plasma  rone  cylinder  heights.  The  mean  thrust  increases  with  the  height  and  the  intensity. 
As  might  be  expected,  a  larger  plasma  rone  height  means  that  the  larger  amount  of  energy  being 
deposited  into  the  thruster  k  thus  producing  a  larger  thrust. 

As  shown  in  Figure  VI-4,  the  PRF  decreases  with  increasing  intensity  at  a  given  plasma 
rone  height.  At  a  given  intensity,  the  PRF  decreases  with  increasing  plasma  none  height,  since  the 
time  required  for  the  rarefaction  wave  to  traverse  from  top  to  bottom  k  longer,  and  the  velocity  of 
the  wave  remains  unchanged.  It  k  noted  that  the  PRF  represented  here  may  be  the  fastest  allowable 
cycle  time  if  the  refresh  time  was  negligible.  By  rotating  the  thruster  surface  or  the  incident  laser 
fluxes  maintaining  a  constant  supply  of  fresh  air  and  area  for  the  LSC  wave,  the  cycle  could  be 
quickly  be  initiated  again  as  soon  as  ambient  conditions  are  reached. 

The  plot  of  coupling  coefficient  versus  intensity  shown  In  figure  VI-5  reveak  some  interesting 
features.  It  does  display  an  extremum,  which  locates  an  optimal  coupling  coefficient  at  about 
10*W  I  cm2,  the  step  in  the  curve  occurs  as  the  program  transitions  from  condition  1  to  condition 
2.  Thk  dip  along  with  other  sudden  changes  in  C.C.,  PRF,  T  and  ro,  may  be  appearing  because 
the  transition  to  LSD  conditions,  which  k  a  complicated  process,  k  occurring.  The  implication  k 
that  condition  2  k  being  chosen  at  thk  intensity  because  the  wave  k  stronger  and  k  reaching  H 
sooner.  Thk  effect  appears  at  later  time  in  the  cycle  implying  that  more  time  k  needed  for  full 
wave  expansion.  Thk  may  support  the  above  hypothesis.  Thk  verifies  that  the  higher  intensities 
produces  higher  coupling  coefficients  due  to  a  longer  thrust  schedule. 

Figure  VI-6  shows  how  much  energy  k  being  accepted  by  the  LSC  wave  from  the  amount 
being  dumped  into  it.  It  also  increases  with  height  and  intensity. 
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The  kirn*  averaged  power  should  show  s  similar  trend  if  than  its  no  lososs  in  the  cyclo. 
However,  around  10 *IV/em*.  figure  VI-7  shows  a  dip  in  the  power  required  by  the  wave  for  the 
cycle. 

Figure  VI-8  simply  shows  that  as  more  energy  is  being  deposited  into  the  wave,  It  takes 
less  time  for  the  wave  to  grow,  hence  the  decrease  of  the  pulse  time  with  the  intensity.  Figure  VI-9 
shows  the  growth  of  the  LSC  wave  within  that  time.  It  depends  on  the  gis  properties  and  not  the 
flux  because  they  should  be  the  ultimate  constraint  on  the  gas  geometry.  This  is  a  nature  of  the  gas 
dynamics  reflected  in  the  equations  (mainly  Sedov’s)  and  in  these  results. 

Figure  VI- 10  shows  the  total  growth  of  the  wave  and  displays  the  effect  mentioned  earlier 
as  it  appears  in  the  complete  wave  expansion. 

C.  Summary 

This  analysis  was  the  remit  of  a  simple  model  of  an  LSC  wave  to  be  used  in  an  ERH  thruster. 
The  results  appear  promising.  A  considerable  amount  of  thrust  seems  achievable  for  hovering  and 
mon.  The  ERH  thruster  itself  appears  to  compare  very  well  with  known  engines  (compare  figure 
VI-5  with  the  40  N/MW  capability  of  existing  ion  thrusters).  However  some  of  the  assumptions 
made  impose  some  restrictions  on  these  results.  The  sise,  shape  and  strength  of  the  wave  may  not 
agree  with  the  use  of  Sedov’s  scaling  laws.  It  is  apparent  that  Sedov’s  scaling  laws  are  questionable 
at  r=t0  and  Infinity.  For  the  low  velocities  of  the  LSC  wave  and  its  high  temperatures,  thermal  and 
viscous  losses  may  be  very  important  and  would  have  to  be  considered  in  further  analyses. 

The  geometry  of  the  LSC  wave  itself  may  even  be  more  complicated  then  it  has  been 
imagined  to  actually  be  in  figure  VI-1.  The  model  used  in  this  analysis  did  not  fully  acknowledge 
the  complexities  of  figure  VM,  though  it  is  believed  to  be  adequate  in  this  first  investigation. 

Pirri’s  pressure  and  velocity  equations  were  derived  for  a  l-D  model  and  was  good  to 
start  with,  but  their  compatibility  with  the  boundary  conditions  of  this  configuration  should  be 
re-analysed  in  the  future.  Pirri  had  also  assumed  that  a  strong  shock  existed  in  front  of  the  LSC 
wave.  This  simplified  his  analysis  but  yielded  those  rather  high  pressures  for  the  wave. 

An  acoustical  analysis  appears  to  be  the  next  best  step  in  the  development  of  any  future 
models.  It  would  be  even  better  if  a  numerical  model  could  be  developed  using  finite  difference 
methods. 
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CHAPTER  VO 


TEMPERATURE  HISTORY  OP  AN  ERH  THRUSTER  SURFACE 

An  ERH  thruster  surface  that  is  exposed  to  a  laser-  supported  combustion  (LSC)  wave  can 
experience  high  temperature  due  to  radiative  heat  transfer  between  the  wave  and  the  surface.  As 
the  wave  propagates  periodically  over  the  surface,  the  surface  temperatures  increase  on  each  cycle. 
An  initially  cool  surface  will,  after  a  time,  attain  the  maximum  allowable  temperature  for  the  surface 
material.  The  time  to  maximum  temperature  can  be  estimated  by  calculating  the  net  radiative  heat 
transfer  to  the  surface  coupled  with  the  transient  thermal  response  of  the  surface. 

The  vehicle  configuration  examined  in  this  study  utilises  an  ERH  thruster,  which  is  sketched 
in  Fig.  VH-1  [10].  The  baa*  of  the  vehicle  is  a  disk  of  radius  1.25  m.  Energy  is  supplied  to  the 
LSC  wave  by  a  laser  beam  that  is  projected  inward  from  the  edges  of  the  disk  (around  the  entire 
circumference).  The  LSC  wave  propagates  radially  outwards  towards  the  edge  of  the  disk.  At  the 
edge,  the  beam  is  extinguished,  and  a  new  wave  is  initiated  at  the  disk  center.  This  process  continues 
periodically  until  adequate  thrust  is  achieved  or  until  maximum  temperatures  are  exceeded.  This 
chapter  analyses  the  static  “hover"  or  subsonic  flight  condition  only,  where  forced  convective  cooling 
of  the  ERH  thruster  surface  cannot  be  invoked.  Forced  convective  heating  of  the  surface  was  not 
analysed. 

Figure  VH-2  shows  an  alternate  view  of  the  propagating  LSC  wave.  As  the  wave  expands 
outward,  a  rarefaction  wave  propagates  inside  the  LSC  wave  cylinder.  The  pressure  difference 
between  the  back  of  the  rarefaction  wave  and  the  lower  vehicle  surface  produces  the  flight  propulsive 
thrust.  The  physical  rise  of  the  two  LSC  wave  configurations  which  are  examined  in  this  analysis  were 
determined  by  the  hovering  thrust  which  is  required  by  a  vehicle  the  sise  of  a  Mercury  capsule.  The 
two  LSC  waves  have  the  requisite  heights  and  pulse  rate  frequencies  to  produce  identical  thrust.  The 
difference  between  these  two  configurations  arises  from  the  assumed  emisrivity  of  the  ERH  thruster 
surface.  The  first  LSC  wave  propagates  over  a  surface  with  an  emisrivity  of  0.1  while  the  second 
has  an  emisrivity  of  0.8.  Note  that  both  ERH  thruster  surfaces  are  basically  horizontal  (l.e.,  a  flat- 
bottomed  vehicle). 

A.  Radiation  Model 

The  LSC  wave  is  essentially  a  thin  layer  of  absorbing  and  emitting  gas  at  an  elevated 
temperature.  In  the  present  study,  the  LSC  wave  is  assumed  to  be  gray  and  isotropic  with  respect 
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to  temperature,  pronto  and  radiative  properties.  No  scattering  of  radiation  ia  considered.  Tho 
thickneaa  of  the  LSC  wave  is  small  relative  to  the  diameter  of  the  ERH  thruster  surface  and  the 
height  of  the  wave.  This  fact  allows  a  significant  simplification  in  the  radiation  modeling.  Energy 
emitted  from  the  LSC  wave  that  arrives  at  a  point  on  the  thruster  surface  will,  in  general,  have 
passed  through  a  short  path  in  the  LSC  wave  and  a  long  path  in  the  transparent  medium  (air) 
between  the  LSC  wave  and  the  surface.  If  all  the  radiation  emitted  by  the  LSC  wave  is  assumed  to 
originate  at  a  plane  in  the  center  of  the  wave,  the  total  path  length  and  orientation  of  the  emitted 
radiation  will  not  be  greatly  in  error.  The  LSC  wave  can  then  be  replaced  by  an  infinitesimally  thin 
LSC  nrftct  which  emits  the  same  total  radiative  flux  as  the  actual  LSC  wave.  To  determine  the 
radiative  flux  emitted  by  the  LSC  wave,  a  mean  beam  length  analysis  [23]w as  used.  The  mean  beam 
length  for  two  infinite  parallel  planes  (the  inner  and  outer  diameter  of  the  LSC  wave), 

L.s-ilni2£a(«,0)l  (V//-1) 

•b 

where  d,  is  the  Planck  mean  absorption  coefficient,  D  is  the  thickness  of  the  LSC  wave,  and  E»  is 
an  exponential  integraL  The  quantity  a ,D  is  the  optical  thickness  of  the  LSC  wave.  For  the  cases 
investigated  here,  the  LSC  temperature  is  16,000  K,  the  pressure  is  3  atm,  and  the  thickness,  D, 
is  1.6  cm.  This  corresponds  to  a  laser  intensity  of  WWfm*.  The  Planck  mean  absorption 
coefficient  for  these  conditions,  from[32^  is  a,  *  0.01cm,  and  the  corresponding  optical  thickness  is 
0.015.  The  mean  beam  length  can  be  used  to  determine  the  total  emittance  of  the  LSC  wave  by: 

s  =  |l-exp(-s,L.)l  {VII-  2) 

In  the  present  study,  the  value  of  the  total  emittance  was  0.03.  The  purpose  of  the  present 
analysis  is  to  predict  temperatures  on  the  thruster  surface.  This  will  be  accomplished  by  solving 
energy  balance  equations  at  the  thruster  surface.  The  present  analysis  does  not  attempt  to  specify 
a  correct  energy  balance  for  the  LSC  wave.  The  temperature  of  the  LSC  wave  is  regarded  as  given 
from  gas  dynamic  considerations.  With  these  stipulations,  it  is  possible  to  replace  the  actual  LSC 
wave  with  an  equivalent  black  surface.  The  black  LSC  surface  emits  the  same  total  energy  as  the 
actual  LSC  wave.  Since  the  actual  LSC  wave  has  an  emittance  of  0.03  and  a  black  surface  has  an 
emittance  of  1.0,  the  black  surface  must  have  a  lower  equivalent  temperature,  specified  by 

oT*=«jTIsc  {VII-  3) 

where  T«  is  the  equivalent  temperature  and  Tisc  is  the  actual  temperature  of  the  LSC  wave.  The 
calculated  equivalent  temperature  of  the  black  LSC  surface  is  6660  K. 

This  equivalent  black  surface  will  produce  the  correct  incident  radiation  on  the  ERH  thruster 
surface.  Radiative  energy  that  originates  at  the  thruster  surface  and  is  incident  on  the  actual  LSC 
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w&v*  will  either  b«  transmitted  through  tho  waveor  be  absorbed  by  the  wave.  There  is  no  scattering 
in  the  LSC  wave,  eo  energy  originating  at  the  vehicle  aui&ce  will  not  be  scattered  back  to  the  vehicle. 
This  is  the  reason  that  the  equivalent  LSC  surface  must  be  black,  so  th&t  no  energy  leaving  the  vehicle 
surface  will  be  reflected  back  to  the  vehicle  surface.  From  the  point  of  view  of  the  vehicle  surface, 
the  actual  LSC  wave  looks  exactly  like  a  black  surface. 

Since  the  equivalent  LSC  wave  surface  has  a  lower  temperature,  the  spectrum  of  emitted 
radiation  will  be  shifted  to  longer  wavelengths  than  exist  in  the  actual  case.  This  is  not  a  matter 
of  concern,  because  the  thruster  surface  will  be  assumed  to  be  gray,  and  only  the  magnitude  of 
the  incident  flux  and  not  the  spectral  distribution  of  energy  will  affect  the  thruster  surface  energy 
balance. 

Replacing  the  actual  LSC  wave  with  an  equivalent  black  surface  will  not  lead  to  a  correct 
energy  balance  for  net  radiation  from  the  LSC  wave.  As  mentioned  earlier,  this  analysis  does  not 
attempt  to  model  the  behavior  of  the  LSC  wave,  and  our  only  concern  is  that  the  net  radiation  to 
the  thruster  surface  can  be  modeled  using  enclosure  theory  for  diffuse  gray  surfaces^}.  The  view 
factors  needed  in  the  radiation  analysis  can  be  obtained  analytically  by  using  the  results  of  Realty 
and  applying  configuration  factor  algebra. 

B.  Conduction  Model 

The  radiation  to  the  thruster  surface  is  coupled  with  conduction  in  the  surface.  Initially  the 
thruster  surface  is  at  atmospheric  temperature.  An  LSC  wave  is  Initiated  at  the  center  of  this  disk¬ 
shaped  surface,  and  propagates  outward  to  the  edges.  The  time  of  propagation  to  the  edge  is  short 
compared  to  the  thermal  time  constant  for  the  thruster  surface;  therefore,  the  surface  temperature 
is  sasumed  to  remain  constant  during  the  cycle.  After  the  first  cycle,  the  net  radiative  flux  to  the 
surface  during  the  cycle  is  calculated  u  a  function  of  radial  position.  This  radiative  flux  is  then  used 
as  input  to  a  transient  conduction  analysis  of  the  surface.  The  result  of  the  conduction  analysts  is 
an  estimate  of  the  thruster  surface  temperature  rise  during  the  first  cycle.  This  temperature  may 
vary  with  radial  position. 

The  proce«  is  then  repeated  for  the  next  cycle.  The  radiati/e  flux  to  the  surface  during 
the  second  cycle  is  calculated  is  a  function  of  the  new  surface  temperature  distribution.  At  the 
end  of  the  second  cycle,  the  conduction  analysis  is  repeated  with  a  new  initial  temperature  and 
heat  flux  distribution.  This  process  of  solving  radiation  and  conduction  alternately  is  continued  in 
each  succeeding  cycle.  It  is  possible  to  solve  coupled  radiation  and  conduction  in  this  way  without 
iterating  on  the  surface  temperatures  because  of  ths  very  short  cycle  time,  and  the  resulting  relatively 
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•low  riss  in  surface  Umpntun.  Using  values  of  tsmpsratun  from  a  provioua  tints  stop  instead  of 
current  values  does  not  result  in  *  significant  error  in  heat  flux  rates. 

To  compute  the  radiative  flux  to  the  thruster  lurfaee,  the  surface  is  divided  into  numerous 
concentric  isothermal  rings.  At  the  beginning  of  the  cycle,  the  LSC  surface  is  located  at  the  boundary 
between  the  innermost  ring  (which  is  actually  a  disk)  and  the  second  concentric  ring.  The  LSC 
surface  resembles  a  very  thin-walled  hollow  cylinder.  During  the  cycle,  the  LSC  surface  steps  toward 
the  disk  perimeter,  with  its  location  always  being  between  concentric  rings.  The  time  which  the 
LSC  wave  remains  at  each  location  is  determined  by  the  velocity  of  the  LSC  wave,  the  laser  pulse 
duration,  and  the  incremental  radial  dimension.  At  each  location  of  the  cylinder,  the  configuration 
factors  must  be  recalculated,  and  the  equations  for  radiative  energy  balance  in  an  enclosure  solved 
again. 

Conduction  in  the  ERH  thruster  surface  will  be  assumed  to  be  one-dimenaional,  with  no 
conduction  in  the  radial  or  circumferential  direction*  Since  the  problem  is  axisymmetric,  then  is 
no  temperature  variation  circumferentially;  however,  a  temperature  gradient  can  exist  in  the  radial 
direction  due  to  the  radiation.  The  one-  dimensional  asmmption  is  justified  because  the  temperature 
gradients  into  the  surface  an  expected  to  be  much  larger  than  the  radial  temperatun  gradients. 

The  one-dimensional  transient  heat  conduction  equation  was  solved  for  each  isothermal  ring 
on  the  vehicle  surface.  The  conduction  equation  was  discretized  using  an  explicit  finite- difference 
method[13).  The  explicit  method  was  superior  to  an  implicit  method  in  this  case  because  of  the 
very  short  time  steps  involved.  The  numerical  solution  was  compand  to  the  analytical  solution  for 
conduction  in  a  semi-  infinite  solid  with  a  constant  surface  heat  flux,  and  agreed  well. 

C.  Analytical  Results 

Two  LSC  wave  ERH  thruster  configurations  wen  analyzed.  Both  thrusters  utilize  &  LSC 
wave  which  is  extinguished  after  traveling  over  the  inner  part  of  the  surface.  The  LSC  wave  effective 
temperatun  is  also  the  same  for  both  designs.  The  diffennce  between  the  two  thrusters  arises  from 
the  diffennt  values  for  the  emissivity  of  the  vehicle  surfaces.  The  parameters  utilized  for  input  to 
the  numerical  solution  an  shown  in  Table  VQ-l. 

A  short  time  exits  between  the  extinction  of  one  LSC  wave  and  the  initiation  of  the  new 
wave  at  the  thruster  center.  This  time  interval  is  npresented  by  the  difference  between  the  thrust 
duration  and  the  laser  pulse  time  in  Table  1.  The  thruster  surface  material  for  the  second  case  was 
graphite,  with  emissivity  of  0.6.  For  the  first  case,  the  assumed  surface  emismvity  wu  0.1.  This 
lower  emissivity  value  could  possibly  be  produced  by  coating  the  graphite  (or  other  base  material) 
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Table  VII 

Pinmten  UtlllMd  far  Input  into  Numerical  Solution 


!■■■■■■■■■■■■■ 1  '  \  ■T fT Mi i .  11  JTn  '-HU  1 

Surface  Emissivity 

0.1  0.6 

LSC  Wave  Height 

0.5m 

Radius  of  Vehicle  Surface 

1.25m 

Radius  of  Propagation 

0.23cm 

Laser  Pulse  Time 

0.36msec 

Thrust  Pulse  Time 

0.40msec 

LSC  Wave  Thickness 

1.5cm 

LSC  Wave  Pressure 

3Atm. 

LSC  Wave  Temperature 

16000K 

LSC  Effective  Temperature 

6660K 

with  a  thin  layer  of  another  substance  with  this  radiative  property  value.  The  base  material  for  the 
first  case  wu  also  graphite. 

For  an  operating  time  of  approximately  0.3  seconds,  the  average  surface  temperature  for  the 
lower  emisslvtty  case  reached  28*<7.  The  average  temperature  for  the  higher  emMvity  case  for  the 
same  operating  time  wu  approximately  66*C.  Longer  operating  times  were  not  investigated  due 
to  the  large  computation  times  which  were  required.  Nevertheless,  it  is  seen  that  this  design  shows 
promise  for  staying  within  permissible  temperature  limits  for  a  total  time  of  a  second  or  more. 

The  radial  distribution  of  the  surface  temperature  for  the  two  emissivity  cases  at  various 
times  are  given  in  Figs.  VII-3a  and  Vn-3b,  respectively.  As  expected,  the  temperature  hu  a 
maximum  value  near  the  thruster  geometric  center,  and  tends  to  decline  near  the  outer  periphery  of 
the  disk.  The  relatively  flat  temperature  profiles  which  exist  over  most  of  the  disk  surface  provide 
justification  for  the  one-dimensional  conduction  analysis.  The  maximum  temperature  after  70  cycles 
(0.3  seconds)  for  both  designs  occurs  at  the  center.  For  the  first  design,  the  center  temperature 
reached  a  maximum  value  of  176*C.  The  second  design,  which  hu  a  higher  emissivity,  had  a  center 
temperature  of  956*(7.  If  two-dimensional  conduction  effects  were  accounted  for  in  the  thruster 
surface,  the  center  temperature  would  be  somewhat  lower.  Based  on  the  difference  in  the  maximum 
values  of  the  two  center  temperatures  of  the  first  and  second  cues,  the  first  cue  is  the  more  promising 
configuration. 

For  the  conduction  analysis,  the  thruster  surface  wu  assumed  to  be  a  semi-infinite  solid. 
The  penetration  depth  of  the  temperature  field  into  the  surface  is  shown  in  Figs.  VII-4a  and  VII-4b 
for  the  first  and  second  cues,  respectively.  Below  this  depth,  there  is  no  significant  rise  in  the  vehicle 
temperature.  At  the  largest  time  studied,  the  penetration  depth  wu  about  2  cm  for  both  designs. 
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Figure  VH-3b  RR  Vereee  T/T„  Surface  Emiaahrity  of  0.0 
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Two  portion*  on  th*  surface,  tho  cantor  and  tho  *dgs,  an  the  location*  of  maximum  and 
minimum  temperature,  respectively.  Figure*  VH-5a  and  VQ-fib  rep  recent  the  surface  temperature 
time  history  of  these  two  locations  for  the  low  and  high  emissitivity  cases.  In  both  of  these  figures, 
the  fine  which  represents  the  center  temperature  is  decreasing  with  time.  This  fact  shows  that  the 
net  heat  fiux  to  the  surface  is  decreasing  with  time  as  the  surface  temperature  increases.  The  surface 
temperature  for  both  cases  may  reach  an  equilibrium  value  if  the  computational  time  is  increased. 
These  two  figures  again  show  that  the  lower  emissivity  case  bolds  more  promise  of  surviving  a  greater 
number  of  LSC  wave  cycles. 

It  should  be  noted  that  the  final  thruster  surface  temperature  is  rather  sensitive  to  the 
assumed  values  of  LSC  wave  temperature,  pressure,  and  thickness.  Since  the  LSC  wave  is  not  truly 
isothermal  and  the  thickness  is  not  sharply  defined,  it  is  difficult  to  specify  these  parameters  exactly. 
Relatively  small  changes  in  the  assumed  average  values  of  these  parameters  can  have  a  large  effect 
on  the  final  results. 

D.  Future  Extension* 

Three  main  refinements  to  the  present  analysis  will  be  made  in  the  fixture.  The  refinements 
to  be  considered  are:  1)  relaxation  of  the  assumption  that  the  absorbtlvity  and  emissivity  of  the 
thruster  surface  are  equal,  2)  two-  dimensional  conduction  effect*  in  the  thruster  surface  (conduction 
in  the  radial  direction),  and  3)  convective  heat  transfer  effects  between  the  plasma  and  thruster 
surface. 

The  radiation  analysis  refinement  will  allow  the  thrater  surface  to  have  different  absorb* 
tivity  and  emissivity  coefficients.  The  absorbtivity  coefficient  will  be  in  the  high  temperature  solar 
range,  while  the  emissivity  coefficient  will  be  in  the  lower  temperature  infrared  range.  This  improve¬ 
ment  in  the  analysis  will  allow  the  surface  to  be  chosen  on  its  selective  radiative  characteristics. 

Two  dimensional  conduction  effects  will  be  utilised  in  the  future  to  better  approximate  the 
conduction  in  the  thruster  surface.  The  inclusion  of  conduction  in  the  radial  direction  will  lower  the 
center  temperature  and  allow  a  longer  total  thrust  time.  The  convective  heat  transfer  effects  will 
also  be  included  in  the  future  when  the  gas  dynamics  of  the  moving  uSC  wave  are  better  understood. 
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CHAPTER  Vm 

PAR-FIELD  SOUND  LEVELS  PRODUCED  BY  ERH  THRUSTER 
A.  Introduction 


The  continuous  series  of  lasers  Induced  impulses  used  as  a  means  of  propulsion  in  vehicles 
powered  by  External  Radiation  Heated  (ERH)  thrusters,  will  generate  high  sound  pressure  levels  in 
the  vicinity  of  the  launch  site.  A  rough  estimate  of  these  levels  can  be  made  under  the  assumption 
that  the  ERH  thruster-propelled  vehicle  can  be  approximated  (for  acoustical  purposes)  as  a  uniformly 
radiating  spherical  source  whose  diameter  equals  that  of  the  vehicle.  The  sound  pressure  level  at 
the  surface  of  this  hypothetical  source  will  be  assumed  to  have  an  amplitude  equal  to  the  average 
pressure  exerted  on  the  vehicle  surface  by  the  laser-induced  explosions. 

For  the  two  types  of  vehicles  considered  in  this  study  -  a  6655- kg  vehicle  6  m  in  diameter 
(Apollo  Lightcraft)  and  a  900-kg  vehicle  2.5  m  in  diameter  (Mercury  Lightcraft)-  the  average  pres¬ 
sures  are  2.775  kPa  and  1.799  kPa.  These  correspond  to  sound  pressure  levels  (re.  20/tPa)  of  162.8 
and  159.1  dB,  respectively,  in  the  immediate  vicinity  of  the  vehicles.  (Sound  pressure  levels  (SPL) 
in  decibels  (dB)  are  given  by 


SPL  =  20  iog( 


Sound  Pressure  in  Pascals 
20pFa 


! 


(VJJJ-1) 


on  this  scale  the  faintest  sound  that  can  be  heard  has  a  SPL  of  0  dB,  while  the  threshold  of  pain  Is 
at  120  dB.) 

The  SPL  of  sound  emanating  from  a  spherical  source  in  the  absence  of  boundaries  will 
decrease  with  distance  r  from  the  source  as  20  log {r/R),  where  R  is  the  radius  of  the  source  (2.5  m 
in  the  case  of  the  5555-kg  vehicle  and  1.25  m  in  the  case  of  the  900-kg  vehicle).  Therefore  the  SPL 
as  a  function  of  distance  can  be  represented  by 


SPL  =  162.8  -  20  log 

SPL  a  170.8  -  20  log  r  (VIII  -  2) 

for  the  5555-kg  vehicle,  and 

SPL  =  159.1  -20^74-  (VIII-  3) 

1.40 

SPL  =  161.0-  20  logr 

for  the  900-kg  vehicle.  These  expressions  hold  for  average  pressure  levels  generated  when  the  vehicles 
are  hovering;  if  they  are  accelerating  at  the  designed  maximum  rate  of  33  gs,  the  SPLs  will  increase 
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by  20  logSS  m  80.4 dB  at  toy  noise  measuremsnt  location.  The  following  Table  show*  the  noiao 
levels  predicted  by  the  above  equations  at  distances  ranging  from  20  m  to  20  km  during  hover 
These  figures  will  be  increased  by  30.4  dB  when  the  respective  vehicles  accelerate  at  the  design 
maximum  of  33  gs. 

Table  vm-1  —  Predicted  noise  Levels  (dB) 


Distance  r 

spl  (6666-kg  vehicle) 

spl  (900-kg  vehicle) 

20  m 

144.8 

135.0 

50  m 

136.8 

127.0 

100  m 

130.8 

121.0 

200  m 

124.8 

115.0 

500  m 

116.8 

107.0 

1  km 

110.8 

101.0 

2  km 

104.8 

95.0 

5  km 

96.8 

87.0 

10  km 

90.8 

81.0 

20  km 

84.8 

75.0 

By  way  of  comparison,  the  SATURN  V  launch  vehicle  used  in  the  APOLLO  program  was 
known  (see,  e.g.,  Peterson  and  Gross,  1974)  to  produce  SPL’s  well  in  excess  of  90  dB  at  a  distance 
of  20  km  from  the  launch  site.  Table  VO-1  shows  that  the  SPL’s  produced  by  the  ERH  thruster 
vehicles  in  hover  will  not  exceed  those  of  the  SATURN  V.  The  additional  30.4  dB  produced  when 
the  vehicles  accelerate  at  33  gs  could  raise  the  SPL’s  above  those  produced  by  the  SATURN  V, 
however. 

The  simplified  picture  of  the  noise  problem  posed  by  ERH  thrusters  in  the  above  discussion 
is  rather  too  crude  to  be  used  as  a  basis  for  a  definitive  estimate  of  the  noise  levels  to  be  expected 
during  hover  and  acceleration.  A  more  refined  picture  will  be  developed  below  from  which  such  an 
estimate  can  be  obtained.  The  amplified  picture  is  useful,  nonetheless,  because  it  demonstrates  that 
the  magnitude  of  the  noise  problem  posed  by  the  launch  of  ERH  thruster  propelled  vehicles  is  not 
so  small  that  it  can  be  dismissed  out  of  hand. 

B.  BLAST  WAVE  CHARACTERISTICS 

A  realistic  acoustic  model  of  an  ERH  thruster  vehicle  must  incorporate  the  nonlinear  char¬ 
acteristics  of  blast  waves,  whose  attenuation  with  distance  are  significantly  different  from  those  of 
ordinary  linear  acoustic  waves.  In  the  discussion  that  follows,  the  spherically  expanding  blast  wave 
will  form  the  basis  for  such  a  model,  notwithstanding  the  fact  that  the  laser-induced  explosions 
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which  propel  the  vehicle  maybe  linear  rather  than  iphericaL  This  approach  is  justified  because  of 
the  fact  that  the  explosions  are  of  finite  length  and  occur  in  close  proximity  to  a  finite  baffle  (the 
lower  vehicle  surface).  The  results  of  a  more  detailed  analysis  incorporating  both  a  cylindrical  near- 
field  blast  wave  and  a  spherical  far-field  blast  wave  would,  to  some  extent,  modify  the  results  of  the 
simplified  analysis  performed  below,  but  would  not  materially  affect  the  conclusions  drawn  from  it. 

The  instantaneous  spatial  pressure  distribution  is  shown  qualitatively  in  Fig.  VIII-1  for  four 
increasing  values  of  elapsed  time  following  an  explosion  (Kinney  Si  Graham,  1985)[ll).  Analysis  of 
the  blast  wave  responsible  for  such  a  pressure  distribution  was  performed  originally  in  order  to 
predict  the  characteristics  of  nuclear  exploeions  (Taylor,  13&0)[27|.  The  equivalent  analysis  of  a 
cylindrical  wave  was  performed  by  Lin  (1954)[l2j.  Both  analyses  introduced  the  assumptions  of 
self-similar  behavior  inherent  in  the  scaling  laws  invoked  by  Reilly  et  aL  (1979)^21]  in  their  study  of 
momentum  transfer  arising  from  laser-induced  explosions  near  non-ablating  surfaces. 

Some  results  of  Taylor’s  analysis  which  are  particularly  appropriate  here  are:  (a)  Blast  wave 
profiles  consist  of  a  strong  shock  front  behind  which  pressure,  temperature,  and  density  rapidly 
decrease;  (b)  Pressure  profiles  at  any  instant  in  time  are  self-similar  -  that  is,  similar  if  the  pressure 
is  normalized  with  respect  to  the  peak  proeure  immediately  behind  the  shock,  and  distance  is 
normalized  with  respect  to  the  shock  radius;  (c)  Peak  pressure  behind  the  shock  front  declines  with 
the  inverse  cube  of  the  shock  radius;  and  (d)  Speed  of  the  shock  front  declines  as  the  inverse  3/2 
power  of  the  shock  radius.  Analogous  results  hold  (except  that  the  respective  powers  of  the  shock 
radius  are  different)  in  the  case  of  cylindrical  shocks. 

The  similarity  rules  governing  blast  waves  cannot  hold  far  from  an  explosion  because:  (a) 
Small-amplitude  (acoustic)  waves  decline  inversely  with  the  first  power  of  the  radius  of  the  wave 
front  (the  equivalent  of  the  shock  front  in  the  case  of  blast  waves);  and  (b)  The  wave  profile  does 
not  change  shape  with  time  or  distance  -  i.e.  the  wave  is  non-  dispersive.  Attention  must  therefore 
be  paid  to  those  circumstances  which  prevail  while  an  explosion-induced  disturbance  undergoes  the 
transition  from  blast  wave  to  sound  wave. 

C.  BLAST  WAVE  APPROXIMATION 

Some  insight  can  be  obtained  from  an  examination  of  Fig.  VIII-2,  in  which  the  solid  line  is  a 
plot  of  data  obtained  from  page  32  of  Kinney  Si  Graham  (1985)[ll|.  The  ratio  of  blast  overpressure 
AP  to  the  ambient  reference  pressure  Po  is  represented  on  the  abscissa,  while  the  ratio  of  shock  radius 
R  to  a  reference  radius  Ro  is  represented  on  the  ordinate  (the  reference  radius  Ro  has  arbitrarily 
been  defined  as  the  value  of  R  at  which  A P  equals  lOOPo).  The  point  corresponding  to  the  actual 
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Figure  Vm-1  Typical  Preeeure  for  Sacceeetve  Time*  After  fixploelon 
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location  and  magnituds  of  peak  overpressure  at  a  given  value  of  time  after  the  exptocion  will  move 
along  the  aolld  line  from  right  to  left  aa  time  increases.  Because  of  this  behavior,  which  is  indicated  in 
Fig.  Vm-2  by  the  arrow  labeled  DIRECTION  OF  WAVE  PROPAGATION,  a  plot  of  peak  pressure 
versus  location  will  hereafter  be  referred  to  as  a  Prttttre-Dirttnce  Trtjeeitrt  (PDT). 

The  dashed  line  tangent  to  the  left  end  of  the  solid  line  has  the  slope  corresponding  to  the 
l/R  pressure  variation  of  linear  acoustic  waves,  while  the  dashed  line  tangent  to  the  right  end  of  the 
solid  line  has  the  slope  corresponding  to  the  1/J2®  pressure  variation  of  strong  blast  waves.  Therefore 
the  right-hand  extremity  of  the  solid  line  lies  in  the  blast  wave  regime  while  the  left-hand  extremity 
lies  in  the  acoustic  regime.  Inspection  of  Fig.  VHI-2  reveals  that,  when  viewed  on  a  logarithmic 
rather  than  a  linear  scale,  the  transition  of  an  explosion-induced  disturbance  from  a  nonlinear  blast 
wave  to  a  linear  acoustic  wave  is  centered  in  the  neighborhood  of  the  point  B,  which  is  the  closest 
point  on  the  solid  line  to  the  point  A  where  the  two  tangent  lines  cross. 

This  suggests  that  the  actual  PDT  can  be  approximated  by  the  left-hand  tangent  line  for 
linear  acoustic  waves,  and  by  the  right-hand  tangent  line  for  blut  waves,  with  transition  occurring 
at  the  single  point  A.  Analysis  would  begin  with  formulation  of  blast  wave  conditions  near  the  source 
explosion,  and  the  blast  wave  scaling  laws  would  then  be  used  to  determine  the  premure  signature 
corresponding  to  point  A.  This  would  then  be  used  to  establish  the  acoustic  farfield  signature  for 
large  values  of  R.  An  estimate  of  the  error  likely  to  be  incurred  by  utilizing  the  approximate  PDT  can 
be  made  by  noting  that  the  ratio  of  the  overpressure  at  the  actual  midpoint  B  of  transition  in  Fig. 
Vm-2  to  that  at  point  A  is  about  4/3  corresponding  to  an  under-  estimation  of  20  Iog(4/3) «  2, 5 dB. 

D.  BLAST  WAVE  SPECTRUM 

Inspection  of  Fig.  VIII-1  reveals  that  a  single  blast  wave  profile  displays,  at  a  given  instant 
in  time,  an  approximately  exponential  variation  of  pressure  with  distance  from  the  blast  center, 
at  least  in  the  neighborhood  of  the  shock  front.  This  suggests  that,  to  a  first  approximation,  the 
pressure-time  signature  observed  at  a  stationary  observation  point  ai  the  blast  wave  passes  by  it 
will  exhibit  a  sharp  rise  (due  to  the  passage  of  the  shock  front)  followed  by  an  initially  (at  least) 
exponential  pressure  drop  with  time  due  to  the  passage  of  the  (spatially)  exponentially  varying  self- 
similar  pressure  profile  behind  the  shock  front.  The  validity  of  this  approximation  was  proved  by 
Taylor  (1950)(27|  for  spherical  blast  waves. 

The  ERH  thruster  will  be  powered  not  by  a  single  laser-  induced  explosion,  however,  but 
by  a  series  of  explosions  occurring  at  a  prescribed  Pulse  Repetition  Frequency  (PRF).  Therefore  the 
resulting  pressure  signature  will  appear,  as  shown  in  Fig.  VIII-3,  in  the  form  of  a  steady-periodic 
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Figure  Vm-2  Plot  of  Peak  Blast  Wave  Pressure  vs.  Distance  from  the  Source. 
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oscillation  each  cycle  of  which  consists  of  a  atop- function  rise  followed  by  an  exponential  decay,  the 
period  T  of  the  cycle  beiaf  equal  to  l/PRF. 

The  prassure-dme  signature  p(t),  bring  periodic,  can  be  represented  u  a  Fourier  series: 

P(0*  s--  -jnu*t)  (VIII-  2) 


where  t  is  the  time,  j  *  y^T,  uq  =  2»PRF,  and 


The  pressure-time  signature  during  one  cycle  (t  =  -r/wo  to  r/wo)  is  given  by 


(VIII -3) 


p{t)*piel-*<‘+*)) 


(V//J-4) 


where  pi  is  the  maximum  pressure  at  the  shock  front.  If  the  indicated  integration  is  carried  out  we 
obtain,  after  some  rearrangement, 

*0 = „[i  -  *.(«-.< +*.)!  [viu- 5) 

where  =  tan'1  JL.. 

In  order  to  obtain  spectral  components  of  p(t)  the  constants  p  and  b  must  be  evaluated.  In 
carrying  out  the  evaluation,  it  is  important  to  note  that  p(t)  in  Eqn.  VQI-4  is  the  total  absolute 
pressure,  not  the  overpressure  AP  of  Fig.  VUI-2.  In  terms  of  the  notation  of  that  figure,  pi  = 
A P  +  P0.  The  ambient  reference  pressure  Po  is  (at  least  in  the  acoustic  far-fleld)  the  time-average 
value  of  p(t)  in  Eqn.  VIU-6  -it. 

Po  =  p,Il-e<^]^  (VIII  -  6) 

At  point  A  in  Fig.  VUI-2,  AP/Po  is  about  0.3,  so  that  pi  =  Po  +  AP  =  1.3Po.  If  we  substitute 
this  relation  for  pi  into  Eqn.  VIU-6  while  defining  the  variable  X  =  in  order  to  simplify  the 
resulting  expression,  we  get  the  following  transcendental  equation: 


Po  =  1.3P0(1  - 


(VIII  -  7) 


This  equation  will  be  satisfied  by  X  »  0.65,  and  therefore  b  will  have  a  range  of  values  for 
the  smaller  vehicle  (900- kg,  2.5  m  diameter)  from  8.85  t~l  under  hovering  conditions  (PRF  *  16.1 
Ha)  to  292  s'1  under  maximum  (33  g)  acceleration  (PRF  =  531  Ha).  For  the  larger  vehicle  (5555- 
kg,  5  m  diameter)  the  corresponding  range  of  values  is  20.2  a-1  under  hovering  conditions  (PRF 


151 


EXPONENTIAL 

DECAY 


=  36 J  Hi)  to  668  s-1  under  maximum  (33  g)  acceleration  (1214  Hi).  The  variable  X  can  also  be 
substituted  into  Eqn.  VHI-5,  thereby  simplifying  the  form: 


p(<)  »  +  fj(-l)»(l  +  tininuot  +  *>)|  (V///-  8) 

where  =  tin-1 

The  RMS  amplitudes  p+  of  the  Fourier  components  of  the  pressure  signature  will  therefore 
be 

=  ^Jj|l  -  exp(-X)|(l+ [VIII -9) 

The  sound  pressure  levels  (SPL%)  due  to  these  amplitudes  it  point  A  in  Fig.  VH1-2  (which  esn  be 
defined  as 

SPL,  *  20  log(P, /Pf./)  [VIII  -  10) 

where  P,,/  is  the  reference  sound  pressure,  20  Pa)  will  have  the  form 

SPL,  =  20  log{-^L-{l  -  e— 1(1  +  ^i)"1/*}  {VIII  -  11) 

as  a  consequence  of  this  result. 

In  the  present  problem,  is  1  atmosphere  (about  101  kPa)  so  that  pi  =  1.3P»  *  131.3fcPa. 
If  this  value,  together  with  the  previously  determined  value  of  0.55  for  X,  is  substituted  into  Eqn. 
Vm-11  the  values  of  the  SPL •  can  be  calculated.  The  following  Table  lists  the  values  for  the  first 
10  spectral  components: 

Table  VDI-J — SPL  Values  for  the  First  Ten  Spectral  Components 


n 

SPL,[db) 

n 

SPL,[db) 

1 

170.6 

6 

155.2 

2 

174.7 

7 

153.8 

3 

161.2 

8 

152.7 

4 

158.7 

9 

151.6 

5 

156.7 

10 

150.7 

For  values  of  n  greater  than  10,  SPL ■  can  be  obtained  from  the  approximate  equation 

SPL,  »  170.7  -  20  logn  [VIII-  12) 

Note  that  the  magnitudes  of  the  spectral  components  are  independent  of  the  PRF.  Note  also  that  if 
an  error  of  0.1  dB  is  acceptable,  then  Eqn.  VIH-12  can  be  used  to  predict  the  values  of  the  spectral 
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components  for  n  in  the  rug*  from  1  to  10.  Aa  free-fleld  measurement*  rarely  are  observed  to 
exhibit  repeatability  within  s  tolerance  of  0.1  dfi,  it  is  reasonable  to  regard  Eqn.  V HI-12  a*  bring 
an  adequate  approximation  (ie.  within  the  limits  of  experimental  error)  of  the  spectral  componta 
at  a  distance  corresponding  to  point  A  in  Fig.  VHI-2. 

E.  FAB-FIELD  SOUND  SPECTRUM  LEVEL 

Estimation  of  far-field  sound  levels  of  the  spectral  components  requires  that  the  distance 
R  from  the  blast  center  corresponding  to  point  A  in  Fig.  VIII-2  be  determined.  For  the  6666- 
kg  vehicle  at  sea  level  conditions  the  maximum  expansion  of  a  typical  plasma  column  is  about  0.4 
meters  in  diameter  and  2.6m  in  length,  enclosing  a  volume  of  about  0.16m3.  This  volume  is  achieved 
assuming  an  expansion  to  ambient  pressure  of  plssma  in  a  half  cylindrical  column  with  maximum 
plasma  pressure  varying  inversely  with  the  square  of  the  column  diameter.  This  same  expansion, 
when  it  has  occurred  only  to  the  extent  when  the  preasun  is  130  percent  of  ambient  (the  condition 
at  point  A  in  Fig.  VIII-2)  will  occupy  a  volume  of  only  0.12m*. 

If  we  now  assume  that,  for  acoustical  purposes  only,  the  plasma  column  is  equivalent  to  a 
spherical  blast  of  such  a  magnitude  that,  when  it  has  expanded  to  reach  point  A  in  Fig.  VIII-2  the 
volume  enclosed  by  the  shock  front  equals  the  pnviously  calculated  value  of  0.12m*  ,  then  a  formula 
for  the  fap-fleld  levels  of  the  spectral  components  can  be  readily  obtained.  The  radius  of  a  0.12m* 
sphen  is  about  0.3  m  so  that,  because  sound  pressure  levels  vary  as  ten  times  the  logarithm  of  the 
squan  of  the  ratio  of  the  distance  R  from  the  blast  center  to  the  sphere  radius  Ao,  a  factor 

10  log(A/0.3)a  =  10.3  +  20  log  A  ( VIII  -  13) 

must  be  subtracted  from  the  components  in  Eqn.  VIII- 12.  Therefore  the  far-fleld  levels  of  the 
spectral  components  will  be  given  by 

SPL,  =  170.7  -  20  logn  -  (10.3  +  20  log  A) 

5PL,  =  160.4 -20  logs -20  log  A  {VIII-  14) 

where  r  is  in  meters. 

This  result  holds  for  the  5666-kg  vehicle,  but  can  also  be  applied  to  the  900-kg  vehicle  if  2 
dB  is  subtracted.  The  need  for  the  correction  arises  from  the  fact  that  in  the  case  of  the  smaller 
vehicle  the  radius  of  the  cylindrical  shock  corresponding  to  a  30%  overpressure  is  approximately  the 
same  as  that  in  the  larger  vehicle,  but  the  shock  length,  and  therefore  the  shock  volume  used  to 
calculate  radius  of  the  equivalent  spherical  blast,  is  half  the  previous  value.  This  means  that  Ao  is 
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reduced  by  the  factor  2*,  twenty  times  the  logarithm  of  which  is  2  dB.  In  addition,  the  proximity 
of  the  vehicle  surface  enhances  the  sound  pressure  levels  achieved  for  the  same  reason  that  a  baffle 
lnhances  the  sound  level  of  a  loudspeaker  driver.  The  enhancement  is  directional,  with  the  greatest 
intensity  being  ichived  underneath  the  vehicle. 

A  conservative  estimate  of  the  enhancement  can  be  obtained  by  »— that  the  vehicle 
surface  functions  as  an  infinite  plane  baffle,  with  a  resulting  enhancement  of  6  dB.  Adding  this 
amount  to  the  spectral  components  obtained  above  yields 

SPL,  =  166.4  -  20  log  n  -  201ogR  (V///-  15) 

for  the  spectral  components  produced  by  the  5566- kg  vehicle.  The  6  dB  enhancement  is  partially 
negated,  in  the  case  of  the  900-kg  vehicle,  by  the  2  dB  scale  reduction  effect  noted  above,  so  that 
for  this  vehicle 

SPL,  s*  160.4  -  20  logn  -  20  logfl  ( VIII  - 16) 

At  large  distances  the  value  of  R  will  be  essentialy  equal  to  r,  the  distance  between  the  vehicle  and 
the  point  at  which  the  noise  level  is  evaluated.  This  fact  will  be  utilized  in  the  following  discussion 
of  overall  noise  level 

P.  OVERALL  FARrPIELD  SOUND  LEVEL 

The  overall  sound  pressure  level  (SPL)  is  obtained  by  the  formula 

oo 

SPL  =  10  log  £  io sPi,lio  {VIII  -  17) 

*5=1 

where  SPL •  is  as  given  in  Eqn.  VHI-16  and  VIO-16.  After  substitution  of  these  equations  into 
VIII-17  followed  by  some  rearrangement,  one  obtains 

SPL  =  166.4+  10  log  £  ^  -20  log  R  (VIII -IS) 

e 

for  the  6666-  kg  vehicle  and 

SPL  =  160.4  +  10  log£  ^  -  20 log/?  ( VIII -  19) 

» 

for  the  900-kg  vehicle.  The  value  of  the  infinite  sum  is  so  that 

101ogI}^  =  10lo*T=2-2  (VIII- 20) 
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and  this,  in  turn  yields 


(V  III  -  21) 


SPL  *  168.6  -  20  log  12 

for  the  5566-kg  vehicle  and 

SPL  =  165.6  -  20  log  J2 


(VIII -22) 


for  the  900-kg  vehicle. 

These  last  equations  permit  the  construction  of  a  somewhat  more  reliable  table  of  SPL 
versus  distance  r  than  Table  Vm-1  (q.v.): 


Table  VHI-3 — SPL  vs.  Distance  From  Vehicle 


Distance  r 

spl  (5665-kg  vehicle) 

spl  (900-kg  vehicle) 

20  m 

142.6 

140.6 

50  m 

134.6 

132.6 

100  m 

128.6 

126.6 

200m 

122.6 

120.6 

500m 

114.6 

112.6 

1  km 

108.6 

106.6 

2  km 

102.6 

100.6 

5  km 

94.6 

92.6 

10  km 

88.6 

86.6 

20  km 

82.6 

80.6 

If  this  set  of  data  is  compared  with  that  In  Table  Vm-1  It  is  seen  that  a  more  detailed  analysis 
results  in  a  5.6  dB  increase  in  predicted  overall  noise  level  above  that  given  by  Eqn.  VIII-3  for  the 
900-kg  vehicle,  and  a  2.2  dB  decrease  in  predicted  overall  noise  level  below  that  given  by  Eqn.  VED-2 
for  the  6566-kg  vehicle.  Note  that  these  results  hold  only  if  laser-induced  explosions  occur  only  one 
at  a  time.  If  m  >  1  explosions  occur  at  one  time,  distributed  over  the  surface  of  the  vehicle,  then 
the  quantity  20  logm  must  be  added  to  the  data  in  Table  VIII-3.  Therefore  minimum  noise  is  likely 
to  be  achieved  when  only  one  explosion  at  a  time  takes  place. 

G.  SUMMARY 

The  preceding  sections  show  that  the  employment  of  a  crude,  simple-source  type  acoustical 
model  of  an  ERH  thruster  propelled  vehicle  in  hover  yields  approximately  the  same  results  (to  within 
a  few  decibels)  as  a  more  detailed  model  employing  concepts  from  the  theory  of  blast  waves  in  the 
derivation  of  a  sound  spectrum.  While  the  assumptions  invoked  in  formulating  the  latter  model  are 
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admittedly  conservative,  it  k  unlikely  that  results  obtained  from  the  analysis  of  a  more  elaborate, 
lean  conservative  model  would  differ  from  thoee  obtained  here  except  in  detail. 

In  particular,  the  above  results  suggest  that  the  overall  noise  level  produced  by  ERH  thruster 
propelled  vehicles  under  uniform,  non- dissipative  atmospheric  conditions  in  the  absence  of  ground 
proximity  effects  can  be  considerable.  The  mitigating  effects  (if  any)  of  atmospheric  and  ground 
absorption  are  too  complicated  and  uncertain  to  be  included  within  the  scope  of  the  present  work, 
but  should  be  considered  in  future  efforts. 

In  one  respect  the  blast  wave  model  does  produce  a  substantial  reduction  in  predicted  noise 
level  in  comparison  with  the  earlier  simple-source  model  It  wu  observed  previously  that  this  model 
predicted  an  increase  in  overall  noise  with  vehicle  acceleration,  the  augmentation  reaching  30.4 
dB  under  33  g  acceleration.  The  more  accurate  blast  wave  model,  on  the  other  hand,  predicts  no 
augmentation  whatsoever  with  acceleration.  This  is  because  thrust  increases  required  for  acceleration 
an  achieved  by  increasing  the  PRF  rather  than  the  blast  intensity.  The  model  also  suggests  that 
one  way  to  minimize  the  overall  noise  level  Is  to  avoid  having  mon  than  one  laser-induced  explosion 
occur  at  a  time  (Len  the  rotary  detonation  wave  engine  approach). 

The  effect  of  vehicle  acceleration  on  perceived  loudness  (the  subjective  level  of  noise)  is 
another  matter,  however,  because  this  quantity  is  frequency  dependent.  In  general  we  can  expect 
Increasing  PRF  to  achieve  acceleration  to  increase  the  perceived  noise  level  because  of  the  increase  of 
sensitivity  with  frequency  of  the  the  human  ear  in  the  range  0-2  kHz.  The  actual  increase  in  perceived 
noise  level  realized  at  large  distances  from  the  vehicle,  where  frequency-dependent  atmospheric  and 
ground  absorption  effects  cannot  be  ignored,  are  beyond  the  scope  of  the  present  work  and  should 
be  considered  later. 

We  therefore  recommend  that  the  work  begun  here  be  continued  in  pursuit  of  two  objec¬ 
tives:  (a)  Assessment  of  the  noise  spectrum  far  from  the  vehicle  including  frequency-dependent 
atmospheric  and  ground  absorption  effects,  (b)  Assessment  of  perceived  noise  levels  far  from  the 
vehicle  under  conditions  ranging  from  hover  to  maximum  acceleration.  As  an  aid  to  verifying  these 
results  we  recommend  that  a  digitally  synthesized  simulation  of  the  anticipated  noise  signature  be 
made.  Such  a  simulation  can  be  accomplished  using  equipment  in  the  Laboratory  for  Noise  and 
Vibration  Control  Research  (LNVCR.  Finally,  we  recommend  that  additional  work  be  carried  out 
to  develop  an  acoustical  model  for  predicting  the  noise  level  and  lift  enhancement  achieved  when 
operating  an  ERH  thruster  vehicle  very  close  to  the  ground. 

Documented  in  the  following  Chapter  is  a  n  account  of  progress  made,  to  date,  and  the 
digitally  synthesized  simulation  of  these  anticipated  noise  signatures. 
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CHAPTER  DC 


SIMULATION  OP  ERH  THRUSTER  NOISE  SIGNATURES 

A.  Introduction 

The  purpoee  of  the  effort  described  here  in  Chapt.  9  is  to  produce  acoustic  signatures  which 
approximate  those  that  would  be  perceived  by  a  stationary  ground-based  observer  during  various 
flight  maneuvers  of  the  reference  vehicle  (Apollo  Lightcraft).  Several  assumptions  are  made  in  order 
to  simplify  both  the  analysis  of  the  engine/vehicle  as  an  acoustic  source,  and  the  creation  of  the 
synthesized  signatures: 

1.  The  engine/vehicle  is  considered  to  be  a  point  acoustical  source.  Directivity  effects  are 
ignored. 

2.  Only  the  direct  acoustic  signature  is  modeled;  no  attempt  is  made  to  include  effects  due  to 
reverberation  (from  the  atmospheric  and  ground  surface). 

3.  The  speed  of  sound  in  air  is  assumed  to  be  constant  and  not  a  function  of  altitude. 

4.  No  attempt  is  made  to  include  sonic  booms  caused  by  the  vehicle  exceeding  Mach  1. 

This  project  essentially  combines  a  model  of  the  acoustic  response  of  the  Apollo  Lightcraft 
ERH  thrusters  to  a  laser  pulse  (i.e.  the  acoustic  impulse  response)  with  various  sets  of  vehicle 
trajectory  and  laser  pulring  (i.e.  forcing  function)  data.  The  acoustic  impulse  response  model 
is  presented  in  Chapter  9,  and  the  vehicle  trajectory /forcing  function  model  is  reviewed  below  in 
Section  G  of  this  chapter. 

B.  Model  Description 

The  computer  program  links  the  acoustic  impulse  response  model,  which  was  developed  for 
a  stationary  source/stationary  observer,  to  the  trajectory/forcing  function  data-in  order  to  generate 
acoustic  signatures  for  a  traveling  source/stationary  observer.  The  effect  of  Doppler  shifting  and 
amplitude  attenuation  of  the  observed  source  signature  are  incorporated  as  a  function  of  the  vehicle 
dynamics.  The  resulting  expressions  for  the  traveling  source/stationary  observer  signatures  form 
the  basis  of  a  computer  program  which  generates  a  sequence  of  numbers  that  represents  digitized 
samples  of  the  analog  signatures.  Next,  these  digitized  samples  are  processed  by  a  digital- to- analog 
converter  (DAC),  then  low  pass  filtered  to  remove  harmonics  of  the  sampling  frequency  produced  by 
the  DAC,  and  finally  recorded  on  magnetic  tape.  The  instrumentation  setup  is  shown  schematically 
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in  riff.  DC-1. 

The  overall  proem  U  regarded  u  a  cascade  of  two  processes,  u  illustrated  in  Fig.  IX-2. 
Hie  flnt  proceae  H\  describes  the  static  acoustic  impulse  response  K«(t)  of  the  ERH  thruster  to  a 
laser  pulse,  where  both  the  vehicle  and  the  observer  are  stationary.  The  acoustic  Impulse  response 
Hi  is  a  function  of  PRF  only.  Figure  IX-3(a)  qualitatively  shows  the  static  acoustic  response  of 
the  thruster  to  a  600  Hi  PRF  pulse  train,  while  Fig.  IX-3(b)  shows  the  response  to  an  1136  Hz 
PRF  pulse  train.  Process  H%  describes  the  time-varying  nature  of  the  transmission  path  from  the 
moving  thruster/vehicle  to  the  fixed  observer;  it  combines  the  output  K,(t)  from  the  static  model 
Hi  together  with  relative  distance  data  x(t)  and  relative  velocity  data  v(t)  to  yield  th*  dynamic 
acoustic  response  Ko(t)  for  a  moving  source  and  fixed  observer. 

C.  Model  Implementation 

The  decision  to  implement  the  synthesis  using  a  digital  system  was  based  largely  upon  three 
considerations: 

1.  The  acoustic  impulse  response  does  not  represent  one  of  the  standard  signal  types  available 
on  hard-wired  analog  function  generators; 

2.  The  acoustic  impulse  response  changes  with  frequency; 

3.  Precise  control  of  the  signature  amplitude  and  frequency  is  required  to  reflect  changes  in 
the  acoustic  transmission  path  due  to  vehicle  dynamics. 

The  natural  choice  for  such  arbitrary  waveform  generation  capability  is  a  digital  computer 
based  system.  This  approach  does,  however,  introduce  several  additional  considerations  which  must 
be  addressed  in  order  to  achieve  the  desired  result.  A  brief  discussion  of  these  considerations  is 
presented  here.  The  goal  is  to  provide  motivation  for,  as  well  as  describe  In  detail,  the  approach 
used  to  digitally  synthesize  the  acoustic  signature. 

Figure  DC-3  depicts  the  nature  of  the  theoretical  static  acoustic  impulse  response  F«(t)  of 
the  Apollo  Lightcraft.  The  figure  was  generated  by  using  a  digital  computer  to  calculate  K«(t)  for  a 
number  of  different  values  of  time  t.  Enough  values  were  calculated  so  that  when  the  corresponding 
points  were  plotted  and  connected  with  straight  line  segments,  the  resulting  figure  closely  approxi¬ 
mates  a  graphical  representation  of  the  actual  analog  signal.  This  is  essentially  the  approach  used 
to  generate  an  analog  signal  from  a  set  of  discrete  data  where  a  DAC  is  used  as  the  output  device  in¬ 
stead  of  a  plotter.  The  main  problem  is  to  determine  how  many  data  points  are  required  to  produce 
a  satisfactory  approximation  to  the  actual  analog  acoustic  response  signature  Y,(t). 

If  the  discrete  data  is  sent  to  a  DAC  instead  of  a  plotter,  an  analog  signal  will  result.  There 
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Figure  DC-3  (a).  Theoretical  Acoustic  Impulse  Response-  PRF  *  600  Ri 


Figure  IX-3(b).  Theoretical  Acoustic  Impulse  Respouse-PRF  =  1135 


•n  a,  number  of  different  varieties  of  DAC  available,  but  the  moat  common  is  the  taro  order  hold 
voltage  DAC.  The  resolution  of  the  DAC  ia  defined  to  be  the  magnitude  of  the  imaiieet  change  in 
output  voltage  that  can  be  attained  under  program  control,  divided  by  the  magnitude  of  the  largest 
change  in  output  voltage  that  can  be  attained  under  program  control,  expressed  in  per  cent.  The 
69720A  DAC  is  a  12-bit,  zero  order  hold  voltage  DAC  with  a  resolution  of  (1/21*)  x  100  s  0.024%. 
The  term  sero  order  hold  refers  to  the  fact  that  if  a  data  value  representing  a  voltage  Vj  is  written 
to  the  DAC,  the  output  of  the  DAC  will  change  (almost)  Instantaneously  to  V\  and  will  remain 
constant  until  a  new  data  value  is  written. 

As  one  might  expect,  the  approximation  of  a  given  analog  signal  by  a  digitally  synthesised 
signal  generally  improved  as  the  frequency  Qf  at  which  the  analog  signal  is  sampled  increases.  Two 
factors  tend  to  limit  the  choice  of  O*  and  hence  limit  the  accuracy  of  the  approximation: 

1.  Each  digital  to  analog  conversion  requires  a  finite  length  of  time  to  occur. 

2.  Finite  machine  memory  limits  the  number  of  data  values  which  can  be  stored  prior  to  bring 

sent  to  the  DAC. 

The  69720A  DAC  can  process  data  at  a  rate  of  up  to  30  kHz  real  time  during  disk  to  DAC  transfers. 
The  236C  computer  system  has  approximately  15  Mbyte  of  free  disk  space  to  hold  data  for  the 
DAC.  We  describe  next  how  these  constraints  Impose  restrictions  on  the  class  of  signals  that  can  be 
satisfactorily  synthesized. 

D.  Digitltal  Simulation  Constraints 

Refer  to  Fig.  XI-4(a)  where  several  cycles  of  the  function  f\(t)  -  sin0.8rt  are  plotted.  This 
function  is  periodic  with  a  period  of  2r/0.8r  =  2.5  seconds.  Suppose  fi{t)  is  sampled  at  one  second 
intervals,  which  is  equivalent  to  saying  that  the  sampling  frequency  w,  is  two  and  one  half  times 
the  frequency  of  ft  (t)  evaluated  at  t  =  n,  n  *  0, 1,2,  —  These  are  the  values  of  f(t)  indicated  by 
the  tic  marks  in  Fig.  IX-4(a).  If  this  sequence  is  subsequently  sent  to  a  zero  order  hold  DAC,  the 
resulting  analog  waveform  will  appear  as  in  Fig.  IX-4(b). 

Next  consider  Fig.  IX-4(c),  where  several  cycles  of  /a(t)  =  sin  2.8*1  are  shown.  If  /»(t) 
is  sampled  at  one  second  intervals  and  the  resulting  sequence  is  sent  to  the  DAC,  the  output  is  as 
shown  in  Fig.  DC-4(d).  Notice  that  the  output  from  the  DAC  is  the  same  for  both  f\(t)  and  h(t). 
This  phenomenon  is  called  aliasing  of  and  occurs  whenever  a  periodic  function  is  sampled  at 
less  than  (or  precisely  at)  twice  its  frequency.  The  signal  in  Fig.  IX-4(b)  may  be  a  rather  crude 
approximation  to  the  original  signal  in  Fig.  IX-4(a),  but  at  least  it  tends  to  track  fi(t)  as  it  oscillates. 
This  is  clearly  not  the  case  for  Figs.  IX- 4(c)  tad  iX-4(d).  Ceashdmthsa  el  this  pbsaMsma  pswriiwi 
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Flgnn  DC-4  (a).  Aliasing  and  Sampling-Induced  Modulatioii-(analog  signal) 


Figure  DC-4(b).  Aliasing  and  Sampling-  Induced  Mod  illation- (sampling  /t(<)  at 
a  rate  of  1  sample/eec) 
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motivation  for  iho  following  criterion: 

Ntiuitt  Criterion.  In  order  to  uniquely  recover  a  band-limited  signal  from  its  sample  sequence,  the 
signal  must  be  sampled  at  a  rate  of  greater  than  twice  the  frequency  of  any  constituent  components. 

Notice  that  if  the  original  function  is  not  band-  limited,  it  is,  strictly  speaking,  impossible  to  recover 
it  exactly  from  its  sample  sequence.  Because  the  69T20A  DAC  has  a  maximum  sampling  rate  of 
30,000  samples  per  second,  the  Nyquist  criterion  demands  that  the  input  data  to  the  DAC  represent 
samples  obtained  from  an  original  signature  having  a  bandwidth  of  less  than  15  kHz.  This  represents 
one  important  constraint  on  the  class  of  signatures  that  can  be  accurately  synthesized.  The  second 
constraint,  to  be  discussed  next,  is  even  more  severe. 

Refer  once  again  to  Fig.  IX-4(b).  Although  f  (t)  is  sampled  in  accordance  with  the  Nyquist 
Criterion,  the  signal  produced  by  the  DAC  exhibits  an  amplitude  modulation  which  is  not  present 
in  the  original  signal.  For  the  purposes  of  synthesizing  signatures  from  sampled  data,  this  modu¬ 
lation  represents  an  unwanted  artifact.  Sampling  of  a  cosinusoid  of  frequency  wj,  using  a  sampling 
frequency  such  as  w,  =  (r/s)<i/i,  where  the  numbers  r  and  s  are  relatively  prime,  results  in  a  periodic 
sample  sequence  of  period  TA.u.  =  s/wi  *  r/w„  For  instance.  Fig.  IX-4(b)  shows  the  result  of 
sampling  /i(t),  which  has  frequency  u>i  *  0.8s /sec,  at  a  sampling  frequency  of  w,  ■  2s /sec.  Then 


MC 
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Hence  the  resulting  amplitude  modulated  sample  sequence  should  have  a  period 

Tam.  =  2Ti  *  -I  =  A  =  5T,  (/*-4) 

wi  u, 

and  this  is  verified  In  Fig.  IX-4(b). 

It  can  be  shown  that  the  frequency  spectrum  of  a  carrier  signal  of  frequency  w«  modulated 
by  a  waveform  of  frequency  um  will  contain  a  component  at  u/„  a  component  at  u,  +  um,  and  a 
component  at  -  um.  The  latter  two  components  are  referred  to  as  sidebands.  The  presence 
of  sidebands  is  a  consequence  of  modulation  caused  by  sampling  and  elimination  of  the  sidebands 
requires  removal  of  the  modulation. 
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Figure  DC-4(c)  Alto— dag  end  Sampling- ( analog  signal  mn(2.8rt)) 


1-0 


Figure  £X-4(d)  Alia— sing  and  Sampling* Induced  Modulation-!— mpling  ft{t) 
a  rate  of  1  sample /sec) 


E.  Redaction  at  Usdidnd  Modulation  Effects 


Modulation  of  the  original  signal  by  the  sampling  function  signal  will  always  occur  to  some 
extent  when  the  sampling  frequency  is  a  non-integer  multiple  of  the  fundamental  frequency  of  the 
original  periodic  signal  The  only  recourse  is  to  reduce  its  effect  to  an  acceptable  level,  if  possible. 
Two  means  of  accomplishing  this  are: 

1.  If  the  signal  is  band-limited,  increasing  the  sampling  frequency  will  reduce  the  amplitude 
of  the  modulation,  and  hence  the  amplitude  of  the  sidebands.  It  can  be  shown  that  the 
amplitude  of  modulation  of  a  cosinusoid  f(t)  =  cos  wt  produced  by  sampling  at  a  frequency 
uit  =  oru/i,  a  >  2,  must  satisfy 

A.M.  <  cos(0)  -  cos(^)  (IX  -  5) 


or 

A.M.<l-cos(^)  (/X-6) 

For  example,  if  /(f)  =  cotuit  is  sampled  at  a  rate  of  u,  =  10wj,  the  maximum  poeaible 
amplitude  modulation  due  to  sampling  is 


Maximum  A.M.  =  1  -  cos(^) 
=  0.191 


=  19.1%  modulation 


(IX -7) 


=  -14 dB 


Given  no  other  constraints,  it  would  seem  that  a  sampling  rate  should  be  chosen  to  ensure 
a  sampling-induced  modulation  of  less  than  -120  dB  for  any  frequency  component  in  the 
acoustic  signature,  since  this  represents  the  dynamic  range  of  the  typical  human  ear.  Not 
only  would  this  be  impractical  given  the  limitations  on  the  sampling  frequency  of  the  69720A 
DAC,  but  it  is  unnecessary  because  the  dynamic  range  of  a  12  bit  DAC  such  as  the  69720A 
is  only  72  dB.  The  ratio  of  sampling  frequency  required  to  ensure  a  modulation  of  amplitude 
of  less  than  -72  dB  is  foowd  as  fclnws: 


-72dfl  =  201og(l-coa(^.)| 


cos(— )»0.9W749 

i.um 

a 

a  a  28 

Thb  results  motivates  the  adoption  of  a  second  sampling  criterion. 

Sampling  Criterion  2.  In  order  to  suppress  sampling-induced  amplitude  modulation  of  the  synthe¬ 
sized  signature,  and  hence  to  suppress  the  associated  sidebands,  the  analog  acoustic  impulse  response 
signature  must  be  sampled  at  a  rate  of  at  least  28  times  w/,  where  w/  is  the  frequency  of  the  highest 
frequency  component  in  the  signature. 

Note  that  this  b  a  conservative  criteria;  thb  criteria  may  be  relaxed  if  additional  information 
regarding  the  frequency  content  of  the  analog  signature  b  available. 

2.  A  second  means  of  dealing  with  the  objectionable  sidebands  b  amply  to  mask  them  by 
adding  white  noise  at  a  level  which  b  sufficient  to  interfere  with  the  sidebands,  but  not  to 
interfere  appreciably  with  the  carrier  tone.  Thb  b  known  as  dithering. 

A  number  of  preliminary  signatures  have  been  synthesized  where  the  criteria  described 
above  were  violated.  In  all  but  one  case  (that  of  a  very  low  PRF),  the  difference  tones  are  dearly 
audible.  On  one  very  high  frequency  trace,  the  signature  b  aliased  from  inaudible  frequencies  down 
to  the  audible  range.  An  updated  version  of  the  software  b  being  generated  which  adheres  to  the 
criteria  that  have  been  discussed  here.  Thb  updated  software  will  utilized  a  truncated  Fourier  series 
approximation  to  the  ideal  static  acoustic  impulse  response  signature  to  assure  adherence  to  the 
criteria  outlined  above.  In  addition,  the  least  significant  bit  of  the  DAC  will  be  dithered. 

The  6M4/69720A  DAC  system  only  supports  continuous  digital  to  analog  conversions  at  a 
rate  of  30  kHz.  However,  it  b  possible  to  achieve  an  effective  sampling  rate  of  250  kHz.  In  order 
to  achieve  an  effective  sampling  rate  of  250  kHz,  the  data  b  sampled  by  the  DAC  at  a  rate  of  25 
kHz  and  recorded  on  the  Nagra  IV-SJ  tape  recorder  which  has  a  10:1  ratio  of  recording  speeds. 
When  the  data  b  recorded  at  1.5  inch/second,  the  effective  250  kHz  sampling  rate  b  achieved.  Thb 
approach  was  used  to  generate  several  of  the  preliminary  traces,  and  it  will  be  used  for  all  future 
work.  Note  that  the  high  sampling  rat.  achievable  with  thb  approach  enables  us  to  retsin  more 
terms  of  the  Fourier  approximation  to  the  ideal  acoustic  impulse  response  than  would  be  allowed  by 
the  sampling  criteria  for  a  30  kHz  sampling  rate.  Hence  a  closer  approximation  to  the  ideal  impulse 
response  b  obtained. 

Having  described  the  major  considerations  for  constructing  the  system  model  for  H\,  the 
static  acoustic  impulse  response  of  the  Apollo  Lightcraft,  the  process  of  modeling  ffj,  the  acoustic 
transmission  path,  b  outlined  next. 


167 


P.  Model  for  Acooatic  TVonwni—lon  Path  Bj 

The  moving  source/stationary  observer  formula  for  calculating  Doppler  shift  of  a  pure  tone 
frequency  source  /,  Is  given  by 

/.  =  |r73ir|/.  («-»> 

where: 


/,=actual  frequency  of  the  source  (Hi), 
v(t)=the  radial  component  of  the  velocity  of  the  source 
relative  to  a  fixed  observer  (ft/sec),  where  velocity 
toward  the  observer  is  considered  positive. 
c»the  speed  of  sound  in  air  (ft/sec) 
fT- the  Doppler  shifted  source  frequency  (Hi). 

Although  the  formula  above  is  written  for  a  source  which  produces  a  single  frequency  /,(t),  the 
formula  is  valid  for  complex  tone  sources  as  well  If  a  source  signature  is  written  as  a  Fourier  series, 
each  component  of  the  series  is  shifted  u  indicated  above. 

The  amplitude  of  the  observed  signature  is  calculated  using 

At  =  A$9  (/X  -  10) 

where: 

A#0=the  peak  amplitude  of  the  acoustic  pressure  observed 
at  a  distance  X  from  the  acoustic  source 
X«=the  distance  in  feet  between  the  source  and  the  observer 
at  which  A  is  measured 

X(t)=the  actual  radial  distance  in  feet  between  the 
source  and  the  observer 

It  should  be  noted  that  for  the  purposes  of  simulating  the  acoustic  signature,  the  sole  criterion  for 
selecting  a  reference  pressure  A«0Xo  is  to  utilise  the  maximum  dynamic  range  of  the  12  bit  69720 A 
DAC.  The  actual  sound  pressure  level  of  the  syntheisised  signature  is  a  function  of  the  gain  (volume) 
on  the  equipment  used  to  play  back  the  recorded  signatures. 

The  parameters  /«((),  v(t),  and  z(t)  used  in  these  formulas  describe  the  PRF  and  vehicle 
dynamics  and  are,  with  the  possible  exception  of  /«(()  in  certain  cases,  continuous  functions  of 
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kim*.  In  «rd«r  to  implement  computer  algorithms  to  calculate  the  Doppler  ehift  end  attenuation 
it  is  necessary  to  choose  discrete  times  at  which  to  evaluate  the  parameters  v(t),  and  z(t). 
We  define  the  parameter  sampling  period  T  to  be  the  time  Interval  in  seconds  between  successive 
evaluations  of  the  parameters  /•(<),  v(t),  and  s(t).  We  define  the  parameter  sampling  frequency  Op 
to  be  the  inverse  of  the  parameter  sampling  period. 

Experience  has  shown  that  the  value  of  Op  must  be  chosen  carefully  to  reduce  quantization 
effects  to  an  unacceptable  level.  A  number  of  signatures  were  processed  for  various  sets  of  vehicle 
flight  data  utilizing  a  parameter  update  frequency  of  30  parameter  samples  per  second.  In  no 
case  were  amplitude  quantization  effects  audibly  detectable;  however,  frequency  quantization  was 
clearly  audible  for  signatures  having  time  rates  of  change  of  frequency  of  greater  than  roughly 
800  Hz/sec.  Research  has  demonstrated  that,  over  a  broad  frequency  range,  an  average  person 
can  discriminate  pitch  changes  of  roughly  0.5%,  whereas  amplitude  discrimination  is  only  12%Q. 
This  implies  that  both  frequency  and  amplitude  quantization  can  be  reduced  to  inaudible  levels  by 
choosing  a  parameter  sampling  frequency  such  that 

0.006  <  -AM  Tt)  [IX  -  11) 

where 

/r(t)=the  Doppler  shifted  source  frequency  in  Hertz, 
evaluated  at  time  t  seconds 
Tp=the  parameter  sampling  period  in  seconds. 

This  will  ensure  that,  for  any  given  time  step  T  ,  the  observed  frequency  change  by  no  more  than 

0.6%. 

G)  Single  Relay  Satellite  Trajectory  Simulation 

Simulation  of  an  accurate  lift-off  trajectory  for  the  Apoilo  Lightcraft  is  an  important  first 
step  in  generating  the  acoustic  signature  representative  of  those  which  would  be  heard  by  a  station¬ 
ary,  ground-based  listener.  A  number  of  simplifying  assumptions  were  made  to  constrain  the  variety 
of  trajectories  possible  for  a  given  set  of  initial  conditions  and  vehicle  flight  characteristics.  These 
assumptions  are  given  below: 

1  The  Lightcraft  is  powered  directly  from  a  single  relay  satellite.  All  laser  power  is  relayed 
from  this  satellite. 
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2.  The  low«r  relay  satellite  ia  in  &  circular  equatorial  orbit  at  100  miles  (163.2  km)  altitude. 
This  fairly  low  altitude  was  considered  to  be  a  lower  bound  for  viable  relay  station  orbits. 

3.  The  Lightcraft  must  always  point  towards  the  relay  satellite,  and  thus  moves  up  the  beam 
emanating  from  it. 

4.  Hie  Lightcraft  lifts  off  with  a  specific  initial  zenith  angle,  which  is  determined  by  the  type 
of  trajectory  to  be  flown.  An  initial  zenith  angle  of  zero  degrees  means  a  vertical  lift-off;  an 
initial  zenith  angle  of  seventy  degrees  implies  a  nearly  horizontal  lift-off.  For  the  trajectory 
simulations,  initial  zenith  angles  of  0  and  45  degrees  were  used. 

5.  The  vehicle  files  either  with  a  constant  acceleration  (as  felt  by  the  pilot),  or  a  constant 
thrust.  In  the  constant  acceleration  cases  thrust  must  vary  due  to  the  changing  drag  forces 
on  the  vehicle  (i.e.,  with  flight  speed). 

6.  The  coefficient  of  drag,  Co,  for  the  Apollo  Lightcraft  was  calculated  as  a  function  of  flight 
mach  number  only.  The  drag  coefficient  data  used  in  the  simulation  is  displayed  in  Fig. 
IX-5. 

7.  The  liftoff  engine  is  an  ERH  thruster,  with  either  a  rotating  line  source,  or  twelve  simultane¬ 
ous  line  sources.  If  the  line  sources  are  simultaneous,  the  PRF  (Pulse  Repetition  Frequency) 
of  the  vehicle  Is  reduced  by  12  x. 

Although  a  large  number  of  trajectories  were  “flown”,  only  three  representative  cases  are 
reported  here:  1)  Horizontal,  constant  thrust  fly  by,  2)  Vertical  takeoff  at  3G’s  acceleration,  3)  45° 
(initial)  takeoff  at  3  G’s  acceleration.  A  cassette  tape  recording  of  these  three  noise  signatures  is 
included  with  this  final  report. 

1)  Horizontal,  Constant-thrust  Flyby 

The  constant  altitude,  constant  velocity  flyby  is  one  of  the  simplest  trajectory /forcing  func¬ 
tion  that  can  be  used  to  generate  acoustic  signatures  for  a  traveling  source/stationary  observer.  It 
was  the  first  case  to  be  simulated,  and  hence  served  to  help  “debug”  the  system. 

Figure  EX-6  portrays  the  geometry  for  this  case.  The  vehicle  is  assumed  to  fly  over  the 
observer  at  an  altitude  of  50  ft,  and  a  velocity  of  100  MPH-much  like  a  helicopter.  The  ERH 
thruster  develops  a  constant  level  of  thrust  at  a  pulse  repetition  frequency  (PRF)  of  41.67  Hz.  This 
is  somewhat  higher  than  that  required  for  “stationary  hover”  (i.e.,  37.8  HZ).  As  indicated  in  Fig 
EX-6,  the  thrust  vector  is  tilted  slightly  forward  (In  the  direction  of  flight)  to  counteract  aerodynamic 
drag  forces  (  ~  56701bs). 

The  trajectory  simulation  routine  computes  the  line-of-sight  distance  from  the  vehicle  to  the 
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Fight  Mach  Number 


Figure  EX-5  Dreg  Coefficient  n  Flight  Mach  Number 


Figure  DC-6  Geometry  for  Horiiontal  Fly-by  at  constant  thrust,  50  ft  altitude, 
100  MPH  and  41.67  Hi  PRF. 


observer,  sod  the  vehicle  velocity  component  radial  to  the  observer.  The  PRF  in  this  esse  it  held 
constant  at  41.67  Hl  The  simulation  begins  when  the  vehicle  is  just  over  1/4  mile  (1467  ft)  away 
from  the  observer. 

The  acoustic  signature  for  this  case  is  displayed  in  Fig.  DC-7.  As  indicated,  the  pulse 
frequency  is  doppier*hifted  from  48  to  36.75  Hz  upon  passage  of  the  vehicle.  Also  the  normalized 
real  amplitude  of  the  ERH  thruster  increases  to  a  maximum  at  the  point  of  closest  approach. 

This  “test*  case  assumed  no  change  in  vehicle  pointing  throughout  the  entire  20  second 
duration  of  the  simulation;  i.e.,  the  relay  satellite  would  have  to  be  at  a  very  high  altitude.  For 
the  following  two  cases,  all  the  afore-mentioned  input  assumptions  are  applied-such  that  the  vehicle 
must  chase  the  low  altitude  relay  satellite  into  orbit. 

2)  Vertical  liftoff  at  3  G’s  acceleration 

To  calculate  the  location  and  flight  direction  of  the  Apollo  Ligbtcraft,  it  is  necessary  to 
determine  four  parameters  which  exactly  describe  the  launch  trajectory.  These  four  parameters  are: 

(i)  vehicle  altitude  -  “h" 

(ii)  zenith  angle  -  V 

(iii)  satellite  angle  -  V 

(iv)  vehicle  angle  -  *5* 

(v)  vehicle  to  satellite  distance  -  *0" 

The  relationship  of  these  key  trajectory  model  parameters  is  displayed  in  Fig.  IX-8,  which 
also  indicates  the  location  of  vehicle  take-off.  Note  the  current  positions  of  the  vehicle  and  relay 
satellite  in  relation  to  the  point  of  liftoff. 

The  initial  conditions  must  be  specified,  of  course,  before  the  trajectory  can  be  analyzed. 
Figure  IX-9  portrays  the  initial  liftoff  geometry,  for  which  the  vehicle  altitude,  h,  is  set  to  zero.  For 
the  case  at  hand,  the  zenith  angle  <j>o  is  also  set  to  zero,  for  the  3  G  constant  acceleration  launch. 
The  observer  is  assumed  to  reside  at  the  launch  site. 

Calculation  of  the  pulse  repetition  frequency,  PRF,  depends  on  the  ERH  thruster  mode  and 
the  level  of  thrust  it  is  generating.  For  the  “simultaneous  twelve-line-source"  model  for  the  ERH 
thruster,  a  PRF  of  37.8  Hz  is  required  to  hover  the  Apollo  Lightcraft.  This  frequency  generates 
about  12^250  lbs.  of  thrust;  hence,  the  PRF  for  a  specific  thrust  output  is 

PRFsl=  37.81^1  {IX  -  12) 

where  T  is  the  desired  thrust  in  pounds. 
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VEHICLE  PULSING  FREQUENCY  t  =  ^O.Sec 


Figure  EX-9  Initial  Liftoff  Geometry 


For  »  “rotating  lint  source”  model,  the  PRF  ia  incr— a»d  by  ft  factor  of  twelve.  Each  line 
source  is  pulsed  in  sequence.  This  has  the  effect  of  multiplying  the  PRF  by  12  x: 

PRFkl  “«4[T^S51  [IX  ~  13) 

The  last  two  3  G  simulations  assume  the  use  of  this  rotating-iine-source  mode  for  ERH  thruster 
operation. 

Figure  IX-10  gives  the  trajectory/forcing  function  data  for  the  constant  3  G  launch  from  an 
initial  zenith  angle  of  zero  degrees.  Figure  IX-  10a  gives  the  altitude  vs.  downrange  distance;  Fig. 
IX-10b,  the  flight  Mach  number  vs.  time;  Fig.  IX- 10c,  the  zenith  angle  vs.  time;  and  Fig.  IX-lOd, 
the  PRF  vs.  time.  Note  that  the  PRF  must  increase  froml380Hz  at  liftoff  to  roughly  4800  Hz  when 
“punching-through”  Mach  1;  at  20  seconds  into  the  simulation,  PRF  has  reached  6500  Hz. 

The  acoustic  signature  for  this  case  is  given  in  Fig.  IX-11.  As  indicated,  the  maximum 
vehicle  distance  from  the  source  is  13,000  ft.  or  about  2.5  miles  altitude.  The  doppler-shifted 
frequency  hits  a  low  of  1043  He  and  a  high  of  3070  Hz. 

3)  45®  Liftoff  at  3  G’t  acceleration 

The  initial  conditions  for  this  trajectory  simulation  are  identical  to  the  previous  one  with 
two  exceptions.  First,  the  vehicle  is  launched  at  a  zenith  angle  of  45®,  instead  of  0®.  Second,  the 
observer  is  placed  one  half  mile  downrange  of  the  vehicle  launch  site,  instead  of  coincident. 

The  trajectory /forcing  function  data  for  the  45®  zenith  angle  case  is  presented  in  Fig.  IX- 12, 
and  the  resulting  acoustic  signature  in  Fig.  IX-13.  Note  that  the  temporal  patterns  for  normalized 
amplitude  and  doppler-  shifted  frequency  are  very  different  from  the  previous  case.  The  differences 
become  especially  pronounced  when  heard  on  the  audio  cassette. 

H)  Summary 

The  acoustic  simulations  are  technically  interesting  in  their  own  right  and  give  a  Arm  tac¬ 
tile  understanding  of  what  this  new  class  of  airbreathing  thrusters  will  sound  like.  However,  one 
additional  noteworthy  realization  came  almost  as  a  by*  product  of  the  study.  During  stationary 
hover,  when  simultaneously  pulsed  with  big  energies  and  low  frequencies,  the  ERH  thruster  can  be 
designed  to  be  inaudible  -  in  the  infrasonic  regime.  Alternately,  when  accelerating  at  high  G  levels 
in  the  rotating-line-source  mode,  the  ERH  thruster  can  be  engineered  to  produce  ultrasonic  noise  - 


Figure  DC*  10  Trajectory  Forcing  Function  for  Constant  3  G  Launch  from  0* 
Initial  Zenith  Angle 
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a)  DISTANCE  FROM  OBSERVER  IN  FEET 


d)  VEHICLE  PULSING  FREQUENCY 
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Figure  DC-13  Trajectory/Forcing  Function  for  Coutut  3  G  Launch  from  45s 
Initial  Zenith  Angle. 
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d)  VEHICLE  PULSING  FREQUENCY 


also  beyond  the  range  of  human  hearing.  For  example,  ueume  that  the  1.6  MJ  pulse  energy  shared 
by  the  12  line  source  is  set  off  sequentially;  then  beyond  Mach  1  (with  a  constant  3  Q  launch),  the 
ERH  thruster  PRF  will  exceed  IS  kHt,  which  Is  Inaudible  by  most  of  the  population. 

Through  proper  engineering,  these  ERH  thrusters  hold  good  potential  for  effectively  “silent” 
operation  -  although  sound  pressure  levels  will  still  reach  dangerous  levels  at  close  proximity  to  the 
vehicle. 


CHAPTER X 


CONCLUSIONS 

Preliminary  investigations  of  the  ERH  thruster  concept  have  demonstrated  that  this  ad¬ 
vanced  airbreathing  beam-powered  engine  has  great  potential  for  application  in  future  SSTO  shut- 
tlecraft.  The  level  of  thrust  produced  by  this  engine,  given  a  certain  level  of  input  thermal  power, 
falls  roughly  between  the  non-afterburning  turbojet  and  the  turbofan.  However,  the  acceleration  per¬ 
formance  can  be  Mtoundiny.  Thrust-to-weight  ratios  for  the  ERH  thruster  can  exceed  2$,  whereas 
conventional  turbomachinery  lies  in  the  range  of  6  and  15. 

Although  the  current  study  concludes  that  the  ERH  thruster  has  great  promise,  further 
theoretical  and  experimental  work  is  needed.  Such  proof-of-concept  research  should  investigate  the 
(i)  “refresh*  process  (i.e.,  exhaust/intake  cycles), (ii)  real  gas  effects,  (iii)  wavelength  scaling  effects 
(incL  the  potential  for  microwave  frequencies),  (iv)  LSD  wave  ignition  energy  losses,  (v)  radiation 
heat  transfer  (i.e.,  energy  losses  out  of  the  plasma),  and  (vi)  flight  dynamic  performance  vs.  quasi- 
static  (hover  and  subsonic)  capability. 

Clearly,  there  can  be  no  substitute  for  experimental  verification  of  these  basic  research  issues. 
Several  pulsed  lasers  with  adequate  characteristics  (e.g.,  kiloJoule-level  pulses,  at  roughly  10  psec 
duration)  currently  exist  at  AVCO  Everett  Research  Laboratory,  and  elsewhere.  Such  lasers  would 
be  ideal  for  small  proof-of-concept  experiments. 

In  summary,  the  ERH  thruster  concept  has  survived  this  first  critical  examination  into 
the  propulsive  physics.  However,  further  substantial  increases  m  the  state  of  knowledge  can  only 
come  with  a  well  conceived  program  involving  both  theoretical  and  experimental  elements  -  pushed 
simultaneously. 
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